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1.0  INTRODUCTION 


Long  length,  high  strength,  low  loss  optical  fibers  are  now  being 
routinely  fabricated  through  the  careful  controlling  of  process 
fabrication  parameters.  However,  one  of  the  key  problems  associ¬ 
ated  with  the  application  of  optical  fibers  is  the  degradation  of 
their  high  strength  properties.  The  degradation  in  strength  usu¬ 
ally  occurs  when  the  silica  fibers  are  subjected,  in  the  presence 
of  moisture,  to  a  load  for  an  extended  period  of  time.  This  phe¬ 
nomenon  is  known  as  stress  corrosion  or  static  fatigue.1 

Several  approaches  have  been  taken  in  an  attempt  to  eliminate  or 

reduce  static  fatigue.  One  method  was  to  hermetically  coat  the 

fiber  with  either  metals  or  ceramics.  Although  significant 

improvements  in  fatigue  resistance  were  reported,  the  initial 

fiber  strength  and  the  optical  properties  of  the  coated  fibers 
2  3 

were  degraded. 

Another  promising  technique  involves  cladding  the  fiber  with  a 

4 

compressive  layer  of  glass.  This  surface  compressive  layer  would 
improve  the  fatigue  resistance  of  the  fiber  as  well  as  increase 
its  overall  strength.5 


The  objective  of  the  high  strength,  long  length,  optical  fiber 
contract  was  to  improve  the  durability  of  a  long  length  optical 
fiber  by  increasing  the  surface  compression  of  that  fiber.  The 
surface  compression  will  be  increased  by  modifying  the  interior 
structure  and  composition  of  optical  fibers.  All  fiber  layers 
(with  the  exception  of  the  thin  surface  layer)  will  have  rela¬ 
tively  low  glass  transition  temperatures  and  high  coefficient  of 
thermal  expansion.  The  surface  layer  will  be  fabricated  from 
glass  having  a  higher  glass  transition  temperature  and  lower  ther¬ 
mal  expansion  characteristics.  The  resulting  surface  compression 
of  the  fiber  should  be  on  the  order  of  0.69  GPa  (100  kpsi).  In 
an  effort  to  achieve  these  goals,  the  project  was  divided  into  the 
following  tasks: 

a.  Task  1  -  To  develop  a  computer  model  for  the  three- 
and  four-layer  structure  fibers  and  to  optimize  the 
variables  for  maximum  surface  compression 

b.  Task  2  -  To  develop  a  fabrication  technique  for  pro¬ 
ducing  high  surface  compression  fibers 

c.  Task  3  -  To  evaluate  mechanical  and  optical  charac¬ 
teristics  of  fibers  having  high  surface  compression 

d.  Task  4  -  To  deliver  the  high  surface  compression  pre¬ 
form  and  fiber  samples  to  Naval  Ocean  Systems  Center 


2.0  SUMMARY  OF  WORK  ACCOMPLISHED 
2.1  Theoretical  Computer  Models 

Theoretical  models  were  developed  on  a  computer  for  both  three- 
and  four-layer  optical  fibers.  This  mathematical  model  was  based 
on  the  axisymmetric  stress/strain  equations  for  cylindrical  multi¬ 
layer  viscoelastic  materials.6 

It  is  well  known  that  the  critical  structural  parameter  in  fabri¬ 
cating  a  durable  fiber  is  the  thickness  of  the  compressive  surface 
layer.  However,  the  effects  which  the  material  properties  of  the 
glasses  in  the  inner  layer  have  on  the  final  surface  compressive 
layer  are  less  obvious.  This  is  due  mainly  to  the  interdependency 
of  the  various  thermal  properties  within  each  layer.  The  theo¬ 
retical  model  developed  was  a  preliminary  attempt  at  trying  to 
establish  and  predict  the  various  stress  relationships  occurring 
within  a  fiber.  The  models  resulted  in  8  simultaneous  equations 
for  the  three-layer  fiber  and  11  equations  for  the  four-layer 
fiber.  (See  Appendix  A  for  Theoretical  Stress  Analysis  and  Appen¬ 
dix  B  for  Computer  Model  Programs.) 

The  boundary  conditions  for  the  mathematical  models  were  based  on 
the  actual  thermal  properties  of  the  glass  deposited  in  the  dif¬ 
ferent  layers.  The  thermal  properties  were  measured  at  Harrop 
Laboratories  from  specially  prepared  specimens.  The  measured 
thermal  properties  are  displayed  in  Table  2-1.  To  check  the 

\ 
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Table  2-1.  Thermal  Properties  Used  in  Theoretical  Stress  Analyses 
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accuracy  of  the  theoretical  model,  data  from  the  ongoing  mechan¬ 
ical  evaluation  were  compared  to  the  theoretical  values. 

Tables  2-2  and  2-3  illustrate  how  the  theoretical  model  compares 
with  measured  data.  In  the  case  of  the  three-layer  fibers,  the 
two  values  generally  compare  quite  favorably.  However,  it  should 
be  noted  that  the  theoretical  values  are  generally  less  than  that 
of  the  values  obtained  from  the  measured  static  fatigue  tests. 

This  small  difference  was  attributed  to  the  fact  that  the  values 
obtained  from-static  fatigue  tests  also  account  for  draw  induced 
stresses  as  well  as  surface  compressive  stresses.  The  theoret¬ 
ically  predicted  stresses  on  four-layer  structure  fibers  are  con¬ 
siderably  lower  than  the  measured  values.  This  was  believed  due 
to  the  slight  difference  between  the  thermal  properties  of  the  two 
outer  layers. 

2.2  Fabrication  of  Durable  Fibers 

As  mentioned  in  the  previous  section,  the  surface  compression  of 
fibers  can  be  increased  by  decreasing  the  surface  layer  thickness. 
The  first  approach  used  in  this  project  was  to  reduce  the  thick¬ 
ness  of  the  silica  outer  layer  of  the  present  low  loss  optical 
fiber.  This  was  achieved  by  employing  one  or  more  of  the  follow¬ 
ing  substrate  removal  techniques:  center  grinding,  flame  milling 
and/or  chemical  etching.  During  the  course  of  this  contract, 
these  removal  techniques  were  extensively  investigated  to 


Table  2-2.  Stress  Computation  Data  of  Durable  Fiber  (Three  Layer)  (continued) 


Stress  Computation  Data  of  Durable  Fibers  (Three  Layer)  (continued) 


characterize  their  process  parameters.  The  partial  removal  of 
substrate  material,  using  these  techniques,  produced  three-layer 
structure  fibers. 

The  second  approach  used  in  the  fabrication  of  compressive  layers 

was  the  outside  deposition  of  low  thermal  expansion  glasses  such 

7 

as  pure  silica  and  Ti-doped  silica.  This  was  performed  on  col¬ 
lapsed  preforms  where  the  substrate  material  was  either  all  or 
partially  removed.  Using  this  approach,  either  three-  or  four- 
layer  structured  fibers  could  be  fabricated.  Figure  2-1  illus¬ 
trates  how  the  first  two  approaches  were  used  to  fabricate  durable 
fibers.  The  surface  compressive  layer  thickness  of  the  actual 
fibers  varies  from  0.5  um  to  9  urn. 

Since  both  the  outside  deposition  and  some  of  the  substrate 
removal  techniques  showed  a  tendency  toward  degradation  of  both 
initial  and  fatigue  strength,  a  third  approach  was  investigated. 
This  approach  required  the  internal  deposition  of  a  Ti-doped 
silica  layer.  This  deposition  was  performed  before  the  sequential 
deposition  of  optical  cladding  and  core  glasses  in  the  substrate 
tube.  The  final  fiber  was  drawn  from  the  preform  whose  substrate 
material  was  completely  removed  by  chemical  etching.  This  etching 
would  expose  the  previously  deposited  thin  Ti-doped  silica  layer. 
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Figure  2-1.  Initial  Approaches  of  Durable  Fiber  Fabrication 


This  approach  is  shown  in  Figure  2-2.  Figure  2-3  shows  the  typ¬ 
ical  microphotographs  of  Ti-doped  silica  layer  fiber  cross  sec¬ 
tion.  The  bright  outer  ring  in  the  microscopic  end  photo  of  the 
fiber  is  the  Ti-doped  silica  layer.  This  particular  fiber  showed 
a  high  fatigue  resistance.  A  high  surface  compression  of  0.641- 
0.744  GPa  (93-108  kpsi)  was  deduced  from  the  static  fatigue  tests 

In  addition  to  the  Ti-doped  silica  outer  layer,  the  concentration 
of  higher  thermal  expansion  materials  (such  as  P-doped  boro- 
silicate  glass  and  high  Ge-dcpant  glass)  were  increased  in  the 
inner  layers.  As  a  result  the  highest  fatigue  resistance  was 
observed  in  the  static  fatigue  test  of  higher  dopant  glass  fiber 
(i.e.,  GS  40492). 

A  total  of  186  preform  starts  was  begun  under  this  program.  Of 
these  starts,  164  were  step  index  preforms,  16  were  graded  index 
preforms,  and  6  were  simple  silica  rod-type  preforms.  Of  the  164 
step  index  preforms,  57  were  fabricated  with  an  internal  compres¬ 
sive  layer.  Once  the  preform  was  fabricated,  the  total  preform, 
or  section  of  the  preform,  was  treated.  Of  the  preforms  fabri¬ 
cated,  31  preforms  were  center  ground,  3  were  flame  milled,  72 
were  chemically  etched,  and  48  underwent  outside  deposition. 


12 
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Durable  Preform  Fabrication  of  Internally  Deposited  Compressive  Layer 


A  total  of  93  fibers  was  drawn  from  preforms  or  parts  of  preforms. 
Of  these  drawn  fibers,  95  fibers  were  mechanically  evaluated.  The 
total  number  of  preforms/fibers  fabricated  for  this  program  is 
listed  in  Appendix  C.  Details  of  the  fabrication  processes  are 
described  in  the  Process  Specification,  Item  A006,  of  this  con¬ 
tract  (see  Appendix  E). 

2.3  Evaluation  of  Surface  Compression  Fibers 

Since  the  major  thrust  of  this  program  was  to  increase  the  dur¬ 
ability  of  fiber  through  the  use  of  a  compressive  cladding,  the 
fibers  were  evaluated  mechanically  in  several  different  ways. 
First,  static  fatigue  tests  were  performed  on  the  fiber  at  ambient 
atmosphere  using  the  classical  mandrel  technique.  The  static 
fatigue  resistance  parameter,  N,  was  determined  from  the  slope  of 
log-log  plot  of  breaking  stress  versus  the  time-to-failure  as  dis¬ 
played  in  Figure  2-4.  To  determine  the  surface  compression,  the 
fibers  were  compared  to  plastic  clad  silica  (PCS)  fibers  which 
have  no  surface  compression. 

Dynamic  tensile  strength  of  2  m  gage  length  fibers  was  evaluated 
using  23%  strain  rate  and  the  strength  distribution  was  presented 
as  a  Weibull  plot.  Further  evaluation  such  as  dynamic  fatigue 
tests  on  selected  fibers  was  performed. 
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Figure  2-4.  Static  Fatigue  of  Typical  Durable  Fiber 


Improvement  in  the  durability  of  surface  compression  fibers  was 
determined  quantitatively  by  comparing  the  N  values  and  the  break¬ 
ing  stresses  at  the  specified  time-to-failure  to  the  data  obtained 
from  PCS  fibers.  Fibers  having  a  thin  surface  compressive  layer 
generally  showed  higher  N  values  and  breaking  stress  values.  An 
increase  in  both  the  N  and  breaking  stress  values  was  considered 
to  be  the  result  of  the  compressive  stress  developed  on  the  fiber. 
The  breaking  stress  level  increases  were  compared  to  the  theoret¬ 
ically  predicted  stresses  and  showed  close  agreements  in  most  of 
the  three-layer  structure  fibers.  To  supplement  the  compressive 
stress  level  measurements  a  photoelastic  technique  was  used  to 
determine  the  surface  compressive  stress.4  This  technique  proved 
not  to  be  suitable  for  measuring  stress  in  a  very  thin  compressive 

g 

layer  due  to  its  limited  resolution.  Abraded  strength  tests  were 
also  performed  on  the  available  compression  fibers  to  evaluate  the 

9 

technique  as  an  alternative  testing  method.  These  test  results 
were  not  promising  again  due  to  very  thin  surface  compressive 
layers. 

Finally,  a  pure  bending  stress  test  was  found  to  be  a  simple  and 
reliable  technique  to  indirectly  measure  the  compressive  stress 
level  on  the  fibers.  A  comparison  between  the  surface  compression 
values  determined  using  static  fatigue  test  and  the  values  mea¬ 
sured  by  bending  test  are  shown  in  Table  2-4;  one  of  the  typical 
results  is  shown  in  Figure  2-5. 
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Figure  2-5.  Bending  Test  Results  Compared  With  Plastic  Clad 
Silica  Fibers. 


The  raw  data  from  mechanical  evaluation  test  were  entered  into  a 


computer.  The  data  were  refined  and  the  results  printed  and 
plotted.  An  example  of  this  program  is  displayed  in  Appendix  D. 
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3.0  SUMMARY  OF  FINDINGS 

3.1  Theoretical  Computer  Model 

As  described  previously,  the  main  concept  employed  to  produce  com¬ 
pressive  cladding  depended  upon  the  development  of  thermally 
induced  stresses  within  the  fiber.  These  stresses  were  basically 
dependent  on  the  mismatch  of  the  thermal  properties  (thermal 
expansion  coefficient  and  glass  transition  temperature)  of  the 
various  glasses  used  in  the  fabrication  of  the  fiber.  The  theo¬ 
retical  approach  used  in  this  program  was  essentially  the  same  as 
that  used  by  Krohn  and  Cooper  in  their  two-layer  fiber  study.4 
Similar  assumptions  and  boundary  conditions  were  also  used  in  the 
three-  and  four-layer  theoretical  models. 

It  should  be  noted  that  the  stresses  induced  by  other  mechanisms 
such  as  quenching  stress,  frozen-in  stress,  and  mechanically 
induced  stress  were  excluded  from  the  analysis.8  These  stresses 
are  generally  considered  to  be  insignificant  in  fibers  drawn  under 
typical  drawing  conditions. 

Since  the  induced  axial  stress  distribution  in  multilayer  struc¬ 
tures  must  be  balanced,  the  compressive  surface  layer  thickness  is 
the  critical  structural  parameter  in  fabricating  a  durable  fiber. 
To  maximize  the  surface  compression,  the  layer  thickness  should  be 
minimized.  The  surface  layer  thickness  of  most  fibers  was  more 
than  1  um.  This  was  determined  from  microscopic  observation  of 


cross  sections  of  the  drawn  fibers.  However,  it  should  be  noted 
that  nonuniforra  layer  thickness  was  often  observed  in  the  compres¬ 
sive  layer  of  the  fiber.  The  nonuniformity  was  attributed  to 
variation  in  the  substrate  tube's  wall  thickness  and  unstable  col¬ 
lapse  conditions. 

In  addition  to  the  structural  parameters,  the  material  properties 
of  the  glasses  in  the  fiber  were  also  very  important  in  determin¬ 
ing  the  amount  of  compression  developed  in  the  outer  layer.  The 
two  most  important  material  parameters  in  development  of  a  com¬ 
pressive  layer  were  the  transition  temperature  and  the  thermal 
expansion  coefficient  of  the  various  glasses.  Through  the  use  of  . 
the  theoretical  model,  it  was  determined  that  the  thermal  expan¬ 
sion  coefficient  was  the  most  influential  material  parameter. 

Using  hypothetical  calculations,  surface  compressions  of  up  to 
0.69  GPa  (100  kpsi)  could  easily  be  achieved  in  the  compressive 
layer  by  varying  the  appropriate  thermal  expansion  coefficients. 
The  results  of  these  theoretical  calculations  are  displayed  in 
Figure  3-1.  However,  attempts  to  fabricate  preforms  based  on 
theoretical  calculations  showed  a  high  incidence  of  failure  during 
the  preform  fabrication  processes.  These  failures  were  attributed 
to  the  high  tensile  stresses  in  the  innermost  layer.  In  an  effort 
to  reduce  the  magnitude  of  the  tensile  stress  in  the  innermost 
layer,  the  composition  of  the  core  and  cladding  glasses  was  fur¬ 
ther  modified  in  such  a  way  as  to  minimize  the  difference  in 
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the  thermal  expansion  coefficients.  This,  along  with  structural 
changes  based  on  theoretical  analysis,  resulted  in  the  fabrication 
of  durable  fibers  with  a  compressive  layer  of  0.69  GPa  (100  kpsi). 

The  thermal  property  data  of  optical  glass  systems  used  in  the 
theoretical  analysis  are  shown  in  Table  2-1.  However,  the  actual 
thermal  properties  of  the  glass  systems  could  be  different  from 
the  measured  values  due  to  the  high  cooling  rates  present  during 
the  fiber  forming  process.10  11 

3.1.1  Three-Layer  Structure  Fiber 

Basically  two  types  of  three-layer  fibers  were  fabricated.  The 
two  types  were  similiar  in  most  respects  except  that  type  I  had  a 
plain  silica  compressive  layer  and  type  II  had  a  titanium  doped 
silica  compressive  layer  (see  Table  3-1  for  typical  compositions). 

Surface  compressive  stresses  up  to  0.24  GPa  (35  kpsi)  were  pre¬ 
dicted  from  the  theoretical  models  of  type  I  three-layer  fibers. 

As  shown  in  Table  3-1,  the  theoretical  stresses  of  most  fibers 
were  in  close  agreement  with  the  stresses  determined  experimen¬ 
tally  from  the  static  fatigue  tests.  However,  these  stress  levels 
have  not  been  directly  verified  due  to  the  lack  of  a  suitable 
method  of  stress  measurement.  (This  subject  will  be  discussed  in 
more  detail  in  the  section  on  evaluation  of  compressive  stress 
fibers. ) 
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Three-layer  type  II  fibers  having  a  Ti-doped  silica  surface  layer 
displayed  higher  theoretical  and  measured  compressive  layers  than 
the  type  I.  This  higher  stress  level  is  generally  attributed  to 
the  lower  thermal  expansion  coefficient  of  the  Ti-doped  silica. 
Theoretical  stresses  up  to  0.31  GPa  (45  kpsi)  were  predicted  for 
the  type  II  fiber.  Table  3-1  shows  that  in  every  case  the  theo¬ 
retical  values  were  less  than  the  actual  measured  values.  This 
difference  was  attributed  to  an  inability  to  properly  measure  the 
thermal  properties  of  the  Ti-doped  silica.  The  thermal  properties 
of  Ti-doped  silica  are  strongly  dependent  on  the  fabrication  pro¬ 
cesses.  The  fabrication  processes  present  in  fiber  drawing  could 
not  be  duplicated  in  the  bulk  samples  required  for  thermal  prop¬ 
erty  measurements. 

3.1.2  Four-Layer  Structure  Fiber 

Type  I  of  four-layer  fibers  having  an  extra  silica  layer  formed  by 
outside  deposition  on  the  untreated  Vycor®  substrate  is  listed  in 
Table  2-3.  Stresses  computed  from  four-layer  fiber  models  are 
considerably  lower  than  the  theoretical  stresses  of  three-layer 
models.  Furthermore,  the  theoretically  predicted  stresses  are 
also  less  than  those  of  the  values  obtained  by  static  fatigue 
tests.  The  possible  sources  of  this  discrepancy  will  be  discussed 
in  the  section  on  evaluation  of  the  compressive  stress. 


Four-layer,  type  II  fibers  have  an  extra  thick  silica  layer  which 
was  not  completely  removed  over  an  internally  deposited  thin  Ti- 
doped  silica  layer  (75  um) .  Type  II  fibers  also  have  an  extra 
high  thermal  expansion  cladding  layer  with  high  Ge-dopant  in  order 
to  increase  the  surface  compression.  Therefore,  the  type  II  fiber 
is  actually  a  five-layer  structure  fiber,  but  the  Ti-doped  layer 
is  very  thin  compared  to  the  outer  silica  layer.  The  theoretical 
stresses  computed  from  a  four-layer  model  could  be  inaccurate  due 
to  a  possible  contribution  of  a  Ti-doped  layer. 

3.2  Fabrication  of  Durable  Fibers 
3.2.1  Preform  Fabrication 

All  preforms  were  fabricated  by  the  modified  chemical  vapor  depo¬ 
sition  (MCVD )  technique.  The  optical  borosilicate  cladding  and 
highly  doped  germanosilicate  core  were  deposited  sequentially 
either  in  a  pure  silica  or  a  high  silica  substrate  tube  by  the 
internal  vapor  phase  oxidation  process.  After  collapsing  the  pre¬ 
form,  the  substrate  material  was  either  completely  or  partially 
removed  in  order  to  form  a  thin  surface  compressive  layer  result¬ 
ing  in  a  three-layer  structure  fiber.  In  the  case  of  four-layer 
structure  fiber  fabrication,  the  extra  silica  layer  as  a  compres¬ 
sive  stress  layer  was  externally  deposited  on  the  Vycor®  substrate 
preform.  The  details  of  the  preform  fabrication  processes  are 
described  in  the  process  specification  report  and  the  findings  of 
the  specific  process  will  be  discussed  in  the  following  para- 


3. 2. 1.1  Internal  Deposition 

The  chlorides  of  germanium,  silica,  phosphorous,  boron,  and  tita- 
niiim  are  quantitatively  dispensed  into  the  glass  substrate  tube 
where  they  are  oxidized  and  form  glass  layers.  The  designed 
geometry  is  achieved  after  the  substrate  is  collapsed  into  a  pre¬ 
form.  Various  commercial  substrate  tubes  such  as  WG®-2,  T0®-8, 
SS®-2,  and  Vycor®  were  used  for  the  fabrication  of  durable  fibers. 

The  process  of  internal  deposition  of  these  optical  oxide  glasses 
has  been  well  established  because  these  materials  except  Ti-doped 
silica  are  customarily  used  for  optical  waveguides.  However,  the 
properties  of  each  layer  had  to  be  modified  by  varying  the  amount 
and  type  of  dopant  in  order  to  increase  the  surface  compression  in 
the  durable  fibers. 

Basically,  optical  cladding  is  formed  by  the  addition  of  boron  as 
®2°3'  which  reduces  refractive  index  (RI),  lowers  Tg,  and 
increases  a,  while  the  optical  core  is  formed  by  adding  GeC>2 
and/or  P2O5  which  increase  RI ,  lower  Tg,  and  increase  a.12  The 
practical  fabrication  limit  to  the  amount  of  dopant  should  be  set 
by  allowable  maximum  thermal  mismatch  between  the  substrate  and 
deposited  glasses.  For  example,  if  the  mismatch  is  too  great,  the 
internal  stress  will  cause  the  preform  to  shatter  either  during  or 
immediately  after  preform  collapsing.  This  problem  was  controlled 
to  some  degree  by  use  of  a  ribbon  burner  and  by  changing  inner 


layer  glass  composition  to  modify  the  stress  profile.  Further¬ 
more,  preform  dimension  controls  such  as  circularity  and  unifor¬ 
mity  were  quite  difficult  to  maintain  if  either  deposited  glass  or 
substrate  glass  (such  as  Pyrex®  glass)  is  extremely  soft  or  has  a 
lower  Tg  compared  to  the  silica  glass. 

3.2.2  Substrate  Removal 

Several  techniques  have  been  employed  to  obtain  the  desired  struc¬ 
tural  layer  thickness  in  the  preforms  for  this  program.  These 
techniques  are  described  in  the  following  paragraphs. 

3. 2. 2.1  Center  Grinding 

This  method  consisted  of  mechanical  grinding  of  the  substrate 
layer  for  either  partial  or  complete  removal.  During  the  early 
stages  of  this  program,  mechanical  grinding  was  an  attractive 
method  for  substrate  removal.  Extensive  preliminary  work  was  done 
on  short  preforms  to  establish  optimum  operating  procedures.  The 
centering  of  the  preform  in  the  grinding  lathe  and  the  circularity 
of  the  substrate  tube  were  the  critical  parameters  controlling 
uniform  substrate  removal. 

The  amount  of  glass  to  be  removed  from  the  collapsed  preform  is 
determined  by  measurements  made  on  an  end  photograph  of  a  cane 
drawn  from  the  initial  preform.  Careful  measurements  of  the  cane 
photolayer  ratios  and  of  the  actual  preform  outer  diameter  (od) 


allow  calculations  of  the  depth  to  which  the  preform  must  be 
ground.  After  completion  of  the  grinding,  the  remaining  preform 
was  rinsed  in  diluted  hydrofluoric  acid  and  fire  polished  to  give 
a  smooth  surface  prior  to  any  subsequent  operation. 

In  the  case  of  partial  substrate  removal,  the  preform  is  drawn 
into  a  fiber.  In  the  case  of  complete  removal,  an  outer  layer  of 
a  selected  composition  is  outside  vapor  deposited  over  the  exposed 
cladding  layer  before  the  preform  is  drawn  into  fiber.  (Outside 
deposition  will  be  discussed  in  a  later  section. ) 

Thirty-one  preforms  were  fabricated  using  the  grinding  technique, 
and  various  glass  compositions  were  outside  deposited.  In  all 
cases,  no  appreciable  increase  in  surface  compression  was  observed 
in  drawn  fibers.  In  many  cases,  the  dynamic  strength  test  showed 
a  multimodal  failure  distribution  in  these  fibers.  A  multimodal 
dynamic  strength  distribution  indicates  that  the  static  fatigue 
analysis  is  unreliable. 

It  was  felt  that  the  grinding  process  introduced  flaws  and  impuri¬ 
ties  into  the  preform  surface  which  were  subsequently  transferred 
to  the  fiber  as  it  was  drawn.  In  order  to  investigate  this,  some 
graded  index  preforms  and  some  solid  synthetic  silica  rods  were 
cut  into  several  sections  each,  and  each  section  was  given  differ¬ 
ent  treatment  prior  to  fiber  drawing.  The  results  are  shown  in 
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In  general,  the  treated  silica  rods  when  compared  to  the  untreated 
silica  rods  produced  fibers  with  as  good  or  better  strength  and 
fatigue  characteristics.  However,  when  the  treated  preform  is 
compared  to  the  untreated  preform,  a  degradation  in  the  dynamic 
and  fatigue  properties  was  observed.  In  the  latter  case,  the  pre¬ 
forms  that  showed  degradation  were  treated  by  grinding.  On  the 
basis  of  this  data,  grinding  of  preforms  to  remove  the  substrate 
layer  was  not  chosen  to  be  the  final  technique  for  the  durable 
fiber  preform  treatment. 


3. 2. 2. 2  Flame  Milling 

The  flame  milling  approach  relies  on  a  high  temperature  heat 
source  such  as  a  CO2  laser  or  H2/O2  torches  to  vaporize  the  sub¬ 
strate  material.  In  this  project  flame  milling  was  accomplished 
by  employing  a  slowly  traversing  H2/O2  torch  flame  to  heat  the 
outer  surface  of  the  preform.  Preform  ES  40287  was  treated  in 
this  fashion.  This  preform  was  fabricated  using  a  substrate  tube 
of  Vycor®  glass.  Vycor®  contains  a  slight  quantity  of  sodium,  and 
the  sodium  gives  the  oxyhydrogen  torch  flame  a  distinctive  yellow 
color.  The  absence  of  sodium  yellow  in  the  torch  flames  was  used 
as  an  indicator  of  complete  substrate  removal  for  the  flame  mill¬ 
ing  experiment. 

The  high  silica  substrate  material  has  a  much  higher  melting  tem¬ 
perature  than  the  materials  used  for  core  and  cladding.  Because 
of  the  high  temperature  required  for  flame  milling  Vycor®,  the 
preform  warped  during  milling.  This  warpage  effect  steadily 
increased  as  the  preform  diameter  decreased  during  each  milling 
pass  of  the  torches.  The  warpage  causes  the  preform  to  axially 
run  out  of  the  milling  flame.  This  run  out,  in  turn,  results  in 
incomplete  and/or  nonuniform  removal  of  the  substrate  glass. 

After  several  experiments,  flame  milling  was  abandoned  as  an 
alternative  removal  technique.  However,  this  warpage  problem  can 
be  either  reduced  or  eliminated  by  using  a  preform  produced  with  a 
thick  Vycor®  substrate  tube  (e.g.,  25-mm  precision  Vycor®  tubing). 


3. 2. 2.3  Chemical  Etching 

The  technique  of  etching  was  evaluated  as  a  potential  method  of 
fabricating  a  surface  compressive  layer  on  preforms.  Preforms 
were  etched  with  hydrofluoric  acid  to  form  the  surface  compressive 
layer  by  either  of  two  procedures: 

a.  A  complete  substrate  removal  was  followed  by  an  out¬ 
side  deposition  of  a  thin  layer;  the  deposited  layer 
subsequently  formed  compressive  cladding  on  the  drawn 
fiber. 

b.  A  complete  substrate  removal  exposed  a  previously 
deposited  thin  layer  inside  the  preform.  The  exposed 
layer  subsequently  formed  compressive  cladding  on  the 
drawn  fiber. 

Four  types  of  substrates,  Vycor®,  T0®-8,  WG®-2,  and  SS®-2  were 
evaluated  in  order  to  compare  their  etching  properties.  Prelimi¬ 
nary  evaluations  showed  that  Vycor®  etched  faster  than  the  other 
substrates.  However,  during  the  preform  fabrication  process,  the 
use  of  Vycor®  as  the  substrate  resulted  in  several  bubbles  at  the 
substrate/cladding  interface  and  also  in  the  substrate  itself. 
Preforms  having  bubbles,  inclusions,  or  other  types  of  localized 
inhomogeneity  are  not  suitable  for  the  etching  process.  These 
localized  inhomogeneities  will  sufficiently  change  the  rate  of 
etching.  The  changing  etching  rate  will  result  in  pits  or  bumps 
along  the  length  of  the  preform.  Further  testing  showed  that  TO®- 
8  substrates  did  not  etch  uniformly  along  the  length  of  the  pre¬ 
form  because  of  these  compositional  nonuniformities.  The  WG®  sub¬ 
strates,  although  yielding  preforms  without  pits,  have  an  etching 


rate  which  is  relatively  slow  when  compared  to  Vycor®.  Preforms 
fabricated  with  SS®-2  substrates  generally  did  not  have  any  pits 
and  etched  very  uniformly  along  the  length  of  the  preforms.  Two 
preforms  having  SS®-2  substrates  were  etched  successfully  to  the 
substrate/cladding  interface  without  any  pit  formation  or  damage 
to  the  preform.  It  should  be  pointed  out  that  substrates  used  for 
preform  fabrication  do  not  form  an  integral  part,  of  the  final  pre¬ 
form  having  compressive  cladding.  Therefore,  it  was  desirable  to 
use  an  inexpensive  substrate  which  can  also  etch  uniformly. 
Unfortunately,  the  SS®-2  substrate  tubes  (when  compared  to  other 
substrates)  were  prohibitively  expensive. 

Aside  from  some  of  the  problems  with  the  etching  technique,  the 
technique  was  judged  to  be  the  most  effective  technique  used  in 
this  program.  Most  of  the  final  deliverables  for  this  program 
were  fabricated  using  WG®  substrate  tubes  followed  by  the  etching 
substrate  removal  technique.  There  are  still  some  fundamental 
problems  with  the  chemical  etching  technique.  At  present,  there 
is  no  control  to  assure  uniformity  prior  to  fabrication  in  the 
etched  layer  nor  to  prevent  formation  of  pits.  Possible  solutions 
were  investigated.  One  solution  called  for  the  deposition  of  a 
barrier  layer,  which  would  prevent  or  retard  the  acid  from  etching 
beyond  a  predetermined  depth.  Another  possible  solution  required 
the  use  of  buffered  acids  as  etchants;  this  acid  would  selectively 
etch  glass  of  a  given  composition.  The  data  indicated  the  need 
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for  a  more  extensive  study  of  this  technique  in  order  to  establish 
optimum  etching  conditions.  A  more  complete  investigation  was 
beyond  the  scope  of  this  project. 


During  the  course  of  this  program  the  shattering  of  preforms 
became  a  problem.  It  is  believed  that  occasional  material  defects 
and  nonuniform  layer  thickness,  caused  by  localized  temperature 
distribution  during  preform  collapsing,  resulted  in  the  formation 
of  residual  stress  concentrators  in  the  preform.13  These  stress 
concentrators  were  potential  sources  of  preform  shattering  during 
etching.  When  detected  through  the  use  of  polarized  light  they 
are  eliminated. 


Another  possible  source  of  shattering  could  be  the  joint  between 
the  preform  and  the  T0®-8  holders.  These  joint  sections,  as  well 
as  the  ends  of  the  preforms,  are  covered  with  paraffin  wax  in 
order  to  block  the  possible  effects  of  hydrofluoric  etchant  in 
these  areas.  As  a  result,  the  incidence  of  shattering  problems 
during  the  chemical  etching  process  has  been  considerably  reduced. 


3.2.3  Outside  Deposition 

Outside  deposition  was  achieved  by  attaching  a  hydrogen  oxygen 
deposition  torch  to  the  moving  carriage  of  a  standard  CVD  lathe. 
While  the  collapsed  preform  rotated,  the  deposition  torch  would 
traverse  the  length  of  the  preform  depositing  a  glass  material. 
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By  connecting  the  torch  to  standard  CVD  bubblers,  it  was  possible 
to  deposit  both  pure  silica  and  titanium  doped  silica.  In  order 
to  determine  the  effect  of  outside  deposition  on  the  resulting 
fibers,  a  series  of  experiments  was  performed.  In  this  experi¬ 
ment,  10-mm  SS®-2  rods  were  cut  into  four  or  five  parts  and  drawn 
into  fibers  after  one  of  the  following  treatments: 

a.  Drawn  as  received 

b.  Outside  deposit  only  (Si02) 

c.  Ground  plus  outside  deposit  (Si02> 

d.  Etched  plus  outside  deposit  (Si02) 

e.  Etched  plus  outside  deposit  (Ti02*Si02) 

Table  3-3  displays  the  results  obtained  from  the  experiment.  The 
SS®-2  preform  rods  with  an  asterisk  following  the  number  are  all 
cut  from  one  rod.  The  rod  numbers  which  are  the  same  (except  for 
the  last  letter)  were  all  cut  from  another  SS®-2  rod. 

Although  scattered,  the  data  indicate  that  outside  deposition  of 

silica  (SiC>2)  does  not  seriously  degrade  the  strength  of  SS®-2 

rods.  Similar  results  have  been  reported  by  other  invest iga- 
1  4 

tors.  When  the  outside  deposited  material  is  Ti-doped  silica, 
as  in  the  case  of  ES  404 28D,  there  appears  to  be  a  definite 
i  improvement  in  the  strength  and  durability  of  the  resulting  fiber 


2 


Treated  SS®-2  Rods 


Cut  from  one  rod 


Throughout  the  course  of  the  program  the  outside  deposition  tech¬ 
nique  tended  to  show  an  increase  in  the  durability  of  the  fiber. 
However,  the  data  were  often  very  scattered.  A  possible  source  of 
this  data  scattering  may  be  in  the  outside  deposition  process  it¬ 
self.  During  the  deposition  it  is  believed  that  impurities  pre¬ 
sent  in  the  air  were  incorporated  into  the  deposition  process. 
These  impurities  were  detected  as  small  bright  flashes  of  light  in 
the  torch  flame  during  the  outside  deposition  process.  Although 
some  physical  techniques  were  employed  in  an  attempt  to  reduce  the 
concentration  of  these  solid  impurities,  the  flashes  persisted. 

3.2.4  Fiber  Fabrication 

Final  preforms  having  either  three-layer  or  four-layer  fiber 
structure  were  drawn  to  maximum  fiber  length  after  establishment 
of  the  stable  draw  conditions.  Most  finished  preforms  except 
four-layer  type  I  fiber  preforms  were  short  and  thin  (i.e.,  <10  in 
length  with  5  to  6  mm  outside  diameter  (od))  because  of  substrate 
layer  removal.  The  resulting  fiber  lengths  were  less  than  500  m; 
these  were  not  long  enough  to  permit  complete  evaluation  of 
mechanical  and  optical  performance.  The  optimum  drawing  condition 
of  each  fiber  draw  was  not  easy  to  establish  with  these  short  and 
thin  preforms. 


All  fibers  were  drawn  using  O2/H2  flame  torches  in  order  to  elimi 
nate  the  possible  effect  of  contamination  from  the  furnace  on  the 
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fiber  strength.  The  closed  loop  feedback  diameter  control  was 
employed  to  minimize  fiber  diameter  fluctuations.  The  diameter 
control  of  ±2%  or  ±3%  for  a  127  ym  fiber  diameter  was  achieved  for 
the  most  fiber  draws  with  occasional  higher  fluctuation  up  to  ±5%. 


The  primary  silicone  RTV  (Sylgard®  184)  coating  was  applied  by 
dip-coating  immediately  below  the  fiber  diameter  control  sensing 
head.  The  secondary  coating  of  copolyester  Hytrel®  was  applied  by 
extrusion  on-line  before  the  fiber  was  taken  up  on  the  spool.  The 
final  jacket  diameter  was  20  ±2  mils  (500  ym  ±50  ym). 


3.3  Evaluation  of  Compressive  Stress  Fibers 

The  first  concern  of  this  program  is  the  improvement  of  mechanical 
performance  rather  than  optical  performance  as  addressed  in  the 
objectives  of  this  program.  First,  all  fibers  fabricated  by  the 
developed  techniques  were  subjected  to  mechanical  evaluation  using 
various  evaluation  techniques  as  presented  in  the  original  pro¬ 
posal.  Then,  the  selected  fibers  longer  than  500  m  were  evaluated 
for  their  optical  performance  such  as  attenuation  loss  and  dis¬ 
persion  using  the  standard  measurement  methods.  There  are  a  few 
fibers  whose  lengths  are  longer  than  500  m  after  the  mechanical 
evaluation  due  to  the  limited  volume  of  the  finished  preform. 
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3.3.1  Evaluation  Techniques  -  Mechanical  Performance 

3.3. 1.1  Static  Fatigue  Test 

The  results  of  static  fatigue  tests  performed  on  the  compressive 
stress  fibers  are  listed  in  Table  3-1.  Table  3-3  also  includes 
the  data  for  plastic  clad  silica  (PCS)  fibers  as  the  references. 
The  fatigue  resistance  parameter,  N,  varies  from  19  to  32,  which 
was  determined  from  the  slope.  The  slope  of  the  static  fatigue 
test  curve  was  determined  by  a  regression  analysis  of  more  than  15 
data  points  with  a  higher  correlation  coefficient  (i.e.  >0.95). 
However,  an  increase  in  the  spread  of  fatigue  data  was  occasion¬ 
ally  observed  due  to  large  flaw  size  distribution.  In  this  case, 
the  quantitative  comparison  of  N  values  seemed  to  be  difficult. 
Furthermore,  the  possible  errors  caused  by  the  large  uncertainty 
in  low  values  of  time-to-failure  (i.e.  >10  s)  with  small  mandrel 
size)  were  considered.  The  good  reproducibility  of  a  static 
fatigue  test,  was  confirmed  by  performing  separated  tests  as  shown 
in  Figure  3-2.  Generally,  the  N  value  is  strongly  dependent  on 
glass  surface  composition,  testing  conditions,  crack  growth  char¬ 
acteristics,  and  coating  materials.15 

In  addition  to  the  increase  in  N  value,  the  stress  level  increase 
of  durable  fibers  when  compared  to  PCS  fibers  was  used  as  an  indi¬ 
cator  of  the  relative  fatigue  resistance  improvement.  Therefore, 
the  breaking  stress  level  increase  at  a  specified  time-to-failure 
may  not  represent  the  actual  value  of  surface  compressive  stress 
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developed  on  the  fiber  even  though  there  is  a  correlation  between 

the  compressive  stress  and  failure  time  as  shown  in  the  following 
16  17 

equation: 


log 


N  log  ( 


(3-1) 


where  t  and  tc  are  the  times-to-failure  of  fibers  having  no  com¬ 
pressive  cladding  and  compressive  stress  (ac)  under  the  same 
applied  stress  (aa),  respectively.  The  surface  compression  stress 
data  deduced  from  the  static  fatigue  data  of  typical  three-and 
four-layer  fibers  having  a  silica  surface  layer  showed  relatively 
good  agreement  with  the  theoretically  predicted  data  from  the 
equation  3-1  if  unrealistic  assumption  of  a  constant  N  is  taken 
into  account. 


These  experimental  data  were  also  compared  to  the  predicted  data 
from  the  theoretical  model.  Stresses  associated  with  three-layer 
fiber  models  are  in  close  agreement  with  the  stresses  determined 
experimentally  from  the  static  fatigue  tests  of  fibers  having  a 
silica  compressive  layer.  However,  the  Ti-doped  compressive  layer 
fibers  showed  that  the  theoretical  data  are  relatively  lower  than 
those  determined  from  the  static  fatigue  data  which  showed  the 
high  durability.  Considerable  discrepancies  were  also  noticed  in 
the  four-layer  models,  as  shown  in  Table  3-1.  These  discrepancies 
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could  be  attributed  to  several  possible  sources.  The  exact  values 

of  glass  properties  especially  above  Tg  are  not  known  and  the 

assumption  of  material  properties  such  as  constant  Young's  modulus 

could  not  be  correct  in  the  theoretical  model  as  reported  by  other 
1  6 

investigators.  The  contribution  of  mechanically  induced 
stresses  should  be  considered  in  the  high  silica  substrate  pre¬ 
forms.  For  example,  four-layer  Vycor®  substrate  preform  prepared 
by  outside  deposition  without  the  milling  process  could  give  quite 
different  drawing  conditions  and  also  resulted  in  different  static 
fatigue  characteristics  because  of  unique  surface  flaw  distribu¬ 
tions  associated  with  preform  fiber  fabrication  processes.  Possi¬ 
ble  surface  modification  by  outside  deposition  of  thin  silica 

layer  on  the  Vycor®  substrate  could  result  in  the  crack  geometry 

1  8 

changes  to  influence  the  static  fatigue  characteristics. 

3.3. 1.2  Dynamic  Tensile  Test 

It  was  necessary  to  select  only  fibers  having  a  unimodal  flaw  dis¬ 
tribution  with  higher  Weibull  parameter,  N,  in  order  to  obtain 
reliable  static  fatigue  data.  However,  most  fibers  fabricated  by 
the  durable  fiber  fabrication  techniques  showed  either  multimodal 
flaw  distribution  or  unimodal  distribution  with  a  tail,  which 
indicates  some  degradation  of  fiber  strength.  To  investigate  the 
effects  of  preform  treatments  such  as  substrate  milling  process 
and  outside  deposition  on  the  strength  degradation,  SS®-2  rods  and 
standard  production  graded  index  preforms  were  cut  into  several 


parts  and  drawn  into  fibers  after  the  preform  treatments  for  fab¬ 
rication  of  durable  fibers.  Although  the  data  are  scattered,  pre¬ 
form  treatment  processes  such  as  grinding  and  outside  deposition 
showed  a  tendency  toward  strength  and  fatigue  degradation  as  shown 
in  Table  3-2. 

3.3. 1.3  Dynamic  Fatigue  Test 

To  confirm  the  N  values  determined  from  the  static  fatigue  data, 
dynamic  fatigue  tests  of  selected  fibers  were  also  performed  at 
different  strain  rates.  One  set  of  data  showed  relatively  good 
agreement  between  the  N  values  obtained  from  dynamic  and  static 
methods.  It  was  apparent  that  more  dynamic  fatigue  tests  are 
needed  at  different  strain  rates  to  properly  determine  the  N 
values,  but  most  fiber  lengths  are  too  short  to  complete  the  test. 
Furthermore,  the  median  strength  determined  by  dynamic  tensile 
strength  is  also  influenced  by  the  flaw  distribution. 

3.3. 1.4  Other  Evaluation  Techniques 

Both  direct  and  indirect  stress  measuring  techniques  were  consid¬ 
ered  to  verify  the  compressive  stress  level  on  the  fibers. 
According  to  literature,  photoelastic  technique  and  average 

strength  measurement  are  the  conventional  methods  to  measure  the 

4  19 

surface  compressive  stress  on  glass  materials.  The  photo¬ 
elastic  technique  was  employed  to  evaluate  the  compressive  stress 
of  the  selected  fibers.  The  results  are  not  conclusive.  This 


measurement  technique  has  certain  limitations.  First,  absolute 
values  are  not  determined  directly  because  they  measure  only  dif¬ 
ferences  in  principal  stresses.  Second,  this  technique  places 
restrictions  on  the  size  and  shape  of  the  glass  sample  because  the 
light  path  length  must  be  measured  accurately.  Moreover,  the 
result  is  not  reliable  if  the  glass  contains  optical  inhomegene- 

ities  in  the  form  of  compositional  striae  paralled  to  the 
2  0  21 

surface.  Therefore,  this  measurement  technique  is  not 

totally  suitable  for  measuring  surface  compressive  stress  on 

fibers  having  a  thin  compressive  layer  (i.e.  <5  tun)  due  to  a 

0 

limited  optical  resolution. 

The  average  measured  strength  method  has  met  with  little  favor, 
mainly  because  of  many  parameters  affecting  the  average  strength 
of  optical  fibers.  For  example,  the  2  m  gage  length  tensile  test 
often  showed  high  data  scatter  due  to  association  of  large  flaw 
distribution  resulting  in  difficulties  in  relating  these  measure¬ 
ments  on  surface-compression  fibers  to  the  no-compression  fibers. 

The  additional  evaluation  technique,  in-line  compressive  stress 

monitoring  technique,  as  proposed  in  the  revised  proposal  R77-043, 

was  investigated.  This  technique  was  based  on  the  theory  of  the 

deflection  of  a  free  end  beam  and  did  not  seem  to  be  feasible  to 

detect  small  deflection  of  about  3%  due  to  surface  compression  on 

the  fiber  with  variables  such  as  changes  in  fiber  diameter  coating 

2  2 

thickness  and  draw  speed. 


The  abraded  strength  test  used  by  Catholic  University  was  also 
used  to  evaluate  the  surface  compression.  The  results  were  not 
promising  mainly  due  to  the  difficulty  in  control  of  degree  of 
abrasion  during  mechanical  stripping  of  coating  material  before 
surface  abrasion  and  testing  of  such  very  thin  surface  layer 
fibers. 

As  discussed  before,  the  probability  of  encountering  a  more  severe 
flaw  increases  with  surface  area  (i.e.,  fiber  length)  since  it 
occurs  with  a  statistical  nature.  Therefore,  the  shortest  gage 
length  tests  are  desirable.  These  tests  should  give  a  sufficient 
narrow  stress  distribution  so  that  the  surface  compression  effect 
could  be  clearly  isolated  allowing  the  comparison  of  strengths  of 
fibers  with  compression  to  those  without  compression. 

2  3 

Pure  bending  stress  was  successfully  applied  to  the  minimum  con¬ 
stant  gage  length  using  a  simple  device  (see  Figure  3-1)  and  this 
approach  is  a  very  promising  technique  to  indirectly  measure  the 
compressive  stress  level  on  the  durable  fibers.  This  technique  is 
simple  and  reliable  by  eliminating  the  uncertainty  of  gage  length 
effect  in  the  conventional  tensile  strength  tests  of  optical 
fibers  and  the  need  of  longer  fiber  to  complete  the  test.  The 
results  of  the  pure  bending  stress  tests  on  the  selected  durable 
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fibers  confirmed  the  results  of  experimental  compression  data 
determined  from  the  static  fatigue  test  results  as  shown  in 
Table  2-4 . 


3.4  Optical  Performance 

It  was  assumed  that  the  optical  properties  of  the  durable  fibers 
remain  unchanged  because  the  basic  structure  and  composition  of 
inner  layers  such  as  optical  cladding  and  core  were  kept  constant 
It  is,  however,  to  be  noted  that  optical  performance  of  these 
fibers  could  be  altered  when  fabricated  with  high  residual 
stresses  developed  by  modifying  the  surface  layer  of  the  fiber. 

To  investigate  this  effect,  several  selected  fibers  longer  than 
1  km  were  evaluated;  their  optical  properties  and  the  results  are 
shown  in  Table  3-4.  No  significant  effect  of  fiber  structure  and 
composition  modification  on  the  optical  performance  was  observed 
when  compared  to  the  standard  step  index  fibers. 


Table  3-4.  Optical  Data  of  Selected  Durable  Fibers 


Preform 

No 

Length 

(m) 

NA 

Dispersion 

(ns/km) 

Attenuation 
Loss  @  0.85  tm 

ES  20763 

1387 

0.30 

7.52 

5.57 

ES  20781 

1032 

0.18 

2.80 

2.84 

4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  Theoretical  Model 

Based  on  the  developed  computer  models  of  stress  analysis  for  the 
three-layer  and  four-layer  structure  fibers,  there  are  three 
critical  structural  parameters  in  fabricating  high  surface  com¬ 
pression  fibers.  The  first  one  is  the  thickness  of  a  surface  com¬ 
pressive  layer  which  must  be  as  thin  as  possible  to  maximize  the 
surface  compression.  It  is  believed  that  the  minimum  thickness  of 
approximately  1  um  can  be  sufficient  to  reduce  the  static  fatique 
since  most  submicron-size  cracks  introduced  during  the  fiber  fab¬ 
rication  process  initiate  the  subcritical  crack  growth  resulting 
in  the  time-dependent  strength  degradation. 

The  other  two  parameters  are  the  thermal  expansion  coefficient  and 
glass  transition  temperature  of  each  glass  layer  which  should  be 
optimized  with  regard  to  maximum  thermal  mismatch  between  the 
inner  and  outer  layers  during  cooling.  Therefore,  the  selection 
of  the  glass  compositions  in  each  layer  is  very  important  for 
these  optimum  thermal  properties  of  glass  layers.  Most  impor¬ 
tantly,  the  inner  layer  glass  properties  are  crucial  to  increase 
the  thermal  stress  by  the  cladding  technique.  For  example,  by 
increasing  the  thermal  expansion  coefficient  of  the  core  and  clad¬ 
ding  layers  to  3.5  x  10”®/°C,  a  surface  compression  of  0.69  GPa 
(100  kpsi)  should  result  in  a  1-ym  thick  silica  surface  layer 
fiber.  However,  such  a  high  surface  compressive  stress  structure 


may  have  a  high  probability  of  failure  during  the  preform  fabrica¬ 
tion  process  because  of  the  highly  tensile  stressed  internal  layer 
developed.  Therefore,  the  selection  of  composition  should  be  com¬ 
promised  with  the  feasibility  of  fabrication  as  well  as  the  opti¬ 
cal  performance  characteristics. 

Accurate  theoretical  prediction  of  surface  compression  from  the 
computer  model  requires  actual  glass  thermal/physical  properties 
in  glass  fiber  form  since  the  properties  of  glass  are  governed  by 
thermal  history  of  glass  forming  processes. 

To  minimize  the  contribution  of  other  stress-induced  mechanisms  as 
discussed  in  the  previous  sections,  fiber  fabrication  process 
parameters  should  be  held  constant  during  each  draw.  For  example, 
a  minimum  drawing  tension  should  be  maintained  during  fiber  draw¬ 
ing. 

The  theoretically  predicted  stress  level  should  be  verified  by 
employing  more  than  one  method  including  either  a  direct  or  an 
indirect  stress  measuring  technique. 

4.2  Preform  Fabrication 

The  major  conclusion  of  this  contract  is  that  durable  fibers  with 
a  very  thin  compressive  cladding  (<5  urn  thickness)  of  pure  silica 
or  Ti-doped  silica  can  and  have  been  fabricated.  Data  from  this 
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project  indicate  that  the  center  grinding  technique  may  introduce 
flaws  and  impurities  into  the  preform  surface.  It  is  believed 
that  these  defects  are  subsequently  transferred  to  the  drawn 
fibers,  resulting  in  fibers  with  both  degraded  strength  and  degra¬ 
ded  durability  characteristics.  It  is  concluded  that  this  sub¬ 
strate  removal  method  by  itself  is  not  suitable  for  the  durable 
fiber  program. 

Another  conclusion  of  this  project  is  that  substrate  removal  by 
flame  milling  is  also  not  a  suitable  technique  for  the  fabrication 
of  durable  fibers.  This  conclusion  was  based  on  the  increased 
warpage  observed  during  the  use  of  ‘•he  technique.  Improvement  of 
the  flame  milling  process  was  fully  addressed  in  paragraph  3. 2. 2. 2 
and  should  be  pursued. 

Substrate  removal  by  chemical  etching  is  a  viable  substrate 
removal  technique  using  the  selected  grade  substrates.  The  etch¬ 
ing  may  be  used  in  conjunction  with  either  inside  or  outside 
deposition  of  a  compressive  layer.  The  most  promising  options 
associated  with  this  technique  are  as  follows: 

a.  Partial  or  complete  etching  of  a  substrate  tube 
(followed  by  outside  deposition,  if  necessary) 

b.  Etching  to  a  barrier  layer  (followed  by  outside 
deposition,  if  necessary) 

c.  Etching  of  precision  tubes  (followed  by  outside 
deposition,  if  necessary) 
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In  most  cases,  WG®-2  would  be  the  recommended  substrate  material 
unless  a  barrier  layer  was  employed.  If  a  barrier  layer  were 
used,  Vycor®  or  T0®-8  might  also  be  used  as  a  substrate  material. 
However,  further  investigations  are  required  in  the  use  of  acids 
and  buffer  solutions  before  the  full  potential  of  the  barrier 
layer  method  can  be  evaluated.  It  is  also  recommended  that  either 
silica  or  titanium-doped  silica  be  the  material  used  in  the  com¬ 
pressive  layer. 

The  fabrication  of  a  compressively  cladded  preform  showed  that 
many  preforms  shattered.  They  would  generally  shatter  either  dur¬ 
ing  internal  MCVD  process  or  during  the  substrate  removal  process. 
However,  the  final  surface  compressive  layer  fibers  seemed  to  have 
relatively  fewer  shattering  problems  because  the  thermal  history 
of  fiber  forming  is  quite  different  from  that  of  preform  fabrica¬ 
tion.  Any  residual  stress  in  the  fiber  can  be  removed  by  high 
temperature  annealing  (i.e.,  above  Tg)  the  fiber  ends,  and  the 
possible  difficulty  in  preparing  smooth  fiber  ends  for  splicing 
can  be  overcome. 

A  process  called  vapor  axial  deposition  (VAD)  would  allow  for  the 
fabrication  of  compressively  clad  preforms  without  going  through 
the  two  shattering  prone  steps.  VAD  is  a  process  by  which  various 
glasses  are  deposited  axially  i'n  the  form  of  soot.  The  soot  is  a 
porous  material  which  can  easily  withstand  the  differences  in 


thermal  expansion  of  the  various  glasses.  The  soot  would  later  be 
consolidated  into  a  preform  with  a  compressive  layer  and  drawn 
into  fiber.  This  process  is  recommended  as  a  possible  method  of 
making  fiber  with  compressive  cladding.  It  should  be  possible  to 
fabricate  compressive  cladding  at  or  higher  than  the  0.69  GPa 
(100  kpsi)  achieved  with  relatively  higher  fabrication  yield  by 
using  this  VAD  process.  In  addition  to  the  VAD  approach,  the 
feasibility  of  the  compressive  cladding  technique  combined  with 
the  hermetic  coating  method  should  be  studied  to  increase  the 
fiber  durability.  The  double  crucible  technique  also  should  be 
considered  as  the  alternative  approach  to  fabricate  a  surface  com¬ 
pressive  fiber  with  low  melting  multiple  component  glasses.  The 
desired  thermal  properties  can  be  obtained  easily  from  the  wide 
variety  of  the  multicomponent  glass  compositions  so  that  the  tri¬ 
ple  crucible  method  instead  of  the  double  crucible  technique 
should  be  developed  for  the  multilayer  durable  fiber  fabrication. 
The  double  crucible  process  is  addressed  by  Beales,  et  al.,  Bri¬ 
tish  Post  Office.24 

4.3  Evaluation  of  Compressive  Layer  Stress 

It  has  been  clearly  shown  that  static  fatigue  resistance  of  most 
fibers  having  a  surface  compressive  layer  was  improved.  This  was 
indicated  by  an  increase  in  the  N  values  and  breaking  strengths 
(at  constant  time-to- failure )  when  it  was  compared  to  the  noncom- 
pressive  cladded  optical  fiber.  However,  studies  of  long-term 


durability  improvement  by  the  compressive  cladding  technique 
should  be  extended  so  as  to  gain  a  better  understanding  of  stress 
corrosion  mechanisms  in  various  environments. 

Surface  compressive  stresses  greater  than  100  kpsi  (0.69  GPa)  were 
deduced  from  the  static  fatigue  test  results  of  the  fibers  having 
higher  thermal  expansion  coefficients  in  inner  layers  with  a  thin 
Ti-doped  silica  surface  layer. 

Since  the  photoelastic  technique  is  not  totally  suitable  for  mea¬ 
suring  surface  compressive  stress  on  fibers  having  a  very  thin 
compressive  layer,  an  alternate  measurement  technique  should  be 
employed  to  determine  the  compressive  stress  level  on  the  fiber 
surface.  The  alternative  measurement  technique  recommended  is  the 
bend  testing  described  in  this  report. 

The  bend  testing  is  a  simple  and  reliable  stress  measuring  tech¬ 
nique  which  supported  the  static  fatique  test  results.  Other  dis¬ 
tinctive  methods  of  compressive  stress  measurement  such  as  the 
indentation  technique  should  be  evaluated  to  confirm  the  results 
obtained  by  other  techniques.  < 

To  achieve  better  results,  in-line  proof  testing  should  be 
performed.  The  proof  testing  would  truncate  the  flaw  size 
distribution,  thus  resulting  in  a  more  nearly  unimodal  strength 


distribution.  Additionally,  shorter  gage  lengths  are  recommended 
for  both  tensile  and  bending  tests.  The  shorter  gage  length  would 
also  tend  to  truncate  the  flaw  size  distribution  for  each  speci¬ 
men. 

More  comprehensive  studies  of  internal  stress  effects  on  the  opti¬ 
cal  performance  of  optical  fibers  fabricated  by  the  compressive 
cladding  technique  should  be  conducted  even  though  the  selected 
compressive  stress  fibers  did  not  show  significant  effects  on  the 
optical  characteristics  of  durable  fibers. 

In  summary,  this  program  has  demonstrated  durability  improvement 
of  optical  fiber  waveguides  by  employing  the  compressive  cladding 
technique.  However,  further  refinement  work  is  still  needed  on 
the  developed  fabrication  processes  of  compressive  layer  fibers. 
The  follow-on  program  should  therefore  include  more  extensive 
stress  analysis  such  as  the  basic  fracture  mechanics  (e.g.,  crack 
growth)  and  stress-induced  mechanisms  in  surface  compressive  layer 
fibers  in  addition  to  the  recommended  tasks  already  discussed  in 
this  section. 
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A- 1 . 0  INTRODUCT ION 


Optical  fiber  waveguides  generally  consist  of  several  concentric 
layers  of  different  glass  materials  having  different  mechanical 
and  thermal  properties.  In  the  typical  CVD  optical  fiber,  these 
layers  consist  mainly  of  optical  core  and  cladding  layer  sur¬ 
rounded  by  another  high  silica  substrate  layer.  When  a  compos¬ 
ite  structure  fiber  of  this  type  cools  from  a  temperature  above 
Tg  of  the  glass  to  ambient  temperature,  stresses  will  be  developed 
in  the  axial  direction  as  well  as  tangential  and  radial  directions 
as  shown  in  the  following  figure. 


The  thermal  stress  levels  that  can  be  achieved  in  multilayer 
optical  fiber  structures  can  be  theoretically  predicted  by  a 
mathematical  analysis  based  on  the  theories  of  axisymmetric 
stress  and  strain  equations  for  a  multilayer  cylinder.  The 
following  stress  and  strain  equations  for  an  nth  layer  fiber  can 
be  written  using  Hook's  Law  and  Poisson's  theorem. 


Where  cra  (ea  )  ,  a  (e  )  ,  and  a  (e0  )  are  the  stress  (strain) 
o  o  r  r  z  o 

n  n  n  n  n  n 

components  in  r,  0,  and  z  direction  in  the  nth  layer,  An,  Bn,  and 

C  are  constants.  E  .  y  ,  and  a  are  Young's  modulus,  Poisson's 
n  n  n  n 

ratio,  and  thermal  expansion  coefficient  in  the  nth  layer  material 
respectively,  and  T  is  the  temperature. 

Using  the  following  simplifying  assumptions  and  boundary  conditions 
depending  on  the  actual  thermal  properties  of  glass  layer  materials 
8  simultaneous  stress  (strain)  equations  for  three-layer  fibers 
and  11  equations  for  four-layer  fibers  were  derived  and  solved 
by  a  computer.  (These  equations  are  listed  in  detail  in  the 
following  section.) 

Assumptions : 

Material  properties:  =  E2  =  .  .  .  •  Er,  P2  =  •  •  •  % 

1-2yn  1  ,  a*  5-3  a,K*  ^0.5  K 

E  3  K 

n  n 

where  K  is  bulk  modulus  and  asterisk  designates  glass  properties 
above  the  glass  transition  temperature  (T  ) . 


General  boundary  conditions: 


Equality  of  c at  layer  boundary 

a  =  aQ  =  a  at  or  above  T 
r  0  z  g 

ar  =  o  at  the  free  boundary 

No  net  force  in  z  direction 

Equality  of  e  at  layer  boundary 
z 

In  order  to  maximize  the  surface  compression,  the  surface  layer 
must  contain  the  glass  with  the  lowest  a  and  the  highest  Tg. 
Several  cases  of  cooling  sequence  can  be  considered  to  analyze 
thermal  stresses,  depending  on  the  Tg  of  inner  layer  glasses. 
However  based  on  the  current  optical  fiber  glass  systems,  only 
one  case  for  each  fiber  structure  was  considered  as  shown  in  the 
following  sections. 

The  total  axial  stress  of  the  surface  layer  is  the  sum  of  each 
thermal  stress  constant  (i.e.,  C3  and  C4  in  three-layer  and  four- 
layer  model,  respectively)  computed  by  solving  the  simultaneous 
equations.  The  examples  of  numerical  computation  are  shown 
in  Appendix  B. 


Three-Layer  Model 


Case: 


T  >  T  >  T 
*3  *1  *2 


A.  Cooling  from  T  to  T  ,  AT  *  T  -  T 


The  stresses  in  the  first  and  second  layers  must  be  hydrostatic 

%  *  °Z1  *  a9j/  C1  -  A1 

C1 

£r1  “  ez1  “  e01"  3Kp"  +  al  T1 


ar2  =  az2  "  a02'  C2  *  A2'  B2  "  0 


£r2  "  ez2  e02  3KX 


Boundary  Conditions: 


1.  a  =0  at  r  *  c, 
r3 


2.  a  *0  at  r  =  b, 
r3  r2 


*  +  a2  ^2 


K;*  =  0.5K 


A3  ‘  ?  =  ° 

C2  "  ^  B3 

c  b 


3.  Ho  net  forces  in  the  axial  direction 
a2 3 4Cx  -  (b2-a2)  C2  +  (c2-b2)  C3  =  0 


4.  e  -  e_  at  r  *  b 
*2  z3 ' 


C2  C3~2yA3 


AcuAT 


Act,  =  a--3a„ 


3K*  “  3!?^  Aa2AT  Aa2AT  -  3(a2~a1) 


B.  Cooling  from  T  to  T  ,  AT  *  T  -  T 

gl  g2  gl  g2 

The  stresses  in  the  second  layer  are  hydrostatic 


Boundary  Conditions: 


1. 

Q 

U 

O 

at  r  * 

c. 

r3 

2. 

a.  -  a 

at  r  = 

a 

rl  r2 

3. 

A2C1  +  (b2 

-a2)  C2 

+  (c 

4. 

ez  =  ez 

Z1  Z2 

at  r  = 

a 

r- '  r-  c2  - 


b2)  C3  =  0 

Aa2*  *  3a2~ai 

C9  -  (^y  -  ^»r)  B,  =  0 


5.  0^ 


a 


at  r  =  b 


at  r  =  b 


6. 


=  e, 


(i-p)  a2  +  B' 

z  bz 


7. 

8. 


e  at  r  =  a 
z2 


e  at  r  =  b 
z3 


E2  =  Aa^  =  a3  ~  a2 

yC2  -  (1-y)  A3  -  -^T  B2  +  yC3  = 

b 

EAa3AT  -  -  -  -  -  (?) 

C1  “  B2  "  C2  =  EAa2  AT - © 

cl 

C2  -  2yA2-C3  +  3yA3  =  EAa-jAT  -  (?) 


Four-Layer  Model 


Case: 


P  >  T  >  T  >  T 
g4  g3  gl  ? 


A.  Cooling  from  T  to  T  ,  AT  =  T  -  T 

g4  g3  g4  g3 

All  layers  except  outer  layer  (4th  layer)  are  hydrostatic 


a  =  a  *  a 
rl  z2  91 


C1  =  A1 


Erl  S  =  £01  =  +  aiTl 


rr  =  az  =  °fl 
r2  z2  9 2 


C2  A2 


b2  =  0 


:r2  *  ez2  *  ee2  "  3K*2  +  a2T2 


°r  "  °z  *  a9 

r3  z3  °3 


C3  "  A3 


b3  -  0 


A- 9 


D.  Cooling  from  T  to  T  ,  AT  =  T  -  T 

g?  o'  g0  o 


All  layers  are  solid  state 


B , 


1.  a  =0  at  r  =  d 


A‘  ’  7  *  0 


B. 


2.  o_  =  a  at  r  =  c 


A3 - 2  +  B4  ( 

c 


3.  a  =  a  at  r  =  b 
r3  r2 


A2  "  A3  + 


B. 


4.  a  =  a  at  r  =  a 
r2  rl 


A1  +  A2  2 
a 


5.  No  net  force  in  Z-direction 


a2C1  +  (b2-a2)  C2  +  (c2-b2)  C3  +  (d2-c2)  C 


6.  efi  =  e  at  r  =  a 
91  02 


(l— y)  ax  -  (l-y)  a2  -  b2  -  uc±  + 


EAa2AT2  where  Aa2  =  (a2»a^) 


7.  e„  =  e0  at  r  =  b 
9  2  9  3 


(l-y)  A2  -  (l-y)  A3  +  B3  -  yC2  +  y^C 

b » 


GAcx-AT  where  AV„  =  V,  -  cu 


cr 


8 


*  ^  A  ^  a 

9 3  94  at  r  =  c 


(1-U)  A,  -  (1-y) 

c 


3  h/  «4  ■  — *i —  -  yC^  ~  yC4  - 


EAa.AT  where  Aa.  =  a*  -  a, 

4  4  4  3 

9.  e _  =  e„  =  at  r  =  a 

Z1  z2 

Bj  +  -  C2  =  EAa2AT  where  Ao^  =  a2  ”  ai 

2L 


.10 .  e  =  e  at  r  =  b 
z2  z3 


-2y  A2  +  2y  A^  +  C2  ~  =  EAa^AT  where  Aa^  = 


-  a. 


11.  e„  =  £  at  r  =  c 
z3  z4 


-2y  A3  +  3yA4  +  C3  “  C4  “  EAa4AT  where  Aa4  - 
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APPENDIX  B 

COMPUTER  MODEL  PROGRAMS 


B-l 


B-1.0  THREE  LAYER  SOLUTION 

The  following  is  an  illustration  of  a  typical  solution  for  a 
theoretical  three-layer  compressive  layer  calculation.  The 
amount  of  compression  in  the  outer  layer  of  this  example  is 
30.71  kpsi  (0.21  GPa) .  This  value  is  arrived  at  by  adding 
the  three  eighth-solution  values  together  (i.e.,  0.196302E+05  + 
0 . 109455E+05  +  O'.  129586E+05  =  0 . 307502E+05)  .  This  solution 
was  performed  on  an  IBM  370/138  computer. 


B-2 


PILE:  TERM2 


FILE 


A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


const456  consteva  datal  ml  cntl  m2  cntl  b3  cntl  t 


EXECUTION  BEGINS... 

ENTEB  A, 3,C, ALPHA  1 , ALPH A2, ALPH A3 , TG 1 , TG 2 , TG3 , TO , E# K,  MUE,ISIZE 

EXECUTION  BEGINS _ 

ENTER  A, 3,C, ALPH A1 , ALPHA 2, ALPH A3 ,TG 1 , TG2 , TG3 , TO , E, K,  MUE,ISIZE 

EXECUTION  BEGINS... 

ENTER  A,3,C, ALPHA1,ALPHA2, ALPHA3,TG1,TG2,TG3,T0,E,K,  MUE,ISIZE 

R; 


soln  ml  cntl  t. 


EXECUTION  BEGINS... 

SOLUTION  VALUES 

1  0.0 

2  0.0 

3  0. 196241E+05 

4  0.0 

5  0. 120689E+00 

6  -0. 204689E*04 

7  -0. 204689S+04 

8  0.  196302E4-05 
EXECUTION  TERMINATED 

R; 


soln  m2  cntl  t 


EXECUTION  BEGINS... 

SOLUTION  VALUES 

1  -0. 169335E+04 

2  -0.  169335E+04 

3  0.  162346E+05 

4  0.0 

5  0. 998438E-01 

6  0. 847809E*03 

7  -0. 169335E+04 

8  0. 109455E+05 
EXECUTION  TERMINATED 

R; 


soln  m3  cntl  t 


EXECUTION  BEGINS. . . 

SOLUTION  VALUES 

1  -0. 338300E+04 

2  0.754102E+02 

3  0. 647890E+04 

4  0.41849  2E-02 

5  0.398457E-01 

6  -0.676553E+04 

7  0.151375E+03 


The  following  files  were  used  in  the  solution  of  the  previously 
mentioned  theoretical  three-layer  compressive  layer  calculation. 
Files  CONST456  EXEC  and  SOLN  EXEC  are  executives  used  to  call 
other  files  and  control  the  solution. 

Files  C4 ,  C5 ,  C6 ,  SOLN,  MATIN,  MXOUT,  SIMQ  and  LOC  are  FORTRAN 
source  listings  use  in  the  execution  of  the  solution. 

File  CONSTEVA  DATA  1  is  the  file  where  the  particular  data  used 
in  the  solution  are  stored. 
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SCONTROL  OFF  30MSG 

GLOBAL  TXTLIB  FORTMOD2  MOD2EEH  HOD2HEEH  CMSLIB 
FILE  FT05F00 1  DISK  SI  S2 
FILE  FT06F001  TERM 

FILE  FT07F00 1  DISK  S3  S4  (LRECL  80  3LOCK  800  RECFM  F3 

LOAD  C4 

START 

FILE  FT05F001  DISK  SI  S2 
FILE  FT06F001  TERM 

FILE  FT07F001  DISK  S5  S6  (LRECL  80  BLOCK  800  RECFM  FB 

LOAD  C5 

START 

FILE  FT05F001  DISK  SI  62 
FILE  PT06FQ0 1  TERM 

FILE  FT07F001  DISK  S7  S8  (LRECL  80  BLOCK  800  RECFfl  FB 

LOAD  C6 

START 
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REAL  M,N,MUE,R 
DIMENSION  N  { 1 1 ) 

DIMENSION  M  (11, 1 1) 

NRITE  (6,2) 

2  FORMAT ( *  ENTER  A, B, C, ALPHA  1 , ALPHA2 , ALPHA 3, TGI , TG2,TG3, TO , E, K 
CMUE, ISIZE  * ) 

READ  (5,*)  A,3,C, ALPHA  1, ALPH A2, ALPHA3,TG1,TG2,TG3, TO,E 

READ  (5,*)  K,MUE, ISIZE 
H  ( 1,  1)  =1.0 
M  (1 ,2)  =0.0 

M  ( 1 , 3)  =0.5 
M (1,4) =0.0 
M(  1,5)  =0.0 
M(1 ,6) =0.0 
M  ( 1 , 7)  =0.0 
M(1,8)=0.0 
M  (2, 1)  =1. 0 
M  (2 , 2)  =1.0 
M  (2, 3)  =0.0 
M (2 ,4) =0.0 
M  ( 2 , 5)  =0.0 
M  (2 , 6)  =0.0 
M  ( 2 , 7)  =0.0 
M  (2,8)  =0.0 
M  (3,  1)  =0.0 
N  (3, 2)  =0.0 
M(3,3)  =  1.0 
M (3 ,4) =0.0 
M (3,5) =-1/C**2 
M  (3,6)  =0.0 
M  (3,7)  =0.0 
M  (3 , 8)  =0.0 
M  (4,  1)  =1.0 
M  (4,2)  =1.0 
M  ( 4 , 3)  =0.0 
M(4,4)  =1.0 
M  (4, 5)  =0.0 
M  (4,6)  =0.0 
M  (4, 7)  =0.0 
H  (4, 8)  =0.0 
M  (5, 1)  =0.  0 
N  (5, 2)  =0.0 
M  (5, 3)  =0.  0 
M (5, 4) =0.0 

M  (5, 5)  =-  (1/C**2-1/B**2) 

M (5,6) =0. 0 
M(5,7)  =1.0 
M  (5, 8)  =0.0 
M  (6, 1)  =0.0 
M  (6 , 2)  =0.0 
M  (6,3)  =0.0 
M  (6, 4)  =0.0 
M (6,5) =0.0 
M (6 , 6) =A**2 
M  (6 ,7)  =  (B**2-A**2) 
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N  (6,8)  =  (C**2-B**2) 

M  (7, 1)  =0.0 
H  (7, 2)  =0. 0 
M(7,3)  =2*MUE/E 
N (7,4) =0.0 
M  (7, 5)  =0.0 
M  (7,6)  =0.0 
H(7,7)  =  (1/(1.5*K)) 

M(7,8)=-1/E 
M  (8,  1)  =0.0 
M  (8, 2)  =0.0 
M  (8 , 3)  =0.0 
M (8,4)  =0.  0 
M (8, 5) =0.0 
M  (8,6) =1.0 
M  ( 8 , 7}  —1.0 
8(8,8)  =0.0 
1=9 

N  (  1 ) =0 . 0 
N  (2) =0.0 
N (3) =0. 0 
N (4) =0.0 
N  (5)  =0.0 
N (6) =0.0 

N  (7)  =  (  (TG3-TG1)  *  (ALPHA 3-3* ALPHA 2)  ) 
N  (8)  =0.0 

WRITE (7, 10)  ISIZE, ISIZE 

10  FORMAT  ( ’  9999*, 214) 

DO  200  J=1, ISIZE 

WRITE  (7,11)  (M  (J,  K  1)  ,  K  1=  1 , 7) 

11  FORMAT (7E10. 4) 

WRITE  (7,12)  (M  (J,K1)  ,K1=8,  11) 

12  FORMAT  (4E10. 4) 

200  CONTINUE 

WRITE (7, 13)  I 

13  FORMAT  (II) 

WRITE  (7,11)  (N  (Kl)  , K 1=  1 , 7) 

WRITE  (7,  12)  (N  (Kl)  ,K1=8,  11) 

WRITE  (7, 14) 

14  FORMAT  (80  (•  •)  ) 
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REAL  M,N,HUE,K 
DIMENSION  N  ( 1 1) 

DIMENSION  M(11,11) 
iRITE  (6,2) 

FORMAT  ( '  ENTER  A , B, C, ALPHA  1 , ALPH A2 , ALPHA3, TG 1 ,TG2 ,TG3 , TO , E , K, 
CMUE,ISIZE' ) 

READ  (5,*)  A,B,C, ALPHA1, ALPHA2, ALPH A3, TG 1 , TG2 , TG3 ,T0 , E 
READ  (5,*)  K,MUE,ISIZ2 
FORMAT  (10E5.5) 

8(1,1)  =0.0 
0(1,2)  =0-0 
M  (  1 , 3)  =1.0 
8(1,4)  =0.0 
M  ( 1 , 5)  =-  1/C+*2 
8(1,6)  =0,0 
M  ( 1 , 7)  =0.0 
M  ( 1 , 8)  =0.0 
8(2,1)  =0.0 
M  (2, 2)  =1.0 
M(2,3)  =0.0 
8(2,4)  =0.0 
M  (2,5)  =0.0 
M  (2,6)  =0.0 
M  ( 2 ,7)  —1.0 
M  (2, 8)  =0.0 
M  (3 , 1)  =  1. 0 
M  (3 , 2)  =0.0 
M  (3 ,3)  =0.  0 
M  (3 , 4)  =0.0 
M  (3 , 5)  =0.0 
M(3, 6)  =0.0 
8  (3,7)  =-1.0 
M  (3 , 8)  =0.0 
M  (4 , 1)  =0.  0 
M  ( 4 , 2)  =0.0 
M  (4,3)  =0.0 
M  (4,4)  =1.0 
M  (4 , 5)  =0.0 
M  (4,6)  =0.0 
M (4,7) =0.0 
M  (4,8)  =0.0 
3(5,1)  =0.  0 
M  (5, 2)  =0.0 
M  (5, 3)  =0.  0 
M (5,4) =0.0 
M(5,5)=0.0 
M  (5,6) =  A**2 
M  (5 , 7)  =  (B**2-A**2) 

H  (5, 8)  =  (C**2-B**2) 

8(6,1)  =-  (2+MUE/E) 

M  (6 , 2)  =0.0 
8(6,3)  =0.0 
M  (6,4)  =0.0 
8(6,5)  =0.  0 
8(6,6)  =♦  (1/E) 
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M(6,7)=-(1/(1.5*K)) 

M(6,8)  =0.0 
M(7,1)=0.0 
M  (7,2)  =0.0 
t!  (7,3)  =0.0 
N(7,4)  =0.0 

M  (7,5)  =-  (1/C** 2-  1/3**2) 

M  (7 , 6)  =0.  0 
M<7,7)  =1.0 
N (7, 8)  =0.0 
H  (8,  1)  =0.0 
M  (8, 2)  =0.0 
M  (8,3)  =  (2*MUE/E) 

H  (8,4)  =0.0 
M  (8,5)  =0.0 
M (8, 6)  =0.0 
M(8,7)  =  (1/(1.5*K)) 

11(8,8)  =-1/2 

1=9 

N  ( 1)  =0.0 
N (2) =0. 0 
N (3) =0. 0 
N  (4) =0. 0 
N (5) =0 . 0 

N  (6)  =  (  (3*ALPHA2-ALPHA  1 )  *  (TG1-TG2)  ) 
N (7) =0.0 

N  (8)  =  (  (AL?HA3-3*ALPHA2)  *(TG1-TG2)  ) 
WRITE (7 , 10)  ISIZE, ISIZE 

10  FORMAT  ( ’  9999', 214) 

DO  200  J=1, ISIZE 

WRITE  (7,11)  (M  (J  ,  K  1)  ,  K 1=  1 , 7) 

11  FORMAT  (7E10. 4) 

WRITE  (7,12)  (M  (J,K1)  ,S1=8,  11) 

12  FORMAT  (4E1  0.4) 

200  CONTINUE 

WRITE (7, 13)  I 

13  FORMAT  (II) 

WRITE  (7,11)  (N  (K1)  ,K  1=  1 ,7) 

WRITE (7,12)  (N ( K 1 ) ,K 1  =  8,  11) 

WRITE (7,14) 

14  FORMAT  (80  (*  «)  ) 
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REAL  M,M,MU3,K 
DIMENSION  N ( 1 1) 

DIMENSION  M  ( 11 , 1 1) 

WR IT3  (6,2) 

2  FORMAT  (‘  ENTER  A, B, C, ALPHA  1 , ALPHA2 , ALPHA3, TG 1 , TG2, TG3 , TO , E, K 
CMOS, ISISE* ) 

READ  (5,*)  A,B,C,ALPHA1,ALPE1A2,  ALPE1A3,  TG 1 ,  TG2  ,  TG3  ,T0,  E 
READ  (5,*)  K , MOE ,  ISIZE 
C5  FORMAT  (10ES.5) 

M(1,  1)  =0.0 

M  ( 1 , 2)  =0.0 

H  ( 1,3)  =1.0 

M  ( 1 , 4)  =0.0 

M(1,5)=-1/C**2 

S  (1,6)  =0.0 

H  f  1,7)  =0.0 

M(  1,8)  =0.0 

M  (2,  1)  =0.  0 

M  (2, 2)  =  1.0 

M  (2,3)  =0.0 

M  ( 2 , 4)  =-  1/P**2 

M  (2,5)  = (1/B**2-1/C**2) 

M  (2 , 6)  =0.0 
M  (2,7)  =0.0 
M  ( 2 , 8)  =0.0 
M  ( 3,  1)  =1.0 
M  ( 3 , 2)  =-1.0 
M  (3, 3)  =0.0 
M (3,4) =1/A**2 
M  (3, 5)  =0.0 
M  ( 3 , 6)  =0.0 
M (3,7) =0.0 
M  (3,8)  =0.0 
M  (4,  1)  =0.0 
M  ( 4 , 2)  =0.0 
M  ( 4 , 3)  =0.0 
M  ( 4 , 4)  =0.0 
M  (4 , 5)  =0.0 
M (4 ,6) = A**2 
M  (4 , 7)  =  (B**2-A**2) 

M  (4, 9)  =  (C**2-B**2) 

M  (5,  1)  =  ( 1-MOE) 

M (5, 2) =—  ( 1-MUE) 

M  (5,3)  =9.0 

M  (5,4)  =-  (1  ♦HUE)  /A* *2 
M  (5,5)  =0.0 
M (5,6) =-MUE 
M (5, 7) =  MOS 
M  (5,8)  =0.0 
M(6,1)=0.0 
M  (6,2)  =  1-MUE 
M  (6 , 3)  =-  (1-MUE) 

M  (6 , 4)  =  (1  +  MUE)  /3**2 
M  (6,5)  =-  ( 1 +MUE)  /B**2 

M (6 , 6) =0. 0  ^ 
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M  (6,7)  =-MUE 

M  (6,8)  =  «1UE 

M  (7,1)  =0.0 

M  (7 , 2)  =0.0 

a(7,3)  =0.0 

M(7,4)  =  +  (2*MUE/A**2) 

M(7,5)=0.0 
a(7,6) =1.0 
M  (7, 7)  =-1.0 
a(7,8)  =0.0 
a  (8,  i)  =o.o 
a  (8, 2)  =-2*aoE 
H  (8, 3)  =2*MUE 
H  ( 8 , 4)  =0.0 
a(8,5)  =0.0 
a(8,6)  =0.0 

M  (  8, 7)  =1.0 
a  (8,8)  =-i.o 
1=9 

N  ( 1)  =0. 0 
N  ( 2)  =0.0 
N  (3)  =0.0 
N (4) =0.0 

H  (5)  =E*  ( 4L  PH A2— ALPHA  1)  *  (TG2-T0) 
N  (6)  =S*  ( ALPHA3-ALPHA2)  *(TG2-T0) 
N  (7)  =E*  ( ALPHA2- ALPHA1)  *  (TG2-T0) 
N  (8)  =  2*  (ALPHA3-ALPHA2)  *  (TG2-T0) 
WRITE  (7, 10)  ISI2E,ISISE 

10  FORMAT  (•  9999', 214) 

DO  200  J= 1 , ISIZE 

HR ITE  (7,11)  (H  (J,K1)  ,K1=1,7) 

11  FORMAT  (7E10. 4) 

WRITE  (7,12)  (M  (J,  K 1)  ,  K  1=8,  11) 

12  FORMAT  (4E10. 4) 

200  CONTINUE 

WRITE (7,13)  I 

13  FORMAT  (II) 

WRITS  (7,11)  (N  ( K 1 )  ,X  1=  1 ,7) 

WRITE  (7,  12)  (N  (K  1)  ,X  1=8,  1 1) 
WRITE  (7, 14) 

14  FORMAT  (80  (•  ')  ) 
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SCONTROL  OFF 

GLOBAL  TXTLI3  FORT  MO  D2  MOD2EEII  MOD2NEFH  CMSLIB 
SIP  S3  =  ?  SGOTO  -A 
SIF  S3  =  T  SGOTO  -B 

-C 

PILEDEP  FT05F00 1  DISK  SI  S2 

LOAD  SOLN 

START 

SEXIT 

-A 

PILEDEP  FT06P001  PRINTER 
SGOTO  -C 
-B 

PILEDEP  FT06F001  TERM 
SGOTO  -C 
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C  PROGRAM  NAME  —  SOLN 

DIMENSION  A  (12  1)  ,B(11) 

10  FORMAT  (1H1 , 34HSOLUTION  OF  SIMULTANEOUS  EQUATIONS) 

11  FORMAT  ( 1  HO, 440 DIME NS 10 NED  AREA  TOO  SMALL  FOR  INPUT  MATRIX, 14) 

12  FORMAT (inO,20HEXECUTION  TERMINATED) 

13  FORMAT  (1HO,47HROW  AND  COLUMN  DIMENSIONS  NOT  EQUAL  FOR  MATRIX, 14) 

14  FORMAT  (1HO,42HINCORRECT  NUMBER  OF  DATA  CARDS  FOR  MATRIX, 14) 

15  FORM  AT  ( 1H0 , 1 8HG0  ON  TO  NEXT  CASE) 

16  FORMAT  (1H0,38HSTRUCTURE  CODE  IS  NOT  ZERO  FOR  MATRIX, 14) 

17  FORMAT  (1U1, 17H0RIGINAL  3  VECTOR,////) 

18  FORMAT  ( 1  HI , 15HS0LUTI0N  VALUES,////) 

19  FORMAT  (1H0, 18IIMATRIX  IS  SINGULA?) 

20  FORMAT  (7F10.0) 

21  FORMAT  (13,  10X,E16- 6) 

22  FORMAT  (1  MO,  1  1HEND  OF  CASE) 

C  WRITE  (6,10) 

25  CALL  M ATIN (ICOD, A, 121 , M,N, MS,IER) 

IF  (N)  30,95,30 
30  I?  (IE3-1)  45,35,40 

35  WRITE (6 , 11)  ICOD 

GO  TO  90 

40  WRITE  (6, 14)  ICOD 

GO  TO  95 

45  IF(N-M)  50,55,50 

50  WRITE  (6, 13)  ICOD 

GO  TO  90 

55  IF  (MS)  60,65,60 

60  WHITE(6, 16) ICOD 

GO  TO  90 

65  CALL  MXOUT  (ICOD, A, N, M, MS, 60,  120,2) 

C 

C  ADDED  END=98  TO  STOP  READING  PAST  END  OF  FILS  CONDITION 

C 

C  READ  (5  ,*,END=98)  (B(I)  ,1=1, N) 

READ(5,20,END=98)  (B(I)  ,1=1, N) 

C  WRITE  (6, 17) 

C  DO  70  1=1, N 

C  WR ITS (6,21)  I, B  (I) 

C  70  CONTINUE 

CALL  SIMQ(A,B, N,XS) 

IF(KS-I)  80,75,80 
75  WRITE (6,19) 

WRITE(6, 15) 

GO  TO  25 

80  WRITE (6,18) 

DO  85  1=1, N 
WRITE  (6,21)  1,0(1) 

85  CONTINUE 

C  WRITE  (6, 22) 

GO  TO  25 
C 

C  ADDED  SMD=98  TO  STOP  READING  PAST  END  OF  FILE  ERROR  CONDITION 

C 

CO  RE  AD  (5  ,  *  ,  END=98)  (B(I),I=1,N) 

90  READ  (5, 20,  END=98)  (B(I)  ,1=1, N) 
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SUBROUTINE  MXOUT  (ICODE, A, N, M, MS, LINS, I POS , IS?) 

DIMENSION  A  ( 1)  ,  B  (3) 

1  FORMAT (1H1,5X,  7HMATRIX  ,I5,6X,I3,5H  ROWS, 6X, 13, 8H  COLUMNS, 
18X, 133STCRAGE  MODE  , I 1 , ST, 5HPAGE  ,12,/) 

2  FORMAT (12X,8HCOLUMH  ,7(3X,I3,1 OX)) 

3  FORMAT  (II!  ) 

4  FO S M AT ( 1 H  ,7X,4HROW  , I  3, 7 ( El  6. 6 ) ) 

5  FORMAT  (130, 7X,4HROW  , 1 3, 7 ( E 1 6. 6)  ) 

J=1 

NEND=I POS/ 16-1 
LEND= (LINS/ISP) -2 
IPAGE^I 
10  LSTRT= 1 

C20  WRITE(6, 1) ICODE, N,M, NS, IPAGE 
20  CONTINUE 

JMT= J+  MEND- 1 
I?AGE=I?AGE+1 

31  IF (JNT-M) 33,33,32 

32  JHT=M 

33  CONTINUE 

C  B3ITE (6,2)  (JCUR, JCUR=J, JNT) 

IF  (ISP-1)  35,35,40 
C35  WRITE  (6, 3) 

35  CONTINUE 

40  LTEWD=LST3T+L2ND- 1 

DO  80  L=LSTBT,LTEND 
DO  55  K= 1 , NEND 
KK  =  K 

JT  =  J  ♦K- 1 

CALL  LOC(L,JT,IJNT,N,M,MS) 

3(K)=0.0 

IF  (IJNT)  50,50,45 
45  B  (K)  =  A  (I  JNT) 

50  CONTINUE 

IF(JT-M)  55,60,60 
55  CONTINUE 

C60  IF  (ISP- T)  65,65,70 

60  IF  (ISP-1) 75,75,75 

65  WRITE  (6,4)  L,  (B  (JW)  ,  JW=  1 , XK) 

GO  TO  75 

70  WRITE  (6,5)  L,  (B(JW)  ,JW=1,KK) 

75  IP  (N-L)  85,85,80 

80  CONTINUE 

LSTRT=LSTFT*LEND 
GO  TO  20 

85  IF  (JT-M)  90,95,95 

90  J=JT+1 

GO  TO  10 
95  RETURN 
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SUBROUTINE  BATIN (ICODE,  A ,ISIZS , IROW, ICCL, IS, IER) 
DIBENSION  A  ( 1) 

DIMENSION  CARD  (8) 

1  FORMAT  (7P10.0) 

2  FORMAT  (16,214,12) 

IDC=7 

IER=0 

READ (  5, 2,END=40)  I CODE, I ROM, I COL, IS 
CALL  LOC (I ROM, ICOL, IC NT, IROW, ICOL, IS) 

IF  (ISIZE-ICNT)  6,7,7 

6  IER= 1 

7  IF  (ICNT)  38,38,8 

8  ICOLT= ICOL 
IBOCR= 1 

11  I3CDS= (TCOLT-1) /IDC* 1 
IF(IS-I)  15,15,12 

12  I3CDS=1 

15  DO  31  K-1 , IRCDS 

READ  (5,1,  END=40)  (CARD(I)  ,1=1, IDC) 

IF  (IER)  16,  16,31 

16  L=0 

JS=(K-1)  *IDC*ICOL-ICOLT+1 
JE= JS+IDC-  1 
IF  (IS-  1)  19,19, 17 

17  JE=JS 

19  DO  30  J=JS,JE 

IF (J— ICOL) 20,20,31 

20  CALL  LOC (I BOCR  , J, I J, I ROM, ICOL, IS) 

L-L*  1 

30  A  (I  J)  =CASD  (L) 

31  CONTINUE 
IROCR=IROCH+1 

IF  (IROM-TROCR)  38,35,35 

35  IF  (IS- 1)  37,36,36 

36  ICOLT=ICOLT— 1 

37  GO  TO  11 

38  READ  (5,1,  END=40)  CAED(1) 

IF  (CASD(1) -9.E9)  39,40,39 

39  IBR=2 

40  RETURN 
END 
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SUBROUTINE  SIMQ 
PURPOSE 

03TAIN  SOLUTION  OF  A  SET  OF  SIMULTANEOUS  LINEAR  EQUATIONS, 
AX=B 


USAGE 

CALL  SIMQ(A,B,N,KS) 

DESCRIPTION  OF  PARAMETERS 

A  -  MATRIX  OF  COEFFICIENTS  STORED  COLUMNWISE.  THESE  ARE 
DESTROYED  IN  THE  COMPUTATION.  THE  SIZE  OF  MATRIX  A  IS 
N  BY  N. 

B  -  VECTOR  OF  ORIGINAL  CONSTANTS  (LENGTH  N) .  THESE  ARE 
REPLACED  BY  FINAL  SOLUTION  VALUES,  VECTOR  X. 

N  -  NUMBER  OF  EQUATIONS  AND  VARIABLES.  N  MUST  BE  .GT.  ONE. 
KS  -  OUTPUT  DIGIT 

0  FOR  A  NORMAL  SOLUTION 
1  FOR  A  SINGULAR  SET  OF  EQUATIONS 


REMARKS 

MATRIX  A  MUST  BE  GENERAL. 

IF  MATRIX  IS  SINGULAR  ,  SOLUTION  VALUES  ARE  MEANINGLESS. 

AN  ALTERNATIVE  SOLUTION  MAY  BE  OBTAINED  BY  USING  MATRIX 
INVERSION  (MINV)  AND  MATRIX  PRODUCT  (GMPRD)  . 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

METHOD  OF  SOLUTION  IS  BY  ELIMINATION  USING  LARGEST  PIVOTAL 
DIVISOR.  EACH  STAGS  OF  ELIMINATION  CONSISTS  OF  INTERCHANGING 
ROWS  WHEN  NECESSARY  TO  AVOID  DIVISION  BY  ZERO  OR  SMALL 
ELEMENTS. 

THE  FORWARD  SOLUTION  TO  OBTAIN  VARIABLE  N  IS  DONE  IN 
N  STAGES.  THE  BACK  SOLUTION  FOR  THE  OTHER  VARIABLES  IS 
CALCULATED  BY  SUCCESSIVE  SUBSTITUTIONS.  FINAL  SOLUTION 
VALUES  ARE  DEVELOPED  IN  VECTOR  3,  WITH  VARIABLE  1  IN  B(1), 

VARIABLE  2  IN  B(2), . .  VARIABLE  N  IN  B  (N)  . 

IF  NO  PIVOT  CAN  BE  FOUND  EXCEEDING  A  TOLERANCE  OF  0.0, 

THE  MATRIX  IS  CONSIDERED  SINGULAR  AND  KS  IS  SET  TO  1.  THIS 
TOLERANCE  CAN  BE  MODIFIED  BY  REPLACING  THE  FIRST  STATEMENT. 


SUBROUTINE  SIMQ  ( A, 3, N, KS) 
DIMENSION  A  ( 1)  ,B(1) 


FORWARD  SOLUTION 
TO L=0. 0 


1 


n  n  nnn  non  nnn  non  nnn 
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JJ=-N 

DO  65  J=  1 ,  N 
JY=J* 1 
JJ=JJ*N*1 
BIG A=0 
IT- JJ- J 
DO  30  I=J,N 

SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 
IJ=IT+I 

IF  (ABS  (BIG A)  -ABS  (A  (IJ)  )  )  20,30,30 
20  BIG A=A  (I  J) 

IMAX=I 
30  CONTINUE 


TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

I?  (ABS  (3IGA)  -TOL)  35,35,40 
35  K5= 1 
RETURN 


INTERCHANGE  ROMS  IF  NECESSARY 

40  11  =  J+M*  (J-2) 

IT=IMAX- J 
DO  50  K=J,N 
11=1 1+  N 
12=1 1+TT 
SAVE=A  (II) 

A  ( 1 1)  =  A  (12) 

A  (12)  =SAVE 

DIVIDE  EQUATION  BY  LEADING  COEFFICIENT 

50  A(I1)=A(T1)/BIGA 
SAVE=B  (IMAX) 

B  (IM  AX)  =3  (  J) 

B ( J) =S  A VE/BIG A 

ELIMINATE  NEXT  VARIABLE 

IP(J-N)  55,70,55 

55  IQS=N*(J-1) 

DO  65  IX=JY,N 
IX J=IQS*IX 
IT= J-IX 
DO  60  JX=JY,N 
IX JX=N*  (JX-1)  +  IX 
JJ  Y=IX JX+IT 

60  A(IXJX) =A (IXJX) - (A  (IXJ) *A(JJX) ) 

65  B  (IX)  =B  (IX)  -  (3  (J)  *  A  (IX  J)  ) 

BACK  SOLUTION 
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FItE:  LOC 


FORTRAN  A 


ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


o. 

c 

c 

■ 

c 

c 

SUBROUTINE  LOC 

• 

V* 

c 

PURPOSE 

c 

COMPUTE  A  VECTOR  SUBSCRIPT  FOR  AN  ELEMENT  IN  A  MATRIX  OP 

c 

SPECIFIED  STORAGE  MODE 

m 

L 

c 

USAGE 

V.  , 

c 

r 

CALL 

LOC  (I, J,IR,N,M,MS) 

J  \ 

c 

DESCRIPTION  OF  PARAMETERS 

c 

I 

-  ROM  NUMBER  OF  ELEMENT 

-*■ . 

c 

J 

-  COLUMN  NUMBER  OF  ELEMENT 

c 

IR 

-  RESULTANT  VECTOR  SUBSCRIPT 

%,  ; 

c 

N 

-  NUMBER  OF  SOUS  IN  MATRIX 

y 

c 

a 

-  NUMBER  OF  COLUMNS  IN  MATRIX 

c 

MS 

-  ONE  DIGIT  NUMBER  FOR  STORAGE  MODE  OF  MATRIX 

c 

0  -  GENERAL 

c 

1  -  SYMMETRIC 

c 

2  -  DIAGONAL 

V. 

c 

REMARKS 

c 

p 

NONE 

c 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

1 

c 

p 

NONE 

w 

.  c 

METHOD 

> 

c 

MS=0 

SUBSCRIPT  IS  COMPUTED  FOR  A  MATRIX  WITH  N*M  ELEMENTS 

c 

IN  STORAGE  (GENERAL  MATRIX) 

c 

MS*  1 

SUBSCRIPT  IS  COMPUTED  FOR  A  MATRIX  KITH  N*(N+1)/2  IN 

■ 

c 

STORAGE  (UPPER  TRIANGLE  OF  SYMMETRIC  MATRIX).  IF 

c 

ELEMENT  IS  IN  LONER  TRIANGULAR  PORTION,  SUBSCRIPT  IS 

c 

CORRESPONDING  ELEMENT  IN  UFPER  TRIANGLE. 

c 

MS=2 

SUBSCRIPT  IS  COMPUTED  FOR  A  MATRIX  WITH  N  ELEMENTS 

c 

IN  STORAGE  (DIAGONAL  ELEMENTS  OF  DIAGONAL  MATRIX). 

v' 

c 

IF  ELEMENT  IS  NOT  ON  DIAGONAL  (AND  THEREFORE  NOT  IN 

;  V 

c 

c 

c 

STORAGE) ,  IR  IS  SET  TO  ZERO. 

c 

SUBROUTINE 

LOC (I, J ,IR, N,M, MS) 

c 

it* 

IX=I 

JX=J 

■,*w 

• 

IP  (MS-  1)  10,20,30 

% 

10 

IRX=N*  (JX- 

1)  *IX 

GO  TO  36 

- 

20 

IF  (IX- JX) 

22,24,24 

22  IRX=IX+  (JX*JX-JX)/2 
GO  TO  36 


24  IHX=JX*  (IX*IX-IT) /2 
GO  TO  36 

!* 

•  * 
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30  151=0 

IF(IX-JX)  36,32,  36 
32  I5X=IX 
36  IF=I5X 
RETURN 


B-2.0  FOUR-LAYER  SOLUTION 

The  following  is  an  illustration  of  a  typical  solution  for  a 
typical  four-layer  compressive  layer  calculation.  The  amount 
of  compression  in  the  outer  layer  of  this  example  is  29.0  kpsi 
(0.2  GPa) .  This  value  is  arrived  at  by  adding  together  all 
the  eleventh- solution  values  (i.e.,  0.129981E+05  +  0.606770E+03  + 

r  0. 250532E+04  +  0.128861E+05  * 

0 . 289963E+05) . 

This  solution  was  performed  on  an  IBM  370/138  computer. 


B'  '  i 
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co891011  co891011  datal  b1  cntl  n 2  cntl  m3  cntl  m4  cntl  t 


EXECUTION  BEGINS.  .  . 

ENTE8  A, B,C, D, ALPHA 1 v ALPHA 2, ALPHA3 , ALPHA4 , 
TG1,TG2,TG3,TG4,T0,E,K,MUE,ISIZE 
EXECUTION  BEGINS... 

ENTER  A,  B,C,  D, ALPH A1 , ALPHA2, ALPH A3 , ALPH A4 , 
TG1,TG2,TG3,TG4,T0,E,K,?1UE,ISIZE 
EXECUTION  BEGINS... 

ENTEB  A, 3,C, D, ALPHA  1 , ALPHA2, ALPH A3 , ALPH A4, 
TGl,TG2,TG3,TG4,T0,E,K,HUE,I5IZS 
EXECUTION  BEGINS... 

ENTEB  A, B,C, D, ALPHA  1 , ALPHA2, ALPH A3,  ALPH A 4  , 
TG 1 ,TG2 ,TG3 ,TG4,T0 ,E, X , HUE , I SIZE 
R; 


soln2  b1  cntl  m2  cntl  m3  cntl  m4  cntl  t 


EXECUTION  BEGINS... 


SOLUTION  VALUES 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 


0.0 

0.0 

0.0 

0. 129976 E+05 
0.0 
0.0 

0. 812351E-01 
-0.542325E+04 
-0.542325E+04 
-0. 542325E+04 
0. 129981E+05 


EXE  _  j  HON  TEBHIN  AT  ED 


EXECUTION  BEGINS. • • 

SOLUTION  VALUES 

1  0.0 

2  0.0 

3  -0. 150393E+03 

4  0. 292390E+03 

5  0.0 

6  -0. 125189E-03 

7  0.  182744E-02 

8  -0. 122000E+03 

9  -0. 122000E+03 

10  -0.  278776E+03 

11  0. 606770E+03 
EXECUTION  TERMINATED 


EXECUTION  BEGINS... 

SOLUTION  VALUES 

1  -0. 131200E+04 

2  0.0 

3  -0.416697E+03 
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4 

5 

6 

7 

8 
9 

10 

11 


0.  13544 1 E+04 

0.0 

0. 65446 1 E-0  3 
0. 8465  06  E— 02 
-0.932306E+03 
-0.  13 1 200E+04 
-0.  1036  86  E+04 
0. 250532E+04 


EXECUTION  TEBMI'l ATED 


The  following  files  were  used  in  the  solution  of  the  theoretical 
four-layer  compressive  layer  calculation.  Files  C08 91011 
EXEC  and  SOLN2  EXEC  are  executives  which  call  other  files  and 
control  the  particular  solution. 

Files  CONSTEV8 ,  CONSTEV9,  CONSTElO,  and  CONSTEll  are  FORTRAN 
source  listings  called  by  executive  CO891011. 

Files  SOLN,  MATIN,  MXOUT,  SIMQ,  and  LOC  are  called  by  executive 
SOLN2  and  are  listed  in  the  three-layer  solution  section  of  this 
appendix . 

File  C08 91011  DATA  1  is  the  file  where  the  particular  data  used 
in  this  solution  are  stored. 


FILE:  C3891011  EXEC 
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SCONTROL  OFF  NOMSG 

GLOBAL  TXTLIB  FORTMOD2  30D2EBH  MOD2NEEH  C3SLIB 


FILE  FT05F001 

DISK 

SI 

82 

FILE  FT06F00 1 

TER3 

FILE  FT37F001 
LOAD  CONSTEV8 

DISK 

S3 

S4 

(LR  ECL 

80 

3LOCK 

800 

RECF3 

F3 

START 

FILE  PT05F00  1 

DISK 

S  1 

52 

FILE  FT06F00  1 

TERM 

FILE  FT07F001 
LOAD  C0NSTEV9 

DISK 

S5 

66 

(LR  ECL 

80 

3LOCK 

800 

RECF3 

F3 

START 

FILE  FT05F00 1 

DISK 

SI 

62 

FILE  FT06F001 

TERN 

FILE  FT07F00 1 
LOAD  CONSTE1 0 

DISK 

67 

88 

(LR  ECL 

80 

3  LOCK 

800 

RECFM 

?B 

START 

FILE  FT05F001 

DISK 

SI 

62 

FILE  FT06F001 

TERM 

FILE  FT07F001 

DISK 

69 

610 (LR ECL 

80 

3L0CK 

800 

RECFM 

FB 

LOAD  CONSTE1 1 
START 
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REAL  M,N,M0B,K 
DIMENSION  N  ( 1 1) 

DIMENSION  M (1 1 , 1 1) 

WR ITE  (6,2) 

2  FORMAT  (>  ENTER  A, B,C, D, ALPHA  1 , ALPH A2 , ALPHA3 , ALPHA4 , 
CTG1,TG2,TG3,TG4,T0,E,K,MUE,ISIZE») 

READ  (5,*)  A,B,C,D, ALPHA1, ALPHA2, \LPHA3, ALPHA4 
READ  (5,*)  TG1,TG2,TG3,TG4,T0,E,K,8rJE,ISIZE 
C5  FORMAT  (10E5.5) 

8(1,1)  =3.0 
M  (  1 , 2)  =0.0 
N  ( 1 , 3)  =0. 0 
M  (  1 , 4)  =1.0 
N  ( 1 , 5)  =0.0 
*1(1,6)  =0.0 
3(1,7)=-(1/(D**2)) 

H(1,8)  =0.0 

a  (1,9)  =0.0 
M  (  1 , 10) =0. 0 
8(1,11) =0. 0 
M  (2,  1)  =0.0 
M  (2, 2)  =0.0 
M  (2,3)  =0.0 
M  ( 2 , 4)  =0.0 
M  (2,5)  =0.0 
M  (2, 6)  =0.0 

H  (2,7)  =-  (1/D**2-1/C**2) 

M  (  2 , 8)  =0.0 
M  (2, 9)  =0.0 
M(2,10)  =  1.0 
M (2, 1 1 ) =0. 0 
M  (3, 1)  -  0.0 
M  (3,2)  =0.0 
M  (3,3)  =0.0 
M  (  3 , 4)  =0.0 
M  (3 ,5)  =0.0 
M  (3, 6)  =0.0 
M  ( 3 , 7)  =0.0 
M  ( 3, 8)  =A**2 
N(3,9)  =  (B**2-A**2) 

M(3,10)  =  (C*  *2- 3**2) 

M(3,11)=(D**2-C**2) 

8(4,1)  =0.0 
M  (4 , 2)  =0.0 
M  (4, 3)  =0.0 
H(4,4)  =  ((2*MUE)/E) 

M  (4 , 5)  =0.0 
M  (4 , 6)  =0.0 
8(4,7)  =0.0 
M  (4, 8)  =0.  0 
M(4,9)  =0.0 
M(4,10)=(1/(1.5*K)  ) 


3  (4,  1 1 )  =—  ( 1/E) 
M  (5, 1)  =0.0 
M  (5, 2)  =0.0 


PILE:  C0NST2V3  FORTH  AN  A  ITT  ELECTRO-OPTICAL  PHODOCTS  DIVISION 


a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 


5.3)  =0.0 

5.4)  =0.0 

5.5)  =0.0 

5.6)  =0.0 

5.7)  =0.0 

5.8)  =0.0 

5.9)  =1.0 

5. 10)  =-1.0 

5.11)  =0. 0 
6,1)  =0.0 
6,2)  =0.0 

6.3)  =0.0 

6.4)  =0.0 

6.5)  =0.0 

6.6)  =0.0 

6.7)  =0.0 

6.8)  =-1.0 

6.9)  =1.0 

6. 10)  -0.0 

6.11)  =0.0 

7,  1)  =1 . 0 

7.2)  =0.0 

7.3)  =0.0 

7.4)  =0.0 

7.5)  =0.0 

7.6)  =0.0 

7.7)  =0.0 

7.8)  =0.0 

7.9)  =0.0 
7,  10)  =0.0 

7.11)  =0.  0 
8, 1)  =1.0 
a,2)=i.o 

3.3)  =0.0 

8.4)  =0.0 

8.5)  =0. 0 

8.6)  =0.0 

8,7)  =0.0 

8,3)  =0.0 

3.9)  =0.0 

8. 10)  =0.0 

8.11)  =0. 0 

9. 1)  =1.0 

9.2)  =1.0 

9.3)  =1.0 

9.4)  =0.0 

9.5)  =0.0 

9.6)  =0.0 

9.7)  =0.0 

9.8)  =0.0 

9.9)  =0.0 

9. 10)  =0.0 

9.11)  =0. 0 
10,1)  =0.0 
10,2) =0. 0 
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51(10,3)  =0.0 
M(10,4)  =0.0 
?f  c  1 0 , 5)  =1.0 
51(10, 6)=0.0 
51  ( 10, 7)  =0.0 
M  ( 1  0, 3) =0. 0 
M  (  1 0, 9)  =0.  0 
M  (  10 , 1 0)  =0.0 
N ( 10, 1 1) =0. 0 
51(11,  1)  =0.0 
51(11, 2)  =0.0 
51(11,3)  =0.  0 
51  ( 1 1 ,4)  =0.  0 
M(1 1,5)  =1. 0 
51(11,6)  =1.0 
51(11,7)  =0.0 
51(1 1,3)  =0.  0 
51(11, 9)  =0.0 
3(11,10) =0.0 
51(1 1,11)  =0.0 
1=9 

N ( 1 ) =0.0 
N  (2)  =0. 0 
N (3) =0.0 

N  (4)  =  ( ALPH A4— 3*ALPHA3)  *  (TG4-TG3) 
N  (5)  =0.0 
N  ( 6)  =0.0 
N (7) =0 .0 
■  (8)  =0.0 
N  (  9)  =0 . 0 
N  ( 1 0)  =0.0 
N  (  1 1)  =0.  0 

WRITE  (7,  10)  ISTZE, ISIZE 

10  FOH51  AT  ('  9999 ',214) 

DO  200  .1=1, ISIZE 

WRITE  (7,11)  (3  ( J,  K 1)  ,  K 1=  1 , 7) 

11  FORMAT  (7E10. 4) 

WRITE  (7,12)  (51(J,K1)  ,K1=8,  11) 

12  FORMAT  (4E1 0. 4) 

200  CONTINUE 

WRITE (7,13)  I 

13  FORMAT  (II) 

WaiTE(7,11)  (N  ( K 1 )  ,K1=  1,7) 

WR ITE  (7,12)  (N  (K1)  ,K1=8, 11) 

WRITE (7,14) 

14  FORMAT  (30  ('  •)) 


FILE:  C0MSTEV9  FORTRAN  A 
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HEAL  M,N,3UE,K 
DIMENSION  N  (1 1) 

DIMENSION  3(11*11) 

WRITE  (6*2) 

2  FORMAT  (»  ENTER  A, B,C, D , ALPHA  1 , ALPHA2 , ALPHA3, ALPHA4, 
CTGl,TG2,TG3,TG4,T0,E,X,aUS,ISIZE») 

READ  (5,*)  A,B,C, D, ALPHA1, ALPHA  2, ALPHA  3, ALPHA4 
READ  (5,*)  TGI ,TG2,TG3 ,TG4 ,T0 , E, K, 3UE, ISIZE 
C5  FORMAT  (10E5.5) 

3(1*1)  =0.0 
M  ( 1  *  2)  =0.0 
3(1/3)  =0.  0 
3  (  1 , 4)  =  l.  0 
M  ( 1  *  5)  =0.0 
3(1,6)  =0.0 
3(1/7)=-  (1/(D**2)  ) 

M(1/3)  =0.0 
M  (  1  r  9)  =0.0 
3(1*10) =0. 0 
3(1*11)  =0.0 
N  (2 *  1)  =0.0 
3(2*2)  =0.0 
3(2,3)  =1.0 
M  (2, 4)  =0.0 
3  (2,5)  =0.0 
3(2,6)  =-  (1/(C**2) ) 

3(2,7)  =-(1/D**2-1/C**2) 

M  (2, 3)  =0.0 
M  (2,9)  =0.0 
M  (2,  10)  =0.0 
M (2, 1 1 ) =0. 0 
M(3,1)=0.0 
3(3,2)  =0.0 
3(3,3)  =-1.0 
3(3,4)  =0.0 
3(3,5)  =0.0 
3(3,6)  =  (1/  (C**2)  ) 

3(3,7)  =0.0 
3(3,9)  =0.0 
M(3,9)=1.0 
3(3,10)  =0.0 
3(3, 11) =0.0 
3(4,1)  =0.0 
3(4,2)  =0.0 
3(4,3) =0.0 
3(4,4)  =0.0 
3(4,5)  =0.0 
3(4,6) =0.0 
3(4,7)  =0.0 
3(4,8)  =1.0 
3(4,9)  =-1.0 
3(4,10) =0. 0 
3(4,  11)  =0.0 
3(5,  1)  =0.0 
3(5,2)  =0.0 
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1(5,3)  =-  (2*1UE) 
1(5,4)  =2*MUE 
1(5,5)  =0.0 
a  (5,6)  =0.0 
M (5,7) =0.0 
1  (5,8)  =0.0 
1(5,9)  =0.0 
1(5,10) =1.0 
1(5, 11)  =-1.0 
1(6,1)  =0.0 
1(6,2)  =0.0 
1(6,3)  =0.0 
1(6,4)  =0.0 
1(6,5)  =0.0 
1(6,6)  =0.0 
1  (6,7)  =0.0 
1(6,8)  =  A**  2 
1(6,9)  =  (B**2-A**2) 
1(6, 10)  =  (C**2-B**2) 
1(6,11)  =  (D**2-C**2) 
1(7,  1)  =0.0 


I 


L2 


k 

k  V 

L.  *  ' 


1(7,2)  = 
1(7,3)  = 
1(7,4)  = 
1(7,5)  = 
1(7,6)  = 
1(7,7)  = 
1(7,8)  = 
1(7,9)  = 
1(7,  10) 
1(7,11) 
1(8,1)  = 
1(8,2)  = 
1(3,3)  = 
1(8,4)  = 
1(8,5)  = 
1(8,6) 
1(8,7)  = 
1(8,8)  = 
1(8,9)  = 
1(8, 10) 
1(8,11) 
1(9,1)  = 
1(9,2)  = 
1(9,3) 
1(9,4)  = 
1(9,5)  = 
1(9,6)  = 
1(9,7)  = 
1(9,8)  = 
1(9,9)  = 
1(9,  10) 
1(9,11) 
1(10,1) 
1(10,2) 


0.0 

-  ( ( 1-1UE)  /E) 

0.0 


0.0 

(  (1+1UE)/(E*(B**2))  ) 
0.0 
0.0 

(1/(1. 5*K)) 

=  (1UE/E) 

=0.0 

0.0 

0.0 

( 1-10E) 

- (1-1UE) 

0.0 

(  (1+1UE)/(C**2)  ) 
-((H1UE)/(C**2)) 

0.0 
0.0 
=-1UE 
=  1UE 
1.0 
0.0 
0.0 
0-0 
0.0 
0.0 
0.0 
0.0 
0.0 
=  0.0 
=  0.0 
=  1.0 
=  1.0 


FILE: 


10 

11 

12 

200 

13 

14 


CONSTEV9  FORTRAN  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


M  { 1 0 ,  3 )  =0.0 
M(10,4)  =0.0 
N ( 10, 5) =0. 0 
M ( 10,6) =0. 0 
M ( 10,7)  =0.0 
f!  (10,8)  =0.  0 
.1(10,9)  =0.0 
N(10,  10)  =0.0 
«  (  10,1  1)  =0.0 
11(11,1)  =  1-0 
M ( 1 1 ,2) =1. 0 

n  (  i i , 3) =0- o 
H  ( 1 1 , 4 )  =0.  0 
11(11, 5)  3 1.0 
M  ( 1 1 , 6) =0. 0 
*1(11,7)  =0.0 
fl( 11,8)  =0.0 
M  ( 1 1 , 9) =0. 0 
M  ( 1 1 ,  1 0)  =0. 0 
M( 11, 11) =0.0 
1=9 

N ( 1) =0. 0 
N  (2)  =0.0 
N  (3)  =0.0 
N (4) =0- 0 

N  (5)  =E*  (ALPHA4-ALPHA3)  *(153-191) 
N (6) =0.0 

If  (7)  =  (ALPS  A3 -3*  ALP  HA  2)  *  (TG3-TG1) 
II  (3)  =E*  ( ALPHA4-ALPHA3)  *(TG3-TGl) 
N  ( 9)  =0.0 
N  ( 1 0)  =0.0 
N(1 1)  =0.0 

WRITE  (7,10)  ISIZE, ISIZE 
FORMAT  ( '  9999  *  ,  214) 

DO  200  J=1, ISIZE 

.WRITE  (7,11)  <a(J,K1)  ,K  1=1,7) 

FORMAT  (7E10-4) 

WRITE(7,12)  (H  (J,K1)  ,K  1=8,  11) 
FORMAT  (4E10.4) 

CONTINUE 
WRITE  (7, 13)  I 
FORMAT  (II) 

WRITE  (7,11)  (N  (Kl)  ,K1=  1 ,7) 

WRITE  (7,12)  (N  (Kl)  ,K1=3,  11) 

WRITE  (7,14) 

FORMAT  (80  ('  »)) 

STOP 

END 


B-34 


FILS 

2 

C5 


CONSTEIO  FORTH  AH  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


* 


REAL  M,N,MUE,K 
DIMENSION  N  ( 1 1 ) 

DIMENSION  a  (11,11) 

WRITE  (6,2) 

FORMAT ( ’  ENTER  A, 3 ,C, D , ALP  HA  1 , ALPH A2 , ALPHA3 , ALPHA4 , 
CTGl,TC.2,TG3,TG4,T0,Z,K,MtI2,ISIZE‘) 

READ  (5,*)  A,B,C,D,ALPHA1, ALPHA2, ALPHA3, ALPHA4 
READ  (5,*)  TGl,TG2,TG3,TG4,T0,E,X,i1UE,I3IZE 
FORMAT  (10E5.5) 

M(1,1)=0.0 
M  ( 1 ,2)  =0.  9 
M  ( 1 , 3)  =0.0 
M  (  t  ,  4)  =1.0 
M  ( 1, 5)  =0.0 
M(  1,6)  =0.0 
M(1,7)=-(1/(D**2)) 
a  (i#8)  =0.0 

M  (  1 , 9)  =0.0 

;i(  i ,  10)  =0.0 
8(1,11) *0. 0 
M  (2,  1)  =0.0 
M  (2, 2)  =0.0 
M  (2 , 3)  =1.0 
M  (2, 4)  =0.0 
M  (2,5)  =0.0 
M  ( 2 , 6)  =-  (1/  (C**2)  ) 

M(2,7)=-  ( 1/D** 2- 1/C**2) 

M  (2, 8)  =0.0 
M  (2, 9)  =0.0 
M(2,10)=0.0 

a ( 2, i i) =o. o 

M (3, 1) =0.0 
N  ( 3 , 2)  =0.0 
M  (3, 3)  =0.0 
M  ( 3 , 4)  =0.0 
M  (  3 , 5)  =0.0 
M  (3,6)  =0.0 
M  ( 3 , 7)  =0.0 
M  (3 , 8)  =  A**2 
M  (3,9)  =  (E**2-A**2) 

M  (3,  10)  =  (C**2-B**2) 

M(3,11)= (D**2-C**2) 

M  (4 , 1)  =0.0 
M  (4 , 2)  =0.0 
M  (4, 3)  =-1.0 
8(4,4)  =0.  0 
M  (4, 5)  =0.0 
8(4,6)  =  (1/(B**2)  ) 

M (4, 7) =0.0 
M  (4 , 8)  =0.0 
M  (4, 9)  =1.0 
M  (4 , 1 0)  =0.0 
8(4,11)  =0. 0 
M  (5,  1)  =-  1.  0 
M  (5, 2)  =0.0 


•i 


.  ,■> 

* 

* 


FILE:  CONSTE10  FOETH  All  A 


ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


1(5,3)  =0.0 
1(5,4)  =0.0 
«  { 5,  5)  =0.0 
1(5,6)  =0.0 
1(5,7)  =0.0 
1(5,3)  =0.0 
1(5,9)  =1.0 
1(5, 10) =0. 0 
1(5,11) =0. 0 
1(6,1)  =0.0 
1(6,2)  =0.0 
1(6,3)  =-2*10E 
M (6,4) =2*10S 
1(6,5)  =0.0 
11(6,6)  =0.0 
1(6,7)  =0.0 
1(6,8)  =0.0 
N (6,9) =0.0 
1(6,  10)  =1.0 
1(6, 11)  =-1.0 
1(7,  1)  =0.0 
1(7,2)  =0.0 
1(7,3)  =  ((2*1UE)/E) 

1(7,4)  =0.0 
1(7,5)  =0.  0 
1(7,6)  =0,0 
1(7,7)  =0.0 
1(7,8)  =0.0 
1(7,9)  =  (1/(1. 5*K)) 

1(7,  10)=- (1/E) 

1(7,  11)  =0.  0 
1(8,1)=-((2*10E)/(E)) 
1(8,2)  =0.0 
1(8,3)  =0.0 
M  (8, 4)  =0.0 
1(8,5)  =0.0 
1(8,6) =0.0 
1(8,7)  =0.0 
1(8,8)  =1/E 
1  (8,9)  =-  (1/(1.  5*K)  ) 
1(8,10) =0. 0 
1(8,11) =0. 0 
1(9,1)  =0.0 
1(9,2)  =0.0 
1  (9 , 3)  =  ( 1-1UE) 

1(9,4) =-  (1-1WE) 

1(9,5)  =0.0 

1(9,6)  =  (  ( 1  ♦1UE)  /  (C**2)  ) 
1(9,7)  =-  (  ( 1  +  1UE)  /  (C**2)  ) 


S  L 


1(9,8)  =0.0 
1(9,9)  =0.0 
1(9,10)  =-1CJE 
1(9,11)  =10E 
1(10,1)  =0.0 
1(10, 2)  =  1.0 


A 


i 


PILE:  CONSTE 10  FORTRAN  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


.1(10,3)  =0.  0 
H(10,4) =0. 0 
1(10,5) =0. 0 
1(10,6)  =0.  0 
1(10,7) =0.0 
1(10,9) =0. 0 
1(10,9)  =0.0 
1(10,10) =0.0 
1(10,11)  =0.0 
1(1 1,1) =0. 0 
1(11,2)  =1.0 
1(11,3) =0. 0 
1(11,4)=0.0 
1(11, 5)  =  1.0 
1(11,6) =0.0 
1(11,7) =0.0 
1 ( 1 1 »  3) =0. 0 
*1(11,0)  =0.  0 
1(11,10) =0.0 
1(11,11) =0.0 
1=9 

11(1)  =0.0 
N (2) =0. 0 
N  ( 3)  =0.0 
N (4) =0.0 
N  (5)  =0.0 

N  (6)  =E*  (ALPHA4-ALPHA3)  *  (TG1-TG2) 
N  (7)  =  ( ALPH  A3-3*ALPH A2)  *  (TG1-TG2) 
N(9)  =  (3*ALPHA2-ALPHA1)  *  (TG1-TG2) 
N(9)=E*  (ALPHA4-ALPHA3)  *  (TG1-TG2) 
N  ( 10)  =0.  0 
N  ( 1 1)  =  0.  0 

WRITE (7, 10)  ISIZE, ISIZE 

10  FOH1AT  ('  9999',  214) 

DO  200  J=1, ISIZE 

WRITE  (7,11)  (1  (J,K1)  ,K1=1,7) 

11  FOR1AT  (7E10.4) 

WRITE  (7,12)  (1  (J,K1)  ,K1=8,  11) 

12  FOR1AT  (4E10.4) 

200  CONTINUE 

WRITE (7,13)  I 

13  POR1AT  (II) 

WRITE(7,11)  (N(K1)  ,K1=1,7) 

WRITE(7, 12)  (N  (Kl)  ,K  1=8, 1 1 ) 

WRITE (7, 14) 

14  FORI  AT  (80  ( '  ')  ) 


PILE:  C0NSTE1 1  FORTRAN  A 


ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


REAL  M,N,MUE,K 
DIMENSION  N  ( 1 1) 

DIMENSION  S  (11* 1 1) 

HRITE(6,2) 

FORMAT  (•  ENTER  A, B,C, D, ALP  HA  1 , ALPHA2, ALPUA3 , ALPHA4 , 
CTGl,TG2,TG3,TG4,T0,E,K,aU2,ISIZE») 

READ  (5,*)  A,B,C, D, ALPHA  1,  ALPBA2, ALPHA3, ALPHA4 
READ  (5,*)  TGl,TG2,TG3,TG4,T0,E,K,MUE,ISIZE 
FORMAT  (10E5.5) 

M(1,1)=0.0 
M  ( 1 ,-2)  =0.0 
M  (1,3)  =0.0 
M(  1,4)  =1.0 
M(  1,5)  =0.0 
M  ( 1 , 6)  =0.0 
M(1,7)  =- (1/ (D**2) ) 

M  ( 1 , 8)  =0.0 
M  ( 1 , 9)  =0.0 
M ( 1 , 10) =0. 0 
M ( 1 , 1 1) =0. 0 
M  (2,  1)  =0.0 
M(2,2)  =0.0 
il  (2,3)  =- 1-  0 
M  (2, 4)  =0.0 
M  (2, 5)  =0.0 
U  (2, 6)  =  ( 1/  (C**2)  ) 

M  (2,7)  =  ( 1/D**2-  1/C**2) 

M{2,8)  =0.0 
M  (2, 9)  =0.0 
M (2, 10) =0. 0 
M(2,  11)  =0.  0 
M  (3, 1)  =0.0 
M  (3, 2)  =-1.0 
H  (3, 3)  =1.0 
M  (3,4)  =0.0 
H(3,5)  =  (1/(B**2)) 

M  (3, 6)  =-  (1/(B**2)  ) 

M  ( 3 , 7)  =0.0 
M  (3 , 8)  =0.0 
M(3,9)=0.0 
M (3, 10) =0. 0 
M (3, 11) =0. 0 
H  (4,  1)  =- 1.  0 
M  (4, 2)  =1  •  0 
M  (4,  3)  =0.0 
M  (4, 4)  =0.0 
«(4,5)=-(1/(A**2)) 

3(4,6) =0.0 
•1(4,7)  =0.0 
M  (4 , 8)  =0.0 
M  (4, 9)  =0.0 
3(4,10)  =0.0 
1(4, 11) =0. 0 
M(5,1)=0.0 
3(5,2)  =0.0 


B-38 


PILE:  CONSTE 1 1  FORTRAN  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


u  % 


*  u 


It  u 

'j 


3(5,3)  =0 
3(5,4)  =0 
3(5,5)  =0 
3(5,6)  =0 
3(5,7)  =0 
3  (5, 3)  =S 
3  (5 , 9)  =  ( 
3(5,10)= 
3(5,11)= 
M(6,1)  =  ( 
3(6,2)  =- 
3(6,3)  =C 
3(6,4) =( 
3(6,5)  =- 
3(6,6)  =C 
3(6,7)  =C 
3(6,9)=- 
3(6,9)  =J 
3(6,10)= 
3(6,11)= 
3  (7,  1)  =( 
3(7,2)  = 
3(7,3)=- 
3(7,4)  =( 
3(7,5)=f 
3(7,6)  = 
3(7,7)=! 
3(7,8)=! 
3(7,9)  = 
3 (7,  1 0) : 
3(7,11) 
3(8,1)  = 
3(8,2)=! 
3(8,3)  = 
3(8,4)  = 
3(8,5)  = 
3(8,6)  = 
3(8,7)  = 
3(8,3)  = 
3(8,9)  = 
3(8,  10) 
3(8,11) 
3(9,1)  = 
3(9,2)  = 
3(9,3)  = 
3(9,4) = 
3(9,5)  = 
3(9,6)  = 
3(9,7)  = 
3(9,8)  = 
3(9,9)  = 
3(9,10) 
3(9,11) 
3(10,1) 
3(10,2) 


0.0 

0.0 

0.0 

0.0 

0.0 

A**2 

(B**2-A**2) 

=  (C**2-B**2) 

=  (D**2-C**2) 
(1-3UE) 

-  (1-3UE) 

0.0 

0.0 

-  (  (1*MOE)/(A**2)  ) 

0.0 

0.0 
-HUE 
=MUE 
I  =0.C 
1=0.0 
=  0.0 

=  (1-3UE) 

=-  (1-3UE) 


0.0 

(  (U3UE)/(B**2)  ) 
0.0 
0.0 

-  (30E) 

=  3US 

=  0.0 
0.0 
0.0 

(1-30E) 

-  (1-3UE) 

0.0 

(  (1  +3UE)/  (C**2)  ) 

-  ( ( 1+3UE)  /  (C**2)  ) 

0.0 

0.0 
=-30E 
=  3CJE 
0.0 
0.0 
0.0 
0.0 

(  (2*30E)/(A**2)  ) 
0.0 
0.0 
1.0 
-1.0 
I  =0.0 
I  =0.0 
1=0.0 
|=-(2*3UE) 
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FILE:  CONSTE 1 1  FORTRAN  A 


ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


M  (10,  3)  =2*MUE 
1(10,4) =0. 0 
.1(10,5)  =0.  0 
M(10,6)=0.0 
1(10,7)  =0-0 
1(10,3) =0. 0 
1(10,9)  =1.0 
1(10,10)=-1.0 
1(10,11)=0.0 
1(11,  1)=0.  0 
1(11,2)  =0.0 
1(11,3)=-  (2*MUE) 

M  ( 1 1 , 4) =2*MUE 
1(11,5) =0. 0 
1(11,6) =0. 0 
1(11,7)  =0.0 
1(11,8)  =0.  0 
1(11,9) =0.0 
1(11, 10)  =1-0 
1(11,11)  =-1.0 
1=9 

N ( 1) =0. 0 
N  ( 2)  =0.0 
N  (3)  =0.0 
N  (4)  =0.0 
N (5) =0. 0 

N  (6)  =  (E* ( ALPH A2- ALPHA 1 )  *  (TG2-T0)  ) 

N  (7)  =  (E*  (AL PH A3— ALPHA 2) * (TG2-T0) ) 

N  (8)  =  (E*  ( ALPH A 4- ALPH A3) * (TG2-T0) ) 

N  (9)  =  ( E*  (ALPHA2-ALPHA1)  *  (TG2-T0)  ) 

N  ( 1 0)  =  (E*  (ALPH A3- ALPHA 2)  *  (TG2-T0) ) 
N  ( 1 1)  =  (E*  (ALPHA4-ALPHA3)  *  (TG2-T0)  ) 
WRITE  (7,10)  ISIZE,ISIZE 

10  FORMAT  ('  9999', 214) 

DO  200  J= 1 ,ISIZE 

WRITE  (7,11)  (1  (J,K1)  ,K1=1,7) 

11  FORMAT  (7E10. 4) 

WRITE  (7,12)  (1  (J,K1)  ,X1=3,  11) 

12  FORMAT  (4E10. 4) 

200  CONTINUE 

WRITE  (7, 13)  I 

13  FORMAT  (II) 

WRITE  (7,11)  (N  ( K 1 )  ,K1=  1 ,7) 

WRITE  (7,  12)  (N(K1)  ,X1=3,11) 

WRITE  (7,  14) 

14  FORMAT  (30  (»  •)  ) 

STOP 


FILE:  S0LN2  EXEC  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


SCONTROL  OFF  NOMSG 

GLOBAL  TXTLTB  FORT MOD 2  MOD2EEH  MOD2NFEH  CMSLIB 
6IF  63  =  P  6  GO  TO  -A 
SIF  83  =  T  5G0T0  -B 
-C 

FILEDEF  FT05P001  DISK  81  62 

LOAD  SOLN 

START 

FILEDEF  FT05F001  DISK  63  64 

LOAD  SOLN 

START 

FILEDEF  FT05F001  DISK  65  66 

LOAD  SOLN 

START 

FILEDEF  FTO5FO0 1  DISK  67  58 

LOAD  SOLN 

START 

5EXIT 

-A 

FILEDEF  FT06F001  PRINTER 
5GOTO  -C 
-B 

FILEDEF  FT06F001  TERM 
6GOTO  -C 


PILE:  CO891011  DA7A1  A  rmT  77. ZC^TO-OPTICAL  PRODOCTS  DIVISION 

6.57E-4,8.55E-4,  2.  102-3,  2.  59  E- 3,  2.  J1Z-6,  1 . 04E-6, 2. 5E-6 , 0. 5E-6 
865,725,900,  1120,25,  1.0577,5*6,.  17,  1  1 
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Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 


Table  C— 1 .  Preforra/Fiber  Fabrication  Data  (continued) 


TJ 

•D 

41 

m 

TJ 

4) 

U 

Vi 

C  CQ 

Ij 

0) 

4) 

•w  0) 

4) 

4J 

4J 

<0  4)  H 

4J 

4J 

4J 

4J  O' ^5 

*J 

AS 

(0 

c  n  .o 

IQ 

J3 

.C 

0  <0  3 

jC 

CO 

(0 

UH£ 

CO 

<u  c 
0  3  0 
•H  O’ 

4J  — i  4J 

fl  u  ® 
■U  IQ  4> 

co  »w  4J 


4J  4) 

i 

a 

m 

»  c 

4)  0 

a 

i 

4> 

X 

1 

1 

V 

X 

-C 
C  4J 
J  O'  — 

<0  c  s 

vj  4>  — 

a  j 


o 

1  CN 


VO  CN 

I  I 

r~  in 

o  CN 

vo  I  p~  | 
o  o 

O'  O' 

I"  I" 


4-1 

c 

, — . 

4) 

•H 

E 

CO 

4J 

RS 

4) 

Q 

Vi 

O 

Pi 

+ 

1  1 

U 

u 

O' 

c 

•*4 

•H 

•M 

H 

(0 

•ri 

£ 

Q 

4) 

•*4 

CO 

o 

E 

+ 

W 

IQ 

<H  1 

a 

U 

Cu 

o 

• 

Vi 

• 

Vi 

CN 

• 

Vi 

• 

Vi 

e 

Vi 

• 

Vi 

e 

Vi 

e 

Vi 

e 

u 

« 

u 

0 

o 

£ 

0 

o 

:> 

• 

£ 

0 

o 

£ 

0 

o 

£ 

0 

o 

> 

0 

o 

£ 

0 

o 

£ 

0 

o 

£ 

o 

o 

:> 

0 

o 

£ 

■<r  oo 
cn  cn 


as 

o 

r> 

VO 

00 

O' 

w 

co 

in 

in 

00 

O' 

O' 

O' 

O 

o 

CM 

CN 

CN 

CN 

CN 

CN 

CN 

co 

ro 

O 

o 

o 

o 

o 

O 

© 

o 

© 

N1 

n* 

n' 

to 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

cn 

Ed 

w 

u 

u 

u 

w 

w 

Ed 

Ed 

n  it  n 

CC  *H  -w 
CJ  CO  E-* 


C-4 


Chemical  Etching  GR  =  Grinding 

silica  FP  =  Fire  Polishing 

Ti-doped  Silica  OD  =  Outside  Deposition 


Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Table  C-l.  Preforra/Piber  Fabrication  Data  (continued) 
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Chemical  Etching  GR  =  Grinding 

Silica  FP  =  Fire  Polishing 

Ti-doped  Silica  OD  =  Outside  Deposition 


Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Chemical  Etching  GR  »  Grinding 

Silica  FP  =  Fire  Polishing 

Ti-doped  Silica  OD  ■  Outside  Deposition 


Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 


Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Chemical  Etching  GR  =  Grinding 

Silica  FP  =  Fire  Polishing 

Ti-doped  Silica  OD  =  Outside  Deposition 


Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Table  C-l.  Preform/Fiber  Fabrication  Data  (continued) 
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Chemical  Etching  GR  =  Grinding 

Silica  FP  =  Fire  Polishing 

Ti-doped  Silica  OD  =  Outside  Deposition  *Under  Chem  Etching 
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Table  C-2.  Graded  Index  Preform/Fiber  Fabrication  Data  (continued) 
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Legends:  0D  *  Outside  Deposition  Si  =  Silica  Layer 

GR  =  Center  Grinding  Ti  =  Ti-Doped  Silica  Layer 

FP  =*  Fire  Polishing  SF  =  Static  Fatigue 

CE  *  Chemical  Etching 


Table  03.  Silica  Rod  Preform/Fiber  Fabrication  Data 


Legends:  OD  «  Outside  Deposition  CE  *  Chemical  Etching 

GR  =  Center  Grinding  Si  *  Silica  Layer 

FP  *  Fire  Polishing  Ti  *  Ti-Doped  Silica  Layer 


D-1.0  DATA  PROCESSING  COMPUTER  OPERATION 

The  following  is  an  illustration  of  a  typical  data  processing 
computer  operation.  This  particular  example  displays  the  process 
of  raw  static  fatigue  data. 


PILE:  TEEM  1  FILE  A  ITT  ELECTRO-OPTIC Al  PRODUCTS  DIVISION 


data 

-If  at  any  tine  daring  the  operation  of  this  program  you  want  to 
leave  simply  type  "OFF"  and  the  hit  the  "ENTER"  key. 

If  at  any  time  while  in  a  section  you  want  to  restart  that 
section  simply  type  "R"  or  "RESTART"  and  hit  the  "ENTER"  key. 

-Is  the  data  in  a  file? 

(Answere  Fes  or  No) 


-What  is  the  name  and  filetype  of  your  data  file? 
static  data 

-Lets  take  one  last  look  at  the  data. 

PLEASE  MAKE  SORE  THE  FOLLOWING  CCNTITIGNS  ARE  MET: 

1)  No  negative  or  zero  data  entered 

2)  If  you  wish  to  include  zero  values  of  enter  a  small  value 
tut  not  zero  in  its  place. 

3)  Enter  no  more  than  100  points. 

4)  There  is  one  and  only  one  comma  between  data  values. 

-To  get  out  of  the  EDGAR  mode  (and  have  this  program  continue  running) 

1)  Hit  ENTER 

2)  Type  the  word  "FILE" 

3)  Hit  ENTER 


HIT  ENTER  WHEN  READY 
(enter) 

ES— 040489E  * 

800123-1 

STF1 

STANDARD  STATIC  FATIQOE  TEST  (EXAMPLE) 


445.0 

4.90 

20.8 

72.0 

385.0 

5.00 

20.6 

72.0 

410.0 

4.90 

20.4 

72.0 

450-0 

4.90 

20.8 

72.0 

390.0 

5.00 

20.8 

72.0 

6000.0 

5.  10 

21.2 

83.0 

107520.0 

5.00 

20.0 

83.0 

600.0 

5.90 

20.6 

83.0 

40440.0 

4.90 

19.8 

83.0 

6120.0 

5.30 

20.4 

83.0 

593460-0 

5.  10 

20.8 

95.0 

402600.0 

4.70 

20.2 

95.0 

239640.0 

5.30 

20.8 

95.0 

163200.0 

5.50 

21.4 

95.0 

file 


PREFORM  NO.  ES-040489B 


FIBER  NO.  800123-1 


FILE:  TE8M1  FILE  A  ITT  ELECTEQ-GPTICAL  PRODUCTS  DIVISION 


STATIC  TEST  DATA 

STANDARD  STATIC  FATIQUE  TEST  (EXAMPLE) 


FD 

OD 

MD 

STBESS 

TIME  TO  FAILURE 

(MIL) 

(NIL) 

(MIL) 

(G  PASCAL) 

(SEC) 

4.90 

20.8 

72.0 

3.8017 

4  45. 

5.00 

20.6 

72.0 

3. 8877 

385. 

4.90 

20.4 

72.0 

3.8182 

410. 

4.90 

20.8 

72.0 

3.8017 

450. 

5.00 

20.8 

72.0 

3.  8  7  S3 

390. 

5.  10 

21.2 

83.  0 

3.5240 

6000. 

5.00 

20.0 

83.0 

3.4951 

107520.. 

5.90 

20.6 

83.0 

4. 1004 

600. 

4.90 

19-8 

83.0 

3.4319 

40440. 

5..  30 

20.4 

33.0 

3.6905 

6120. 

5.  10 

20.8 

95.0 

3.  1710 

593460. 

4.70 

20.2 

95.0 

2.9375 

402600. 

5..  30 

20.3 

95.0 

3.  2953 

239640. 

5.,  50 

21.4 

95.0 

3.4021 

163200. 

you  want 

furter  data 

reduction? 

PREFORM  NO.  ES-040489B  FIBEB  NO.  800123-1 

STATIC  TEST  DATA 

STANDABD  STATIC  FATIQUE  TEST  (EXAMPLE) 


COEFFICIENT 

OF 

DETERMINATION 


LINEAR 

T  = 

. 1778E+07 

♦ 

-. 4645E+06 

S 

0.664 

EXPONENTIAL 

T  » 

.  1457E  +  18 

EXP 

(  -8.491 

S) 

0.848 

LOGARITHMIC 

T  * 

. 2216E+07 

♦  -. 

1652E+07  LN 

s 

0.685 

LOG-LOG 

LOG  T  = 

-29.58 

S  + 

20.  29 

0.839 

USING  LOG-LCG  EQUATION 

STRESS  AT  TIME  TO  FAILURE  OF  10,000  SEC.  =  3.5555  GIGA  PASCAL 
STRESS  AT  TIME  TO  FAILURE  OF  1,000  SEC.  =  3. 8434  GIGA  PASCAL 


N  VALUE  =  29.5776 


PILE:  TEBMl  FILE  A  ITT  ELECTRO-OPTICAL  PBODUCTS  DIVISION 


-Do  you  want  that  sent  to  the  printer? 

y 


-Please  enter  name  of  the  person  to  yet  or  receive  all 
the  printed  material 


final 


-Do  you  want  a  plot  of  that  data? 

y 


-Do  you  have  more  data  to  enter? 
n 


0001  FILE  CHANGED 

-There  are  1  sets  of  multiple  plots  awaiting  output 


H; 

■sg  op  please  attach  plotter  tape.... thank  you. ..pete 
B; 

TAPE  181  ATTACHED 

08:23:49  8SG  FBOH  OPEB  :  TAPE  IS  READY 

move 

B; 


tape  run 
B; 


asg  op  thank  you 


The  following  are  the  computer  source  listings  used  in  the  previous 
example.  The  listings  include: 

DATA  EXEC  -  An  executive  that  interacts  between  the 
laboratory  technician  and  the  various 
load  modules  and  data  files 

DATA  FORTRAN  -  A  FORTRAN  source  listing  of  the  program 
used  in  processing  the  data 

STATIC  DATA  -  An  example  of  a  typical  static  fatique 
raw  date  file 


D-6 


PILE:  DATA  EXEC  A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


SCONTRCL  OFF  NOMSG 

* 

* 

*  GENERAL  PURPOSE  EXECUTIVE  FOR  CALCULATION,  PRINTING  AND  PLOTTING 

*  CALCULATED  RESULTS- ... P. H.  PBIDEAUX  APRIL  22,  1980 

*  UPDATED  MAY  22,  1960 

* 

* 

GLOBAL  TXTLIB  FORTMOD2  MOD2EEH  MQ02SEEH  CMSLIB 

*  EXEC  FCBTH  DATA 

ERASE  DATA  PLOT 

6P  =  0 

SO  =  0 

-DATA1 

FILE  FT01F001  DISK  CCPLOT  FILE  (LRECL  80  BLOCK  800  PERM  RSCFM  FB 

FILE  FT02P00 1  DISK  PLOT  DATA  (LRECL  80  BLOCK  800  PERM  BECFM  FB 

FILE  FT03F001  DISK  DEBUG  TEMP  (LRECL  132  BLOCK  792  PERM  BECFM  F3 

FILE  FT07F001  DISK  DATA1  DATA  (LRECL  132  BLOCK  792  PERM  RECFM  FB 

FILE  FT08F00 1  DISK  DAT A2  DATA  (LRECL  132  BLOCK  792  PERM  BECFM  FB 

SV  -  1 
CLEAR 
SBEGTYPE 

-If  at  any  time  daring  the  operation  of  this  program  you  want  to 
leave  simply  type  "CPF"  and  the  hit  the  "ENTER"  key. 

If  at  any  time  while  in  a  section  you  want  to  restart  that 
section  simply  type  "R"  or  "RESTART"  and  hit  the  "ENTER"  key- 

SEND 

SBEGTYPE 

-Is  the  data  in  a  file? 

SEND 

-DATA4 

SBEGTYPE 

(Answers  Yes  or  No) 

SEND 

SREAD  VARS  61 
SR  a  SI 

6IF  .SR  a  .  SGOTO  -DATA4 
SIF  SR  =  RESTART  SGOTO  -DATA  1 
SIF  SR  a  R  SGOTO  -DAT A 1 
SIF  SR  =  NO  SGOTO  -ERASE 
SIF  SB  =  N  SGOTO  -ERASE 
SIF  SR  =  YES  SGOTO  -FILE 
SIF  SR  =  Y  SGOTO  -FILE 
SIF  SR  =  CFF  SGOTO  -EXIT 
SGOTO  -DATA4 
-ERASE 
CLEAR 

ERASE  DATA  DATA 
SBEGTYPE 

After  you  get  done  with  this  page  you  will  he  put  into  the  edit  mode 
of  a  new  file.  First  enter  the  Preform  ID  Number  on  the  first  line 
then  the  Fiber  ID  number  on  the  second  line.  The  type  of  data  in  the 


FILE:  DATA 


EXEC 


A 


ITT  ELECT EQ-QPTI CAL  PRODUCTS  DIVISION 


file  * DYN 1 1  for  standard  dynamic  data,  ' BND 1  *  for  bending  dynaaic 
data  1 STF 1 '  for  standard  static  fatigue  data  in  the  next  line.  In  the 
forth  line  type  in  remarks.  Follow  this  with  the  data. 


SEND 

SBEGTYPE 

ES— 040489B 
800123-1 

STF1 

STANDARD  STATIC 

FATIQUE 

TEST 

(preform  number. .. first  12  characters) 
(fiber  number.. .. first  12  characters) 
(data  type  STF  1 ,  DYN  1 ,  or  BND  1) 
(remarks. -.first  40  chracters) 

4  45 

4.90 

20.8 

72.0 

385 

5.00 

20.6 

72.0 

410 

4.90 

20.4 

72.0 

450 

4.90 

20.8 

72.0 

(data. ..as  described  below) 

390 

5.00 

20.8 

72-0 

-for  static  fatigue  results  enter  in  this  order 

time  to  failure...  fiber  diameter. ..  coating  diaaeter.. 
and  mandrel  diameter,  (all  units  to  be  in  ails) 

-for  dynamic  strength  results  enter 

total  load  on  fiber...  and  fiber  diameter. 

-for  bending  test  results  enter 

micrometer  reading...  coating  diameter...  and  fiber  diameter 


SEND 
—  ERASE  1 
S  A  =  DATA 
SB  -  DATA 
SGOTO  -EDIT  1 
-FILE 
SBEGTYPE 

-Hhat  is  the  name  and  filetype  of  your  data  file? 
SEND 

SBEAD  VA3S  Si  S2 
SA  -  SI 
SB  =  62 

SIP  -SA  *  .  SGOTO  -FILE 
SIP  SA  =  CFF  SGOTO  -EXIT 
SIP  SA  =  B  SGOTO  -CATA1 
SIP  SA  -=  BESTABT  SGOTO  -DATA  1 
SIP  .SB  =  .  SGOTO  -FILE 
-EDIT2 
SV  =  2 
SBEGTYPE 

-Lets  take  one  last  look  at  the  data. 

SEND 
—EDIT  1 

PILE  FT04F00 1  DISK  SA  SB 
S3EGTYPE 

PLEASE  HAKE  SURE  THE  FOLLOWING  CCNTITIONS  ABE  NET: 


PILE:  DATA  EXEC  A  ITT  ELECTRO— OPTICAL  PRODUCTS  DIVISION 


1)  No  negative  or  zero  data  entered 

2)  If  you  wish  to  include  zero  values  of  enter  a  small  value 
but  not  zero  in  its  place. 

3)  Enter  no  more  than  100  points. 

4)  There  is  one  and  only  one  comma  between  data  values. 

-To  get  out  of  the  EDGAR  mode  (and  have  this  program  continue  running) 

1)  Hit  ENTER 

2)  Type  the  word  "PILE" 

3)  Hit  ENTER 

SEND 
—LOOP  1 
6BEGTYPE 

HIT  ENTER  WHEN  HEADY 
SEND 

SREAD  VARS  51 
SR  =  SI 

SIP  .SR  *  .  SGOTO  -EDIT3 
SIP  SR  =  OFF  SGOTO  -EXIT 
SIP  SR  -  B  SGOTO  -DATA  1 
SIP  SR  =  RESTART  SGOTO  —DATA  1 
SGOTO  -LOOP  1 
-EDIT3 

EDGAR  SA  SE  A 

SIP  SV  =  1  SGOTO  -EDIT2 

CLEAR 

PILE  FT04P001  DISK  SA  SB  A  (LRECL  80  BLOCK  800  PERM  RECFM  FB 

*  LOAD  DATA 

*  GENMOD  DATA 

*  START 
DATA 

TYPE  DAT A1  DATA 

-LOOP2 

SBEGTYPE 

-Do  you  want  furter  data  reduction? 

SEND 

BREAD  VARS  SI 
SR  -  SI 

SIP  -SR  =  .  SGOTO  — LOOP2 
SIF  SR  *  NC  SGOTO  -LOOP3 
SIP  SB  *  N  SGOTO  -LOCP3 
SIP  SR  =  YES  SGOTO  -CUBVE 
SIP  SR  *  Y  SGOTO  -CURVE 
SIP  SR  =  OFF  SGOTO  -EXIT 
SIP  SB  *  B  SGOTO  -BATA  1 
SIP  SB  =  BESTART  SGOTO  -DATA  1 
SGOTO  -LCOP2 
-CURVE 
CLEAR 


PILE:  DATA 


EXEC 


A  ITT  ELECTRO-OPTICAL  PRODUCTS  DIVISION 


TYPE  DATA2  DATA 

-LOOP3 

SBEGIYPE 

-Do  you  want  that  sent  to  the  printer? 

SEND 

SBEAD  VARS  £1 

CLEAR 

SB  =  SI 

SIF  .SB  =  .  SGOTC  -L00P3 
SIF  SB  =  OFF  SGOTO  -EXIT 
SIF  SB  =  Y  SGOTC  -PBINT1 
SIF  SB  =  YES  SGOTO  -PRINT1 
SIF  SR  =  NC  SGOTO  -LOOP4 
SIF  SB  =  N  SGOTO  -LOOP4 

SIF  SB  =  B  SGOTO  -DATA  1 

SIF  SB  =  EESTABT  SGOTO  -DATA  1 
SGOTC  -LOCP3 
-PRINT1 

SIF  SC  NE  0  SGOTC  -PBINT3 
SO  =  1 

CP  SET  aSG  OFF 

-PBINT2 

SBEGTYPE 

-Please  enter  name  of  the  persoa  to  get  or  receive  all 
the  printed  material 


SEND 

SBEAD  VAHS  SI 
6D  =  SI 

SIF  -SD  =  .  SGOTC  -PBINT2 
SIF  SD  =  OFF  SGOTO  -EXIT 
.SIF  SD  »  R  SGOTO  -DATA  1 
SIF  SD  =  BESTABT  SGOTO  -DAT A 1 
-PBINT3 

CP  SPOOL  FBINTEB  COHT  HOLD 

PRINT  DATA  1  DATA 

CP  SPOOL  PRINTER  NOCONT  BOLD 

PBINT  DATA2  DATA 

EBASE  DATA  1  DATA 

ERASE  DATA2  DATA 

-LOOP4 

SBEGTYPE 

-Do  you  want  a  plot  of  that  data? 
SEND 

SBEAD  VABS  SI 
SB  ■  SI 

SIF  .SB  =  .  SGOTO  -LCOP4 
SIF  SB  =  OFF  SGOTO  -EXIT 
SIF  SB  *  Y  SGOTO  -PLOT 
6 IF  SB  *  YES  SGOTO  -PLOT 
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FILE:  DATA 

EXEC 

A 

ITT  ELECT BO— OPTICAL  PRODUCTS  DIVISION 

SIF  SB  =  NC 

SGOTC 

-LOOPS 

SIF  SB  *  N 

SGOTO 

-LOOP5 

SIF  SR  =  B 

SGOTO 

-DATA  1 

SIF  SB  =  EESTAST  SGOTO  -DATA  1 

SGOTC  -LOOP4 

-PLOT 

SP  =  6P  ♦  1 

FILE  FT05F001  DISK  JUNK5  TEMP  (LBECL  80  BLOCK  800  PEBM  BECFM  FB 
FILE  OUTBOF  DISK  DATA  PLOT  (DISP  MOD  LBECL  508  BLOCK  508  BECFM  F 
PLOT 15 

FILE  INMOVE  DISK  DATA  PLOT  (LRECL  508  BLOCK  508  PEBM  BECFM  F 
FILE  OUTMOVE  TAPI  (LBECL  508  BLOCK  508  PERM  BECFM  F 
-LOOP5 
SBEGTYPE 

-Do  you  have  aore  data  to  eater? 

SEND 

6BEAD  VABS  SI 
SB  =  SI 

SIF  -SB  *  .  SGOTC  -LOOPS 
SIF  SB  =  YES  SGOTO  -DAT A 1 
SIF  SB  =  Y  SGOTO  -DATA  1 
SIF  SR  =  NC  SGOTO  -EXIT 
SIF  S3  =  N  SGOTO  -EXIT 
SIF  SB  =  OFF  SGOTO  -EXIT 
SIF  SR  »  B  SGOTO  —BATA  1 
SIF  SR  *  RESTART  SGOTO  -EATA1 
SGOTC  -LGCP5 
-EXIT 

SIF  SO  =  0  SGOTO  -EXIT2 

CP  CHANGE  PRINTER  ALL  DIST  SD  NOHCLC 

-EXIT2 

SIF  6P  HE  0  STYPE  -There  are  SP  sets  of  multiple  plots  awaiting  output 
ERASE  PLOT  DATA 
CP  SET  MSG  c;; 

SEXII 

SCONTBOL  ALL 


D-ll 


FILE:  DATA 


FORTRAN  A  ITT  ELECTRO-OPTICAI  PRODUCTS  DIVISION 


C 

c 

C  PROGRAM  NAME—  DATA  (TEST  CALCULATIONS- DYNAMIC,  BENDING,  STATIC) 

C  WRITTEN  BY  P.B.  PRIDEAUX  MAY  6,  1580 

C  MODIFIED  MAY  22,  1980 

C 

c 

C  VARIABLES . 

C 

c 

C  COAT  -COATING  DIAMETER 
C  DIAM  -MANDBEL  DIAMETEB 

C  F(  ,  )  -DATA  DISPLAY  VARIABLE...  '****•  INDICATES  THAT  THIS  DATA  WAS 
C  SELECTED  AND  USED  IN  THE  DETERMINATION  OP  EQUATION  CONSTANTS  AND 
C  HAS  WITHIN  2  STD  OF  THE  RESULTING  EQUATION..  •**  •  INDICATES  THAT 

C  THIS  DATA  WAS  SELECTED  IN  THE  DETERMINATION  OF  THE  EQUATION  CONSTANTS 
C  BUT  WAS  NCT  WITHIN  2  STD  CF  THE  EQUATION.  •  ***  INDICATES  THAT  THIS 

C  DATA  HAS  NOT  SELECTED  BUT  HAS  WITHIN  2  STD  CF  THE  RESULTING  EQUATION. 
C  *  '  INDICATES  THAT  THIS  DATA  HAS  NOT  SELECTED  AND  WAS  NOT  WITHIN 

C  2  SID  OF  THE  GENERATED  EQUATICN. 

C  FIBER (  )  — PBEFOBM  ID  NUMBEB  STORED  HERE 
C  FI HR  -FIEEB  DIAMETEB 

C  RIND  — KIND= 1  STATIC  FATIQUE;  KIND=2  DYNAMIC;  KIND=3  BENDING 
C  LSTRESS  -LOG  OF  STBESS  VALUE 
C  OD  -OVERALL  DIAMETER 

C  PBEFBM  (  )  -PREFORM  ID  NUMBER  STORED  HERE 
C  TFAIL  -TIME  OF  FAILURE 

C  TL  -TOTAL  LOAD  AS  BEAD  FROM  INDICATOR  (POUNDS) 

C  S#  — VAHIOUS  SUMMATIONS 
C 

c 

C 

REAL  INT,LTIBE,LTBESS, LTR2S2 , LTIME2 , M, MD,MEAN 
DIMENSION  TIME  (100), D(100)  ,LTIME(100) 

DIMENSION  PPCDUC  (100)  ,  ZZZ&52  ( 130)  ,  ..TIME 2  ( 100) 

DIMENSION  NVL (5) ,NVU(5) , STRESS (100) ,F( 100,5) ,Y  (100) ,LTBESS  (100) 
DIMENSION  STRES2  (100)  ,  DATA  (50)  ,  ALIN  (5)  ,  BLIN  (5)  ,R2LIN  (5) 

DIMENSION  PBEFBM  (3) , FIBER (3) , ST D  (5) , TYPE  (3) ,PLOT  (2,3) ,B£(10) 

DATA  BLANK, AST, FHAIF,SHALF/'  »,'**♦*•,•**  *♦'/ 

DATA  TYPE/'STF1','DYN1','BNDl'/ 

DATA  PLOT/'  STA* , ' TIC  ', 'STAN DASD BEND' ,» ING  •/ 

READ  (4,4Q0,EHD=99)  (PBEFBM  ( I)  ,  1=  1 , 3)  ,  (FIBER  (J)  ,  J=  1 , 3) 

400  FORMAT  (3A4,/,3A4) 

WRITE  (2,700)  (PBEFBM  (J)  ,J=1 , 3)  ,  (FIBER  (K)  ,K=  1 , 3) 

READ  (4,410,END=99)  ANS 
KI N  D=0 

IF  (ANS.EQ.TYPE  (1)  )  KIND*1 
IF  (ANS.EQ.TYPE  (2)  )  KIND=2 
IF  (ANS.EQ.TYPE  (3)  )  KIND=3 
IF  (XIND.EQ.O)  GO  TO  90 
410  FORMAT  (1A4) 

READ  (4 ,420, END*99)  (BE  (I)  ,1*1 ,  10) 

420  FORMAT  (10A4) 

WHITE  (7,700)  (PREFRM(J)  ,J=1,3)  ,  (FIBER  (K)  ,K=1,3) 

WRITE  (8,700)  (PBEFBM  (J)  ,  J=  1  ,  3)  ,  (FIBER  (K)  ,  K=  1 , 3) 
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FILS:  DATA 


FOBTB  N  A 


ITT  ELECT HO-QPTICAI  PBODUCTS  DIVISION 


700  FOBNAT  (2X,'PBEFOBN  NO.  • , 3A4 , 1  OX , • FIBEH  NO.  • ,3A4) 
iBITE  (7,710)  (FL0T(I,  KIND)  ,1=1,2)  ,  (BE(J),J=1, 10) 

BBITE  (3,710)  (PLCT(I, KIND)  ,1=1,2)  ,  (HE(J),J=1, 10) 

710  FOBNAT  (16X,2A4,»  TEST  DATA* ,/, 2QX, 1 0A4,/) 

NT=0 

GO  TO  (10,50,60) ,  KIND 
10  iBITE  (7,720) 

720  FOBNAT  (////, 9X, • FD • ,71, »OE •, 8X ,' ND •, 1 21 ,» STBESS* , 5X, » TINE  TO  FAIL 
COBE  ' ,8X,  '  (NIL)  • , 4X, *  (NIL)  ',5X,»  (NIL)  *,9X,»  (G  PASCAL)  ',6X,*  (SEC) 
C») 

K=0 

L=C 

S1=0 

S2=0 

S3-Q 

S4=0 

S5=0 

S6=0 

S7=0 

sa=o 

S9=0 
S10=0 
S1 1=0 
SI  2=0 

BEADING  IN  THIS  OB CEB  TINE  TO  FAILUEE,  FIBEB  DIANETEB, 

COATING  DIANETEB  AND  NAHDBEL  DIANETEB 

DO  26  1=1,100 

BEAD  (4 , *, EBB=9 2, END=28)  TF AIL , EIBB , COAT, Cl AN 
LTIHE  (I)=A10G10  (TFAI1) 

CALCULATING  IN  GIGA  PASCALS 

STBESS  (I)  =  (  (7.  2E04*0.5*FIBB)/  (0.5*  (CO  AT*  El  AN)  ))/1000.  0 
LTBESS  (I)  =ALOG  10  (SIBESS  (I)  ) 

BBITE  (7,730)  FIBB, COAT, DI AN, STBESS (I) ,TEAIL 
730  FOBNAT  (9X , F4. 2, 5X,F4. 1 , 5X ,F5. 1 , 5X , 5X , F6.4, 3X, F7. 0) 

BBITE  (2,210)  LTBESS (I) , LTINE  (I) , SIBESS  (I) ,TF AIL 
210  FOBNAT  (5X , E12. 6, 5X, El  2. 6, 5X, E 1 2. 6 , 5 X, E 12. 6) 

THESE  EQUATIONS  GENEBATE  THE  SUNNATIONS  OF  THE  VABIABLES  FOB  USE  IN 
CONPUTING  THE  FOUB  BEGBESSIONS. 

51  =  SI  +  STBESS  (I) 

52  *  S2  ♦  TFAIL 

53  =  S3  ♦  SIBESS  (I) **2 

54  =  S4  ♦  TFAIL**2 

55  -  S5  ♦  STBESS  (I) *TFAIL 

S12*  S12*  (LOG  (STBESS  (I)  ))  *LCG  (TFAIL) 

IF  (STBESS  (I))  20,20,12 
12  S7  ■*  S7  ♦  LOG  (STBESS  (I)) 

S9  =  S9  ♦  IFAIl*L0G  (STBESS  (I) ) 

S10*  S10*  (LOG  (STBESS  (I)  ) )  +*2 
14  IF  (TFAIL)  22,22,  16 
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16  S6  *  S6  ♦  LOG  (TF AIL) 

sa  =  S3  ♦  STRESS  (I) *LQG  (TFAIL) 
S11=S1U  (LOG  (TFAIL)  )  **2 
GOTO  24 
20  K  -  1 
GOTO  14 
22  L  *  1 
24  NT=*NT+ 1 
26  CONTINUE 
28  CONTINUE 


OK.  HEBE  3E  GO  INTO  THE  CALCULATIONS  OF  THE  CUBVE  FITTING  CONSTANTS. 

LETS  USE  THE  CONSTANTS  ALIN,  BLIN,  B2LIN  FOB  THE  LINEAB  FIT.  LIKE- 
HISE  LETS  USE  AEXP,  BEXP,  B2EX P ,  ALU ,  3LN,  B2LN ,  APSiB  ,  BPMB  AND 
B2PHT  FOB  THE  CONSTANTS  ASSOCIATED  WITH  EXPONENTIAL,  NATUBAL  LOG  AND 
PC  WE  E  CUBVE  FITS. 

30  BLIN  (1)  =  (S5-S1*S2/NT)/  (S3-S 1  **2/NT) 

ALIN  (1)  =S2/NT-ELIN  (1)  *S1/NT 

B2LIN  ( 1)  =  (  (S5— S 1  *S2/NT)  **2)  /  (  (S3-S1**2/NT)  *  (S4-S2**2/NT)  ) 

I=A*EX 


IF  (I)  31,31,32 

3 1  EEXP=  (S8-S  1*S6/NT) /  (S3-S1**2/NT) 

AEXP=EX?  (S6/NT-BEXP*S1/NT) 

B2EXP=  (  (S8-S1*S6/NT)  **2) /  {  (S3-S  1  **2/NT)  *  (S  1 1-S6»*2/NT)  ) 
Y=  A  EXP  (EX) 


32  IF  (K)  33,33,35 

33  BLN=  (S9-S7*S2/NT) /  (S  10-S7**2/NT) 

ALN=S2/NT-BLN*S7/NT 

R2LN=  (  (S9-S2*S7/NT) **2) / ( (SI 0-S7**2/NT) *  (S4-S2**2/NT) ) 
I*  A*B  LN  ( X) 


IF  (L)  34,34,35 

34  BPHB-  (S12-S7*S6/NT)/(S10-S7**2/NT) 

APBB-EXP  (S6/NT-BPBF*S7/NT) 

B2PHB*  ( (S12-S7*S6/NT) **2)/  (  (S 10-S7**2/NT) *  (S11-S6**2/NT) ) 
ALFHE= ALOG 10  (APHB) 

T*A  X**E 


35  BBITE  (8,810)  ALIN(1),  BLIN(1),  B2LIN(1) 

810  FOBNAT  (59X, 'COEFFICIENT*  ,/ ,  63X,  '  OF  '  ,/,58X,  '  DETEBfilNATION '//, 

C' LINEAB ',  1 2X, ' T  *  ',G10.4,'  ♦  ',G10.4,'  S' , 1 3X, P6. 3 ,//) 

IF  (L)  36,36,37 

36  H3ITE  (8,820)  AEXP,  BEXP,  B2EXP 
820  FOBSAT  ('EXPONENTIAL',  7X,*T  =  ',G10.4,'  EXP  (  ',G10.4,'  S)  • , 

C11X,F6.3,//) 

37  IF  (K)  38,38,99  i 
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38  HRITE  (8, 330}  ALN,  BLN ,  R2LN 

830  FOBMAT  ('LOGARITHMIC' ,  7X,'T  =  *,310.4,'  ♦  *,G10.4,»  LN  S', 

C13X,  F6.3,//) 

IF  (L)  39,39,99 

39  HRITE  (8,840)  BPWB  ,  ALPWR,  B2PWE 

840  FOB NAT  ( *  LOG-LCG ' , 7X , ' LOG  1  =  ',G10„4,'  S  +',2X,G10.4, 

C  15X,F6. 3,//) 

TF3=  (10.0**  (  (3.0-ALPBR)/BPWR)  ) 

TF4=  ( 10. 0** ( (4. 0-ALPNB) /BP HR)  ) 

VN=AES  (3PHB) 

iBITING  CONTBOL  FILE  FOB  PLOT15  (PHP'S  GENERAL  DATA  PLOTTING  SOUTINE) 

WRITE  (8,850)  TF4,TF3,VN 
850  FORMAT  (//, 'USING  LCG-LCG  EQUATION • , //, 

C 'STRESS  AT  TIME  TO  FAILUBE  OF  10,000  SEC.  =  ',F7.4,»  GIGA  PASCAL', 
C//, 'STRESS  AT  TIME  TO  FAILURE  OF  1,000  SEC.  =  ',F7.4, 

C'  GIGA  PASCAL'  ,/*/,33X,  *«  VALUE  *  ',F7..4) 

HRITE(1,  100)  NT,  (PBEFBM(J)  ,J=1,3)  ,  (FI3EB  (K)  ,K=1,3) 

100  FORMAT  ('01, 100142'  ,/,  •  2  ',115,/,'  5.50  7. 50 • ,/, ' LNLG' ,/, 

C' STRESS  (G  PASCAL) ',/, 'TIME  TC  FAILURE  (SEC)',/, 

C'STATIC  FATIQUE  E ATA ',/,' PREFORM- ', 3 A4 , '  FIBER-* , 3A4,/, 

C' 1-42, 8. 0,0- 2, 7.65',/, • (39 X, E 1 2. 6, 5X,E 12.6) ') 

HRITE  (1,110)  NT,  (PBEFBM(J)  ,J=1,3)  ,  (FIBER  (K)  ,K=1,3)  ,TF4,TF3,VN 
110  FORMAT  ('01,100146',/,'  2  ',1l5,/,»  5.50  7. 50* , /, ' LGLG • ,/, 

C' STRESS  (G  PASCAL) ',/, 'TIME  TC  FAILUBE  (SEC)',/, 

C'STATI-C  FATIQUE  DATA ',/,' PREFORM-' ,  3  A4 ,  •  FIBER-' , 3A4,/, 

C'PAILORS  STRESS  AT',/,'10000  S EC= • , P5. 3,/,  '  1000  SEC=* ,F5. 3,/, 
C'N=',F 5.  1,/,  *1-42,8-0,0.2,7.  65,0.  16,1.3,0.3,  1-0, 0.3, 0.8,  1.28,0.6* 
C,/,'  (39X,E12.6,5X,E12.  6)  •) 

GO  TO  99 
50  HRITE  (7,740) 

740  FOBMAT  (//,10X,'  LOAD  ',5X,'  DIAMETER  ',5X,'  STRESS  • ,/, 

C  10X,'  (LBS)  ' , 5X, '  (MILS)  ',5X,'  (G  PASCAL)',/) 

DO  52  1-1,100 

BEAD (4 ,*, EBR=92, ENE=64)  FORCE, FIBR 
CALCULATE  STRESS  FROM  LOAD  &  DIAMETER 
FORCE-FORCE/2.0 

STRESS  (I)  =  ((14.0*FCRCE)/(11.  0*  (FIBR**2))  )*6.895 
LTBESS  (I)  =  ALOG  10  (STRESS  (I)  ) 

HRITE  (7,750)  FORCE, FIEB, STRESS  (I) 

750  FORMAT  (13X,F5- 2,9X,F5- 2,9X,F6.3) 

F  (I,  1)  =  AST 
F  (1,2) =AST 
F  (1,3)  =AST 
F  (1 , 4)  =AST 
F(I,5)=AST 
NT=NI* 1 
52  CONTINUE 
GO  TC  64 
60  HRITE  (7,760) 

760  FOBMAT  (//,'  DIAL  BENDING  COATING  FIBER  ', 

C*  STRESS ' ,/, '  READING  DIAMETER  DIAMETER  DIAMETER  ', 
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C*G  PASCAL*,/) 

DO  62  1=1,100 

BEAD  (4,*,E2B=92,  END=64)  HEAD, COATD, FIBB 
CALCULATE  STRESS  FHCN  LOAD  S  DIAMET2H 
BAD=B£  AD-562.  5 

STBESS  {I)  =  ( (1.  05E*04*FIBB)  /(BAD-CGATD)  )  *0.006895 
LTBESS  (I)-ALOGIO  (STRESS  (I)  ) 

SHITE  (7,770)  BEAD, BAD, CCATD , FIBB, STRESS  (I) 

770  FOHBAT  (5X,  F6.2 , 5X ,  F6.  2 , 5X,  F6.  2 , 5X,  F6.  2,5X,  F6. 3) 

F  (1 , 1)  =  AST 
F  (I, 2)  -AST 
F (I, 3) = AST 
F (I, 4) = AST 
F(I,5)=AST 
NT=NI+1 
62  CONTINUE 

64  N=NT 

ABBANGE  STBESS  VALUES  IN  DECBEAS2NG  ORDEB 

65  H=0 

DO  66  1=2, NT 

IF  (STBESS  (I).  LE.STBESS  (1-1)  )  GO  TO  66 
T=  STBESS  (1-1) 

STBESS  (I—  1 )  =STRESS  (I) 

STBESS  (I)  =T 
B=LIBESS  (I-  1) 

LTBESS  (I-1)=LTHESS  (I) 

LTBESS  (I)=H 

n=i 

66  CONTINUE 

IF  (H.  NE-0)  GO  TO  65 
JT=  5 

IF  (NT.LE.  6)  JT=  1 

NVL  (1)  =  1 

NVL  (2)  =1 

NVL  (3)  =1 

NVL  (4)=HT/3 

NVL  (5)  =  (NT*2)/3 

NVU  (1)  =NI 

NVU  (2) -  (NT*2)/3 

NVU  (3)  =NT/3 

NVU  (4)  =NT 

NVU  (5)  =NT 

I  (1)=  (LOG  (LOG  (1/(0.0001)))  ) 

DO  78  J= 1 , JI 

IF  (NT.EQ.O)  GO  TO  78 

X=0 

K=  0 

S1=0 

S2  =  0 

S3=0 

S4=0 
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S5=0 

NV=UNVU  (J)— NVL  (J) 

H?L0=HVL(J) 

NV00=NVU  (J) 

8L-NVL0 

ANI=(1.0/NV) 

ASSIGN  PROBABILITY  OF  FAILURE  GIVEN  THE  NO.  CF  DATA  POINTS 
DO  68  1=2, NV 

I  (I)  =  (LOG  (LOG  (1/  (  (  (1-1)  *ANT)  )  ) )  ) 

68  CONTINUE 

DO  70  1=1, NV 

51  =  SI  ♦  LTRESS  (NL) 

52  =  S2  *1(1) 

53  =  S3  ♦  LTBESS  (NL)  **2 

54  =  S4  ♦  Y  (I)  **2 

55  =  S5  ♦  LTBESS  (NL) *Y  (T) 

NL=NL* 1 

70  CONTINUE 

Y=A+B*LOG  (X) 

BLIN(J)  =  (S5-(S1*S2/NV)  )/(S3-  (S1**2/NV)  ) 

ALIN  (J)=  (S2/NV)-  (BLIN  (J)  ♦S1/NV) 

R2LIN  (J)  =  (  (S5-  (S1  +  S2/NV)  )  **2 )/ (  (S3-  (S1**2/NV) )  ♦  (S4-  (S2**2/NV)  )  ) 

STANDARD  ERROR  OF  ESTIMATE 

NL=NVL0 
DO  74  1=1, NV 

X=X*  (  (LTRESS  (NL)-  (  (I  (I)  -ALIN  (J)  )  /BLIN  (J)  ) )  **2) 

NL=NL+1 

74  CONIINUE 

IF  (IV.  EQ.  2.0)  GO  TO  7  8 
STD  (J)  =3QHT  (X/  (NV-2)  ) 

ST AN=  (2.  0*  (STD  (J)  )  ) 

ANT= (1.0/NT) 

X* AES  (LTRESS  (1)  -((1(1)  -ALIN  (J) )  /BLIN  (J)  )  ) 

IF  (X.  LB.STAN)  GO  TO  75 

F(1,J)=FHALF 

IF  (1.EQ.NVL0)  GO  TO  76 

F(1,J)=BLANK 

GO  TC  76 

75  IF  (1.EQ.NVL0)  GO  TO  76 
F(1,J)=SHALF 

76  DO  78  1=2, NT 

T(I)=(LOG(LOG(1/((  (1-1)  *ANT)  )))  ) 

X=  AES  (LTBESS  (I)  -  (  (Y  (I)  -ALIN  (J)  )  /BLIN  (J)  )  ) 

IF  (X.LE.STAN)  GO  TO  77 
F  (I,  J)  =FHALF 

IF  (I.LE.NVUO.  AND.I.GE.NVLO)  GO  10  78 
F  (I, J) = BLANK 
GO  TC  73 

77  IF  (I.LE.NVUO.  AND.I.GE.NVLO)  GO  TO  78 


m 


■  i 
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EJ 


F  (I,  J) =SHALF 
78  CONTINUE 
NT=N 

PRINT  CALCULATED  VALUES 

WRITE  (8,890)  (ALIN  (J)  ,  J=1  ,5)  ,  (5LIN  (J)  ,  J=1,5)  ,  (R2LIN  (J)  ,  J=  1,5)  , 

C  (STB  (J)  ,J=1,5) 

890  FORMAT  (‘LOG  (LOG  (1/ (1-F) ) )  =  A  ♦  H*  (LOG  (STRESS)  )*,//,  4X, 

C  ‘INTERCEPT  -  A  -  • , 5  (E  10. 4, 2X)  ,/, 8X , • SLOPE  -  M  =  • , 5  (E 10. 4 , 21) , 

C  /, 5X, ‘ COEFF  OE  DET  =  •  ,5  (E 1  0.  4 , 2X)  ,/, 

C14X,  ‘SID  =  • ,5  (E10.4,2X) ,//, 

C  12X, ‘STRESS  VALUES  MITHIH  2.0  STD  OF  ABOVE  EQUATION*) 

A=0. 0 

DO  80  1=1, NT 

WRITE  (8,898)  I  (I)  , STRESS  (I)  ,  (F  (I ,  J)  ,  J=  1 , 5) 

898  FORMAT  (1  X,  F7.3 ,3X , F7.  3 ,4X,  5  ( 2X ,  1 A4 , 6X)  ) 

WRITE  (2,220)  STRESS  (I) ,X (I) 

220  FORMAT  (2X, E  12.  6 , 2X, E 1 2.  6) 

A= A ♦STRESS  (I) 

80  CONTINUE 
A=  A/NT 

SX=STBE5S  ( 1) 

SN=  STRESS  (NT) 

NT=NT+ 1 

WRITING  CONTROL  FILE  FOR  PLOT  15  (PH?‘S  GENERAL  DATA  PLOTTING  ROUTINE) 

BRITE  ( 1 ,  120)  NT,  (BE  (I)  ,  1=  1 ,  10)  ,  ( PLOT  (J,  KIND)  ,J=1,2)  , 

C  (PHEFBM  (K)  ,  K=1 ,3)  ,  (FIBER  (L)  ,  L=1, 3)  ,SX,SN,  A,  BLIN  (1) 

120  FORMAT  (‘01, 100147*  ,/,  *  2  ‘,  115,/,*  5.50  7.50* ,/, ‘LGLN* ,/, 
C*STBESS  (G  PASCAL)  *,/, ‘LOG  (LOG  (1/1-F) )»,/, 10A4,/, 

C ‘DYNAMIC  TEST  DATA  • 2A4,/, • PBEFORM- • ,3 A4, •  FIBER-*, 3A4,/, 

C  *  MAX= ‘ , F7. 3, •  G  PA‘,/,*MIN=‘,F7.3,‘  G  PA‘,/, 

C  •  AVG=‘,F7.3,»  G  PA»,/,‘  M=*,F7.3,/, 

C‘0. 0,0. 0,0. 86, 8. 0,0. 2, 7. 65,0. 0, 0. 0, 0. 0, 0. 0 , 0. 0 ,0. 0, 0- 0 ,0. 0* ,/, 

C»  (2X,E12.  6,2X,E12.  6)  •) 

GO  IC  99 

SIMPLE  ERROR  HANDLING  ROUTINES 

90  WRITE  (7,790) 

WRITE  (8,790) 

790  FORMAT  (//, *  UNKNOWN  TYPE  OF  DATA  FILE  TYPE  ON  LINE  3  OF  DATA  FILE*, 
C/,  ‘DATA  TYPE  MUST  BE  “STF1“,  “DYN1“,  OB  *  *  BND 1  •  *  ') 

GO  TC  99 

92  BRITE  (7,792)  NT 
WRITE  (8,792)  NT 

792  FORMAT  (//, 'IMPROPER  DATA  IN  THE  ',113,  •  BECORD  OF  THE  DATA  FILE') 
99  STOP 
END 


D-18 


FILE:  STATIC  DATA 


A  ITT  ELECTHO-OPTICAL  PRODUCTS  DI7ISI0N 


ES-040489B  * 

800123-1 

STP1 

STAHDABD  STATIC  PATIQOE  TEST  (EIAHPLE) 


445.  0 

4.90 

20.8 

72.0 

385.0 

5.00 

20.6 

72.0 

410.0 

4.90 

20.4 

72.0 

450.0 

4.90 

20.  8 

72.0 

390.0 

5.00 

20.8 

72.0 

6000.0 

5.10 

21.2 

83.0 

1075  20.0 

5.00 

20.0 

83.0 

600.0 

5.90 

20.6 

83.0 

40440.0 

4.90 

19.8 

83.0 

6120.0 

5.30 

20.4 

83.0 

593460.0 

5.10 

20.8 

95.0 

402600.0 

4.70 

20.2 

95.0 

239640.0 

5.30 

20.8 

95.0 

163200. 0 

5.50 

21.4 

95.0 

LURE  (SEC) 


STATIC  FRT ! CUE  GRTR 


STRESS  (G  PfiSCRU 


STATIC  FATIQUE  DATA 
PAEFQRM-ES-040489B  FIBER-  300 123-1 
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HIGH  STRENGTH,  LONG  LENGTH 
OPTICAL  FIBERS 
PROCESS  SPECIFICATION 
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1.0  SCOPE 

The  process  to  be  defined  is  the  fabrication  of  a  high  strength, 
long  length  optical  fiber  with  a  compressive  cladding.  This 
specification  will  document  the  fabrication  process  in  suffi¬ 
cient  detail  in  order  that  anyone  with  similar  equipment  and 
knowledge  in  the  art  of  optical  fiber  fabrication  can  duplicate 
these  techniques. 


2.0.  APPLICABLE  DOCUMENTS 

The  Technical  Proposal  No  R77-043  submitted  by  ITT  Electro- 
Optical  Products  Division  on  November  1,  1977  and  the  Eighth 
Bimonthly  Progress  Report  written  under  this  contract 
(contract  N00123-79-C-0301)  are  the  only  two  documents 
referenced  to  in  this  text. 


3.0  REQUIREMENTS 

The  following  are  the  minimum  requirements  necessary  to  produce 
long  length  optical  fiber  with  a  highly  compressed  surface 
layer.  The  surface  layer,  which  was  produced  under  this 
contract,  increases  the  durability  of  a  fiber. 

3 . 1  Equipment 

3.1.1  Glass  Cleaning  Station 

This  station  should  be  capable  of  properly  handling  hydro¬ 
fluoric  and  nitric  acids  (liquid  and  fumes)  which  are  present 
during  glass  cleaning.  It  should  be  of  adequate  area  to  enable 
the  cleaning  of  approximately  three  foot  silica  tubes  while 
providing  a  safe  environment  for  the  technician. 

3.1.2  Glassworking  Lathe 

This  lathe  must  be  capable  of  supporting  glass  sizes  of  up 
to  37  mm  in  diameter.  Support  should  consist  of  chucks  on 
both  head  and  tail  stocks  and  these  units  must  rotate  up 
to  150  rpm  while  maintaining  the  set  speed  within  ±2%. 

The  fire  carriage  should  be  capable  of  traversing  at  least 
70  cm  between  the  head  and  tail  stocks.  The  fire  carriage 


must  maintain  a  smooth  traverse  over  a  range  of  speeds 
from  1  to  20  cm/min  within  ±1%  of  the  set  value.  A  switching 
system  must  be  present  to  change  the  fire  carriage  direction, 
while  the  return  speed  should  be  fast  enough  to  prevent  delays 
and  damage  to  the  preforms.  A  traverse  delay  at  the  intake 
and  exhaust  ends  is  also  required.  The  delay  should  be  capable 
of  being  varied  up  to  10  s  depending  upon  process  conditions. 

The  fire  carriage  will  be  equipped  with  an  oxyhydrogen  torch 
which  has  an  adjustable  height  as  well  as  individual  torch  adjust 
ments  for  the  hydrogen  and  oxygen  so  that  a  proper  deposition 
temperature  can  be  attained.  An  optical  pyrometer  should  be  use 
to  continuously  monitor  the  deposition  temperature  throughout 
the  process.  Oxyhydrogen  hand  torches  must  be  available  at 
both  head  and  tail  stocks  to  aid  in  preform  preparation. 

3.1.3  Chemical  Delivery  System 

This  system  will  be  used  to  deliver  gaseous  chemicals  to  the 
deposition  torch  and/or  substrate  tube  on  the  glass  fabrication 
lathe.  The  unit  consists  of  a  constant  temperature  enclosure 
capable  of  maintaining  ±0.5°C  while  the  delivery  system  is  in 


use.  Contained  in  this  enclosure  are  the  chemical  bubbler  for 
the  SiCl^,  T^Cl^,  and  POCl^.  The  bubbler  should  be  equipped 
to  deliver  a  controlled  flow  of  oxygen  through  the  bubblers 
to  the  process  line.  Chemicals  such  as  oxygen,  helium,  and 
BClj  are  delivered  through  gas  controllers  and  into  the  process 
lines.  All  gas  flows  must  be  controlled  to  ±21  of  designated 
flow  as  specified  in  Tables  3-1,  3-2,  and  3-3. 

Each  of  these  chemical  vapors  and  gases  is  valved;  this  allows 
only  the  specific  chemicals  called  for  to  be  delivered  to  the 
process  at  any  given  time.  The  process  line  consists  of  a 
transfer  line  from  the  chemical  delivery  system  to  the  glass¬ 
working  lathe. 

The  process  line  is  connected  to  either  the  deposition  torch 
or  to  the  intake  end  of  the  tube  assembly  via  a  rotary  coupling 
and  a  filter. 

This  entire  system  should  be  free  from  leaks,  not  permeable 
to  moisture,  and  free  from  metal  impurities  which  will  adversely 
affect  the  quality  of  the  deposited  glass.  The  system  should 
be  purged  with  dry  nitrogen  when  not  in  use. 


Table  3-1.  Preform  Fabrication  Parameters  for  the  Standard  Three- Layered,  Internally 
Deposited  Compressive  Cladded  Fiber. 
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Table  3-2.  Preform  Fabrication  Parameters  for  the  High  Phosphorous  Three-Layered 
Internally  Deposited  Compressive  Cladded  Fiber. 


of  chemical 


Table  3-3.  Preform  Fabrication  Parameters  for  the  Typical  Four-Layered,  Externally 
Deposited,  Compressive  Layer  Fiber. 
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3.1.4  Exhaust  System 

An  exhaust  system  on  or  about  the  lathe  is  necessary  in  order 
to  properly  remove  the  excess  heat  and  spent  process  chemical 
which  develops  during  the  fabrication  of  the  preform.  Proper 
handling  and/or  treatment  is  necessary  to  keep  the  amount  of 
exhausted  doped  silica  soot  and  chlorine  gases  below  the  threshold 
limit  value  (TLV)  set  by  the  American  Council  of  Government 
Industrial  Hygienists  (ACGIH) . 

3.1.5  Clean  Area  for  Deposition  Lathe 

The  area  in  which  the  glassworking  lathe  will  be  located 
should  be  kept  as  clean  and  dust  free  as  possible  in  order 
to  prevent  degradation  of  the  preform  surface  which  may  adversely 
affect  the  fiber  strength  and  durability. 

3.1.6  Chemical  Etching  Station 

The  station  consists  of  an  exhaust  hood  with  a  lamp  and  an 
etching  tank.  The  facility  must  be  capable  of  handling 
hydrofluoric  acid.  The  etching  tank  provides  a  constant 
flow  of  acid  over  and  around  the  immersed  preforms.  Details  of 
the  setup  are  shown  in  Figure  3*1. 


m  u:«i 


Figure  3-1.  Chemical  Etching  Facilities 


Preforms  shorter  than  23  cm  were  usually  etched  in  a  magnetically 
stirred  etching  bath.  The  etching  bath  could  simultaneously 
etch  nine  preforms  at  various  temperatures  and  stirring  rates. 
Details  of  the  etch  bath  are  shown  in  Figure  3-2. 

3.1.7  Fiber  Draw  Tower 

The  equipment  necessary  to  draw  a  preform  into  a  high  strength, 
long  length  optical  fiber  consists  of  the  following  components: 

a.  Down-feed  system 

b.  Hj/Oj  torc^es 

c.  Diameter  control  device 

d.  Primary  coating  applicator 
€.  Curing  oven 

f.  Drive  capstan 

g.  Extruder 

h.  Proof  tester 

i.  Spooler 


A  basic  diagram  of  the  draw  tower  is  presented  in  Figure  3-3 
which  shows  the  equipment  and  the  sequence  in  which  it  is  used 
during  the  fabrication  process. 


Figure  3-2.  Chemical  Etching  Setup  for  Short  Preforms. 


PREFORM  FEEO 


Figure  3-3.  Optical  Fiber  Drawing  Instrument  Diagram. 


3. 1.7.1  Down-Feed  System 

This  unit  attaches  to  the  preform  at  the  butt  end  and  lowers 
the  preform  into  the  flame  at  a  constant  rate  to  supply 
glass  for  the  draw  process.  The  feed  rate  should  be  between 
0.3  and  0.8  cm/min  with  a  linearity  of  ±0.5%. 

3. 1.7. 2  Hydrogen  Oxygen  Torches 

These  torches  must  be  mounted  in  the  tower  in  such  a  way 
as  to  allow  for  x,  y,  and  z  positioning.  The  torches  must 
also  be  capable  of  obtaining  temperatures  of  2000°C. 

3. 1.7. 3  Diameter  Control  Device 

This  unit  will  measure  the  diameter  of  the  fiber  as  it  exits 
the  flames.  The  unit  supplies  fiber  diameter  information 
via  a  feedback  loop  to  the  drive  capstan.  By  precisely 
controlling  the  draw  speed,  the  proper  diameter  of  fiber  is 
maintained.  This  entire  system  must  be  able  to  produce 
fibers  of  diameters  of  125  ym  ±4.0%.  The  detection  window 
in  this  device  should  be  sufficiently  large  to  allow  for  some 


E-15 


lateral  motion  of  the  fiber  during  draw  (2  mm  x  2  mm).  The 
diameter  sampling  rate  should  be  1000  samples/s  to  allow  for 
proper  diameter  measurement  and  control. 

3. 1.7. 4  Primary  Coating  Applicator 

The  Applicator  is  a  dipcoater  in  which  the  silicone  RTV  elastomer 
will  be  applied  to  provide  a  coating  on  the  fiber.  The  applicator 
must  be  mounted  in  the  tower  to  allow  for  on-line  dipcoating 
of  the  fiber.  The  fiber  must  be  thoroughly  encapsulated  by 
the  silicone  elastomer  to  an  outer  diameter  of  12.5  mils 
±2  mils  (317  ym  ±50  ym) .  After  application  the  silicone 
elastomer  is  thermally  cured. 

3 . 1 . 7 . 5  Curing  Oven 

The  primary  coating  thermal  curing  oven  should  be  capable  of 
maintaining  a  temperature  of  500°C  ±5°C.  The  curing  oven 
should  be  located  directly  below  the  primary  coating  appli¬ 
cator  to  assure  a  uniform,  well-cured  coating.  The  curing 
oven  should  be  37"  (94  cm)  long  for  draw  speeds  up  to  32  m/min. 

3. 1.7. 6  Drive  Capstan 

The  drive  capstan  provides  the  pulling  force  during  fiber 
draw.  It  must  be  capable  of  maintaining  no-slip  contact 
with  the  fiber  at  all  times  without  degrading  the  coating 


surface.  The  capstan  must  be  able  to  attain  draw  speeds  of 
32  m/min  with  sufficient  sensitivity  and  control  to  respond 
to  the  diameter  monitoring  device.  The  capstan  should  have 
both  automatic  and  manual  control  capability. 

3 . 1 . 7 . 7  Extruder 

The  extruder  should  develop  sufficient  pressure  and  tempera¬ 
ture  to  extrude  polyester  elastomer  materials  on  the  optical 
fiber  at  diameters  up  to  20  mils  (308  urn) .  The  coating 
thickness  is  a  nominal  4  mils  (102  rm)  over  the  primary  coating 
The  extruder  must  be  capable  of  applying  the  polyester  elastome 
in-line  at  a  draw  speed  up  to  32  m/min.  The  extruder  must 
also  be  able  to  maintain  the  temperature  specified  in  Table  3-4 

3. 1.7. 8  Proof  Tester 

The  proof  tester  should  be  capable  of  providing  an  in-line 

2 

proof  stress  level  of  690  N/mm  (100  kpsi)  on  the  fiber. 

3. 1.7. 9  Spooler 

The  spooler  should  be  capable  of  taking  up  fiber  onto  a  spool 
at  a  rate  equal  to  the  draw  speed  and  at  low  tension  (50  grams 


Table  3-4.  Extruder  Setup  and  Draw  Parameters. 


Temperature  in:  neutral 

Temperature  out:  °C 

Zone  1:  200°C 

Zone  2:  230°C 

Zone  3:  232°C 

Adapter:  235°C 

Die:  238°C 

TIP  inner  diameter: 

70  mils 

TIP  outer  diameter: 

100  mils 

Die  inner  diameter: 

230  mils 


Screens:  stainless  35  mesh 
stainless  60  mesh 
stainless  80  mesh 
Material  drier:  80°C 
Line  speed:  32  m/min 
Screw  speed:  1.5  rpm 
Screw  set:  14 

Material  use:  1.5  lb/h 

Motor  amperage:  2.0 

Melt  temperature:  424°F 

Final  outer  diameter  of 
fiber:  20.5  mils 


3.1.8  Clean  Room  Facility  for  Fiber  Draw  Tower 
The  entire  draw  system  should  be  kept  as  clean  and  dust  free 
as  possible.  This  prevents  degradation  in  the  strength  of 
the  fibers. 

3.2  Materials 

The  following  is  a  list  of  materials  used  in  the  fabrica¬ 
tion  of  high  strength,  long  length  optical  fibers: 

a.  Natural  fused  quartz  tubing:  25-mm  x  21.6-mm 
waveguide  grade 

b.  Silicon  tetrachloride:  electronic  grade  -  resistivity 
to  be  >200  ohm* cm 

c.  Germanium  tetrachloride:  electronic  grade 

d.  Phosphorus  oxychloride:  99.99%  purity 

e.  Boron  trichloride:  99.95%  purity 

f.  Titanium  tetrachloride:  99.95%  purity 

g.  Oxygen:  99.99%  purity 

h.  Helium:  99.995%  purity 

i.  Hydrogen:  bulk  tank  grade 

j.  Oxygen:  bulk  tank  grade 


k.  Nitrogen:  bulk  tank  grade 

l.  Natural  fused  quartz  rod:  12 -mm  outer  diameter  (od) 

m.  Natural  fused  quartz  tube:  36  ram  x  32  mm 

n.  Natural  fused  quartz  tube:  25  mm  x  21.6  mm 

o.  Hydrofluoric  acid:  reagent  grade 

p.  Nitric  acid:  reagent  grade 

q.  Deionized  water 

r.  Filtered  dry  air 

s.  Paraffin  wax 

t.  Sylgard®184  silicone  RTV 

.  u.  Hytrel ^  7246  polyester  elastomer 

3.3  Required  Procedure  and  Operations 

3.3.1  Cleaning  of  Glassware 

The  substrate  tube  must  be  cleaned  to  remove  any  impurities 
from  the  glass  which  can  degrade  both  fiber  strength  and 
attenuation. 


This  operation  is  performed  m  the  glass  cleaning  area  with 
the  necessary  protective  equipment  being  used  by  the  operator. 
The  substrate  tube  is  placed  in  a  bath  of  50%  hydrofluoric 
acid  or  rinsed  with  50%  hydrofluoric  acid  for  10  minutes 
to  remove  the  outer  layers  of  glass.  Subsequently,  a  deionized 
water  rinse  is  administered.  The  glass  is  then  placed  in  a  40% 
nitric  acid  bath  or  rinsed  with  nitric  acid  for  10  minutes  to 
remove  other  foreign  materials.  Again,  a  deionized  water  rinse 
is  administered  and,  after  this  material  is  blown  dry  with 
filtered  nitrogen.  Proper  precautions  must  be  taken  throughout 
the  whole  process  to  prevent  scratching  or  recontaminating 
the  glass. 


3.3.2  Tube  Mounting 

In  tube  mounting,  both  the  intake  tube  and  substrate  tubes  are 
supported  snugly  in  the  lathe  chucks  in  such  a  way  that 
the  unchucked  ends  are  facing  each  other.  Prior  to  sealing 
the  intake  and  substrate  tubes,  the  two  pieces  are  straightened 
to  ensure  a  proper  joint.  This  process  is  done  by  starting  from 
the  lathe  chuck  and  straightening  towards  the  free  end  of  each  of 
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the  tubes.  Afterwards,  both  ends  are  moved  to  approximately  1 
apart,  heated  to  their  softening  point,  and  then  sealed. 

During  this  process  it  is  important  to  keep  a  large  flow 
(100  sccm/mm)  of  filtered  gas  (helium  or  nitrogen)  moving 
through  the  tube  to  prevent  dust  particles  from  entering 
the  tube  during  the  sealing  operation.  It  is  also  important 
that  the  tube's  diameter  be  maintained  throughout  the 
sealing  operation.  Tube  dimension  can  be  maintained  by  control 
backpressure  during  the  sealing  operation. 

The  exhaust  tube  is  straightened  as  described  above  but  then 
tapered  down  using  a  torch  and  paddle  until  it  matches  the 
substrate  tube  diameter.  It  is  then  joined  to  the  substrate 
using  the  aforementioned  techniques. 

3.3.3  Internal  Deposition 

Internal  deposition  is  started  immediately  following  a  fire 
polishing  pass.  The  proper  chemical  gas  flows  are  set  up 
and  allowed  to  enter  the  substrate  tube  assembly  while  the 


traversing  burner  starts  at  the  intake  end  of  the  tube. 

Depending  on  torch  design,  the  flow  rates  of  the  hydrogen 
and  oxygen  are  established  using  the  deposition  temperatures 
specified.  These  flows  and  temperatures  are  listed  in  Tables  3-1 
through  3-3. 

The  glass  layer  that  is  deposited' should  be  clear  and  free 
of  bubbles.  An  oxyhydrogen  hand  torch  should  be  placed  so 
that  its  light  flame  is  heating  the  sealed  area  on  the  exhaust 
tube.  Also,  a  ribbon  burner  should  be  mounted  along  the 
entire  length  of  the  substrate.  This  assembly  should  be 
used  from  mid-cladding  deposition  to  the  end  of  the  final 
collapse  stage. 

During  deposition  silica  soot  will  build  up  in  the  exhaust 
tube.  This  soot  must  be  removed  to  prevent  preform  contamina¬ 
tion  or  clogging. 


3.3.4  Collapse 

Collapse  is  either  a  three-  or  four-pass  operation  depending 
upon  the  amount  of  open  area  remaining  inside  the  tube  after 
deposition.  The  collapse  is  performed  in  a  traverse  direction 
opposite  to  that  of  deposition. 

The  collapse  takes  place  at  temperatures  of  2250°C  with 
the  deposition  torch  traversing  at  a  rate  of  1.5  cm/min  or 
slower  (depending  on  the  open  area) .  The  lathe  spindle 
rotates  at  137  rpm  and  close  attention  should  be  given  so 
that  the  consolidation  of  the  tube  is  complete  and  bubble 
free.  During  the  collapse  a  white  film  will  develop  on  the 
outside  of  the  preform.  This  film  can  be  removed  with  subsequent 
acid  treatment  or  fire  polishing.  Depending  on  the  require¬ 
ments  of  the  draw  tower  feed  system,  the  preform  should  be 
mounted  to  a  section  of  support  rod  which  is  of  a  compatible 
length  and  diameter.  During  this  mounting  procedure  a 
preform  cane  sample  should  be  generated  to  check  the  preforms* 
properties . 


3.3.5  Etching 

The  amount  of  silica  to  be  removed  from  the  preform  will 

depend  on  the  diameter  of  the  preform  and  the  depth  to  which 

etching  is  to  be  performed.  This  depth  is  usually  determined  by 

examining  photographs  of  cross  sections  of  canes  drawn  from 

the  preform.  The  actual  diameter  of  the  preform  is  carefully 

% 

measured  at  1-cm  intervals  to  determine  the  mean  preform 
diameter.  The  target  diameter,  the  diajneter  after  chemical 
etching,  is  computed  using  the  cane  photograph  and  the  mean 
diameter. 

Before  starting  the  etching  process,  the  tip  of  the  preform 
and  the  joint  between  the  preform  and  the  supporting  rod  is 
covered  with  hot  paraffin  wax.  The  preform  is  immersed  in  the 
etchant  bath  with  acid.  The  frequency  of  preform  inspections 
should  increase  as  etching  time  increases. 

After  etching  is  finished  the  preform  is  removed  and  thoroughly 
washed  and  dried. 


3.3.6  External  Deposition 

The  T08  end  of  the  etched  preform  is  firmly  attached  into  one  of 
the  lathe’s  chucks.  A  piece  of  T08  is  placed  in  the  .other  chuck, 
held  firmly,  and  attached  to  the  preform.  To  avoid  contaminating 
the  preform,  straightening  is  performed  in  the  T08  regions  only. 

The  preform  is  then  fire  polished  with  the  deposition  torch 
mounted  on  the  lathe's  carriage  traverse.  Following  the  fire 
polishing  procedure,  the  compressive  layer  is  externally  deposited 
under  the  parameters  detailed  in  Table  3-3. 

3.3.7  Mounting  and  Aligning  Preform 

The  T08  butt  end  of  the  preform  is  clamped  to  the  downfeed 
device.  The  preform  is  then  aligned  so  that  its  axial  centerline 
does  not  deviate  more  than  40  yra  from  the  center  line  of  the 
drawing  process. 


3.3.8  Preliminary  Fiber  Draw 

The  hydrogen/oxygen  torch  is  ignited  in  a  safe  manner  and  the 
flow  rate  adjusted  to  obtain  the  2000°C  drawing  temperature. 

After  flames  have  been  adjusted,  the  torch  alignment  should  be 
checked  with  preform.  This  is  accomplished  by  lowering  the  tip 
of  the  preform  into  the  flame  and  adjusting  x,  y,  and  z  position 
on  torches  until  alignment  is  complete.  The  draw  process  is 
now  ready  to  commence. 

3.3.9  Fiber  Draw  Start  Up 

Set  the  preform  down- feed  speed  at  0.35  cm/min  and  initiate  the 
feed.  Insert  glass  starting  rod  until  contact  with  the  preform 
is  made.  Slowly  remove  glass  starting  rod  from  the  preform.  The 
rod  should  have  the  small  end  of  the  preform  attached  with  fiber 
trailing.  Break  the  starting  rod  from  the  fiber  and  continue 
pulling  the  fiber  manually.  The  fiber  must  then  be  passed  through 
the  dipcoater  assembly,  through  the  curing  oven,  and  around  a 
sheave.  Adjust  the  x-y  position  of  the  dipcoater  for  a  concentric 
coating.  After  adjusting  the  dipcoater,  fill  the  dipcoater 
assembly  with  coating  resin.  When  the  fiber  is  coated,  it  can 


then  be  passed  through  the  extruder.  The  extruder  should  be  set 
up  per  Table  3-1.  After  the  fiber  has  been  passed  through  the 
extruder  and  the  secondary  coating  applied,  the  fiber  is  then 
placed  around  the  drive  capstan.  The  diameter  of  the  fiber  is 
checked.  Adjust  the  draw  rate  (^30  m/min)  to  obtain  the 
desired  fiber. diameter.  When  the  fiber  diameter  has  been 
established,  the  fiber  can  be  put  on  the  proof  tester  wheel. 
Measure  the  proof  test  tension  with  the  tension  meter  and  adjust 
the  tension  to  obtain  the  desired  proof  test  level.  String  the 
fiber  over  the  pulley  system  on  the  spooler  dancer  arm.  The 
dancer  controls  spooling  tension  and  may  be  adjusted  as 
needed.  Pass  the  fiber  from  the  dancer  onto  the  take-up  spool 
and  switch  the  spooler  to  "auto"  to  control  spooling  tension 
for  the  remainder  of  the  draw. 

3.3.10  Fiber  Draw 

Several  times  during  draw,  check  the  level  of  the  coating 
material  in  the  dipcoater,  the  diameter  of  the  secondary 
coating,  the  proof  test  level,  the  spooler  traverse,  and  the 
spooling  tension. 
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3.4  Recommended  Procedures  and  Operations 

3.4.1  Chemical  Handling 

The  chemicals  used  in  these  processes  are  of  high  purity  and 
steps  should  be  taken  to  ensure  their  quality.  By  placing  in 
line  filters  on  all  gas  cylinders  the  contamination  due  to 
cylinder  changes  will  be  reduced. 

3.4.2  Maintenance  of  Glass  Surface  Quality. 

The  preform  fabrication  station  should  be  located  in  a  clean 
room  to  prevent  contamination  of  the  glass  surface. 

3.4.3  Fiber  Draw  Tower  Room 

It  may  be  advantageous  to  have  the  entire  draw  process  in  a 
clean  room.  This  should  improve  the  strength  characteristics 
of  the  drawn  fiber. 
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4.0  QUALITY  ASSURANCE  PROVISION 

Preform  qualities  are  ultimately  determined  through  the  fiber 
testing  procedures  outlined  in  the  Prelimanary  First  Article 
Inspection  Report  submitted  under  this  contract  on  February 
13,  1980.  However,  some  basic  testing  is  undertaken  prior  to 
and  following  the  fabrication  of  the  preform  to  ensure  good 
quality  at  this  early  stage  of  the  process. 

4.1  Responsibility  for  Inspection 

Inspection  of  the  preforms,  components  of  the  preforms,  and 
vendors'  lot  analyses  (which  are  supplied  with  materials)  is 
the  responsibility  of  the  lathe  operator.  Standard  laboratory 
procedures  are  required  to  inspect  preform  and  glass  components. 
High  purity  laboratory  procedures  are  required  to  determine 
purity  of  materials  and  only  used  when  necessary. 

4.2  Monitoring  Procedures  for  Equipment  Used  in  the  Process 

The  following  is  a  list  of  routine  maintenance  and  inspection 

items  which  should  be  previous  to  every  fabrication  start: 

a.  Fill  SiCl^  bubbler  to  a  predetermined  level  prior 
to  deposition  of  each  run 
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b.  Time  traverse  system  and  lathe  spindle  rpm 

and  adjust  if  necessary  to  ensure  reproducibility 

c.  Monitor  chemical  bubbler  enclosure  temperature 
and  adjust  to  nominal  figure 

d.  Clean  rotary  coupling  filter  and  traversing 
burner  torch  tips  to  remove  foreign  matter 

e.  Check  exhaust  system  to  ensure  that  it  is 
functional 

During  chemical  etching  the  preform  is  routinely  inspected 
for  visible  defects.  Stress  concentrations  are  detected  by 
using  polarized  light.  If  the  stress  level  becomes  too  acute 
the  area  of  high  stress  is  cut  out  of  the  preform  using  a 
oxyhydrogen  torch. 

Prior  to  the  start  of  the  draw  process  all  equipment  will  be 
checked  and  calibrated  to  assure  performance  within  specifications. 

4.3  Monitoring  Procedures  for  Materials 

Materials  used  in  the  preform  deposition  area  and  fiber  draw 
area  are  obtained  to  exact  specifications  and  many  are  accomp¬ 
anied  by  a  vendor's  certification  of  analysis.  All  materials 
that  can  be  checked  by  using  simple  laboratory  practices  will  be 


inspected  for  compliance  with  specifications.  All  materials 
used  will  be  stored  and  handled  in  such  a  manner  as  to  not 
alter  their  characteristics. 

4.4  Test  Methods 

Some  basic  tests  are  performed  to  check  the  preform  quality- 
prior  to  draw.  Any  preforms  exhibiting  properties  other  than 
those  specified  within  a  degree  of  error  are  not  drawn  into 
fibers.  These  tests  are  outlined  in  paragraphs  4.4.1  through 
4.4.3. 

4.4.1  Visual  Inspection 

The  preform  is  thoroughly  inspected  for  metal  impurities  on 
the  preform  surface  and  for  bubbles  included  within  the  preform 
glass;  observation  of  such  defects  requires  the  rejection  of  that 
preform  portion. 

4.4.2  Preform  Diameter  Ratios 

A  cane  sample,  pulled  from  the  collapsed  preform,  is  inspected 
under  a  microscope.  A  photograph  of  the  cross-section  of  this 
cane  is  made  and  the  diameters  are  measured  to  check  for 
correct  ratios. 


4.4.3  Mechanical  Test 

All  fibers  are  mechanically  tested  to  assure  that  they  meet 
contract  specifications.  Static  fatigue,  dynamic  and  bend 
testing,  are  performed.  The  first  two  tests  are  described  in 
the  Technical  Proposal  No  R77-043  submitted  by  ITT  Electro- 
Optical  Products  Division  on  November  1,  1977.  The  bend  test 
is  described  in  the  Eighth  Bimonthly  Progress  Report  submitted 
by  ITT  Electro-Optical  Products  Division  on  April  12,  1980 
under  this  contract. 


E-33 


S.O  NOTES 


S.l  Intended  Use 

The  intended  use  of  this  process  is  to  produce  a  durable 
high  strength,  long  length  optical  fiber.  The  durability  of 
this  fiber  is  enhanced  by  an  outer  compressive  layer  fabricated 
according  to  the  process  specification  described  herein. 


